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Abstract 


In this thesis, we study properties of many-body fermionic systems in a nonperturbative method 
of quantum field theories. Especially we will concentrate on thermodynamic properties of the 
BCS-BEC crossover and those of dipolar fermionic systems. These systems are recently studied 
from both the experimental and theoretical aspects, and thus we need to invent a theory to 
describe them quantitatively. For this purpose, we use the functional renormalization group 
(FRG) method and develop it in fermionic systems. It is characteristic to this study discussing 
superfluid phase transitions without introducing Hubbard-Stratonovich fields, and it enables us 
to study those systems in a less biased and nonperturbative way. As a preparation for studying 
the BCS-BEC crossover, we discuss the vacuum physics in detail, especially the atom-dimer 
scattering problem. We also derive some formulae to estimate the number density of fermionic 
atoms, and rederive the Thouless criterion using the Ward-Takahashi identity. Applying FRG 
to the BCS-BEC crossover, we have reproduced the NSR theory in a systematic analysis and 
discussed possible corrections to it. In dipolar fermionic systems, we established the Landau 
Fermi liquid theory with the dipole-dipole interaction using FRG and revealed that the system 
shows °P, superfluid at sufficiently low temperatures even in weakly coupling regions. 
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Chapter 1 


Introduction 


In this thesis, I will discuss many-body problems of non-relativistic fermionic theories using 
the functional renormalization group (FRG). Before going into the details, I would like to give 
brief introductions in order to motivate us to have interests in these topics and to study them 
via nonperturbative methods of quantum field theory (QFT). I hope that I can tell not only 
how interesting physics these topics contain but how wide range of physics shares the similar 
kinds of problems. 


About QFT 


QFT is now a standard tool for broad areas of physics, especially for condensed matter physics, 
nuclear physics, and particle physics. In the case of condensed matter physics, suppose that we 
would like to reveal thermodynamic properties of matter. With usual quantum and statistical 
mechanics we must solve the Schrodinger equation with too many variables. In many interesting 
cases, we must cease to expect that we can solve such problems in this way. In the case of 
particle physics, we also encounter the problem which may seem to be more fundamental. If 
we require the relativistic description of quantum mechanics, antiparticles should appear and 
thus processes with annihilation and with creation of particles appear from this fundamental 
requirement TAR Historically this fact posed some paradoxes when one interpreted relativistic 
wavefunction as a single-particle one (e.g. see the Klein paradox (31/4]). 

Surprisingly, answers for these questions posed in different contexts are the same: We should 
use the field description of quantum physics, i.e. QFT. In my understandings, we must rely 
on the renormalization group (RG) and its application to effective field theories in order to 
understand why the same method applies to such different situations. In QFT, every physical 
quantity is written in terms of fields, and a particle is regarded just as a cluster of a given energy 
and momentum. Only when the perturbation theory in QFT works well, the elementary field 
is composed of creation and annihilation operators of the corresponding particle, which itself is 
not required from any principles of QFT. Indeed, the ergodic or stable property of the vacuum, 
or the cluster decomposition property, allows us to introduce creation/annihilation operators 
for particles (12) 5, but they are not needed for the elementary ones. 

As long as low-energy physics is concerned, short-distance behavior is washed out and 
physics can be described with an effective theory related to the original one via RG transfor- 
mations (6). That is, we can expect universality when physical scales such as the correlation 
length are much larger than microscopic scales. At this stage, we already need not require the 
naive renormalizability with power counting, and the only important thing is that the effective 
couplings are controlled under coarse-graining with the aid of the RG flow. 
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If we use this idea to the get an effective field theory, we at first have to write down any 
possible terms allowed by the symmetry with given fields, and then calculate Feynman diagrams 
with a given cutoff. In this case, the bare Lagrangian does not take any simple forms and it 
contains many non-renormalizable terms. This was pointed out by S. Weinberg |7| and it was 
used to describe soft pions without current algebra. The methods and ideas of an effective field 
theory are now widely used. 

These experiences tell us an important lesson: Even if the high-energy theory may not be 
suitably described with QFT, QFT is still a powerful tool to describe low energy physics and 
provides a systematic expansion in terms of the energy under the name of RG. To see this, let 
us explain some classic examples of effective field theories taken from condensed matter physics 
and from particle physics. 

At first, let us consider the theory of superconductivity. In many metals, it is well described 
with the Bardeen-Cooper-Schrieffer (BCS) theory (3), that is, the electrons are weakly attracted 
via phonon exchange and it leads to the Cooper instability of the Fermi surface at sufficiently low 
temperatures. However, we must point out that the attractive interaction energy via phonon 
exchange is typically of the order of ImeV and that it is much weaker than the Coulomb 
repulsion which energy scale is about the order of 10eV. At this point, the notion of effective 
field theories is very important. After coarse graining of the full theory containing the bare 
Coulomb repulsion, we can obtain the Landau Fermi liquid theory as its low-energy effective 
theory. Its typical cutoff scale is given by the Debye temperature, and for such low energy 
processes, the Coulomb repulsion is screened and weaken. As a result, the phonon attraction and 
the effective Coulomb repulsion can compete. When the total interaction becomes attractive, 
the condensation of the Cooper pairs is formed. 

In particle physics, theory of the Fermi weak interaction is regarded as a low energy effective 
description of the Weinberg-Salam theory. Its typical cutoff scale is given, for example, by the 
W-boson mass my œ~ 80GeV. Above that scale, the Fermi weak interaction theory breaks 
down and many important properties such as the unitarity does not hold. As long as we only 
consider low energy phenomena such as (6-decay, then the Fermi four-point coupling gives a 
convenient description. 

Even for the standard model of particle physics, we believe that it should be a low-energy 
effective description of some UV complete theory. If we adopt some grand unification scenario, 
the standard model with the gauge group SU(3) x SU(2) x U(1)y is an effective theory with 
a cutoff scale about Agur œ 10!°GeV. Since our usual experiments search physics up to about 
1TeV, those two energy scales are well-separated. Therefore, as long as we believe that the 
perturbation theory works well in between those energy scales, the non-renormalizable couplings 
are suppressed with a factor 1TeV/Acgur œ 10713 « 1. It explains why we can describe particle 
physics using a renormalizable theory, and as a result we can conclude that the baryon number 
conservation, for example, is a good conservation law as a result of an accidental symmetry 
since it can be violated only through higher-dimensional operators. 


About FRG 


In this way, the notion of RG pioneered by K. G. Wilson is applicable widely, and the description 
with effective field theories provides a systematic expansion in terms of energy scales. In the 
context of gauge theories, the inconsistency between the gauge symmetry and separation of 
energy scale poses a serious problem, and the nonperturbative and practical realization of the 
Wilsonian RG is still open in this area. On the other hand, in condensed matter physics such 


problem associated with the non-linear symmetry transformation is often absent. The FRG is 
a realization of the Wilsonian RG in a suitable form for practical computations. 

In 1984, J. Polchinski gave a version of the FRG in order to give a simple proof of perturba- 
tive renormalizability of the ¢j-theory (9). He separated the propagators into low-energy modes 
and high-energy modes, and derived a differential equation of RG transformations by gradually 
integrating out high-energy modes. At each step, the modes lying only in the infinitesimal 
shell of Euclidean momenta are integrated out, and then the differential equation consists only 
of tree and 1-loop diagrams. Of course, the diagrammatic expansion of correlation functions 
solves the RG differential equation, but those Feynman diagrams has forest structures [10]. 
Therefore, the subtraction of their subdivergences is necessary for the renormalization, and the 
proof of Dyson’s counting rule becomes awful, because we have to take care about symmetry 
factors for correct cancellations of each subdivergence [10]. If we use the RG differential equa- 
tion given by J. Polchinski, we are free from those subdivergences and the proof of perturbative 
renormalizability is given with a simple mathematical induction. 

After about a decade, a suitable form of FRG for practical applications was established. 
This is called the Wetterich formalism, and it gives RG transformations for the 1PI effective 
action [H3]. If the tree structure appears in the vertex function, some drastic momentum 
dependence can appear, because if the sum of the external momenta exceeds the cutoff it 
suddenly flows into the internal tree leg [14]. Since the resultant vertex function sticks to the 
1PI structure in the Wetterich formalism, such drastic change would not appear and it is easier 
to treat in practical applications. At the end of the RG flow, it indeed gives the 1PI effective 
action of the original theory. 

In a bosonic system, local potential approximation (LPA), which is the lowest order ap- 
proximation in the derivative expansion, already gives a powerful tool to study nonperturbative 
physics. In this case, the notion of the optimized regulator is well established and also many 
other techniques are developed (15/16). On the other hand, in the case of fermionic field theories 
there are many problems associated with the Fermi surface even in the weak coupling region. 
Technically, due to the anticommuting behavior of the fermionic fields the naive application 
of LPA is meaningless and it makes the nonperturbative problem of fermionic theories more 
difficult [17]. In the case of strong coupling systems, the identification of correct low-energy 
degrees of freedom is a difficult problem and in most cases we need to put some ansatz to treat 
such systems. This is true also for the bosonic theories. Indeed, the formation of a bound state, 
for example, changes the scaling behavior drastically and the counting rule of the expansion 
will be totally changed from that in the weak coupling region. 


Problems tackled in this thesis 


In this thesis we will attack problems of the Bardeen-Cooper-Schrieffer to Bose-Einstein conden- 
sate (BCS-BEC) crossover and of dipolar Fermi systems. Both of them now become important 
topics in the field of ultracold atoms. I will give an introductory review for each topic at the 
beginning of the appropriate chapter. Here I would like to discuss why they are interesting and 
also why they are difficult from a theoretical point of view. 

At first, let us consider the BCS-BEC crossover. Imagine that we have a many-body two- 
component fermionic system with a short-range attractive interaction. If the coupling is very 
weak, the system shows superfluidity described by the BCS theory. On the other hand, if 
the attraction is very strong, fermionic atoms form a bound state (often called a dimer) and 
the ground state is a superfluid with the Bose-Einstein condensation (BEC) of those bosonic 
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dimers. If we change the coupling in between those two limits, the ground state changes 
smoothly without phase transitions and this is called the BCS-BEC crossover. This system is 
now experimentally accessible thanks to developments of ultracold atomic experiments. 

The BCS-BEC crossover gives a unified description of two different kinds of superfluidity, 
which provides important theoretical challenges. If we could give a clear and quantitatively good 
prescription for this problem, we would acquire a new insight for the superfluid phase transition 
especially for the strongly coupled systems. We can expect that we will learn very important 
lessons which may open ways for understanding of high-temperature superconductivity and for 
other strongly coupled systems associated with phase transitions such as the chiral symmetry 
breaking in QCD. However, the quantitative description of the BCS-BEC crossover remains to 
be intractable, although the system is defined with a very simple classical action. 

We will analyze the BCS-BEC crossover applying FRG within a purely fermionic framework. 
That is, we do not introduce auxiliary bosonic fields using the Hubbard-Stratonovich transfor- 
mation in order to describe Cooper pairs in the BCS regime and bosonic molecules in 
the BEC regime. Of course, the Hubbard-Stratonovich transformation provides a powerful tool 
to study nonperturbative physics when we already know the correct channel of instabilities, 
since it describes the condensates already in the mean field approximation. However, it always 
contains some ambiguities associated with the Fierz transformation so that different 
Hubbard-Stratonovich transformations give different results within the mean field approxima- 
tion. and for more difficult systems the correct way to introduce Hubbard-Stratonovich fields 
is often nontrivial. We shall consider the BCS-BEC crossover with a less biased method and 
see how instabilities emerge in the purely fermionic language using FRG. 

Next we will consider a dipolar fermionic system. Compared with the BCS-BEC crossover, 
this system contains many open problems, and due to the complexities of interactions, we 
can expect the occurrence of many unprecedented phenomena. The dipole-dipole interaction 
is anisotropic and long-ranged interaction, so scattering processes with higher harmonics will 
become dominant even if we consider low energy phenomena. 

Furthermore, what is interesting is that those systems may be helpful to understand prop- 
erties of high density nuclear matter. In a normal nuclear matter, the interparticle distance is 
very long and anisotropic natures of nucleon-nucleon interactions does not affect many-body 
properties. However, nucleon-nucleon interactions contain anisotropy coming from the LS in- 
teraction, tensor force of pion exchanges, and so on. In high density nuclear matters, these 
interactions play an important role and the mean field calculation predict interesting physics 
such as neutron °P, superfluid, meson-condensation, etc. Those phenomena are considered to 
occur inside neutron stars and to affect cooling mechanism. It is still difficult to observe such 
phenomena but some counterparts in cold atoms may give some insights in a near future. 


Composition of the thesis 


Before closing the introduction, let us comment on composition of this thesis. 

In chap }2| we will review the method of FRG. At first, we introduce the language of field 
theories using functional integrations, and we derive the Wetterich formalism, which will be 
used in this thesis for practical computations of 1PI effective actions. To get a physical insight 
of FRG, we compare it with another version of FRG proposed by J. Polchinski. 

In chap}3} we introduce notations for fermionic systems and calculate some analytic expres- 
sions of the RG flow, which will be frequently used later. We also review the Landau Fermi 
liquid theory from the perspective of FRG, and understand it as an IR effective theory of many 


body fermionic system. The idea established here will play an important role, especially in the 
analysis of the BCS regime in the BCS-BEC crossover and in the analysis of dipolar fermions. 
We also marshal problems to be considered in following chapters from lessons learned through 
the RG flow of Landau Fermi liquid theory. 

In chap {| we attack the BCS-BEC crossover using FRG with purely fermionic description. 
Starting from scattering problem of the model in the vacuum, we consider superfluid phase 
transition using the RG flow in normal phases. We derive the number equation and the Thouless 
criterion from the RG point of view, and consider possible corrections to the Nozieres-Schmitt- 
Rink theory in a systematic way. 

In chap {5} we discuss possible instabilities of dipolar fermionic systems using the scaling 
ansatz of the Landau Fermi liquid theory. We will find that the ansatz is justified in weak 
coupling regions and that 3P, superfluid is predicted from our calculations. 

We summarize results in chap [6] and discuss some future perspectives there. 

The rest of the thesis is composed of appendices. In Appendix[A] we denote some computa- 
tional techniques to derive diagrammatic expressions of the RG flow. We also discuss the way 
to compute composite operators there. In Appendix |B} we review the proof of perturbative 
renormalizability using FRG, which was originally proposed by J. Polchinski. It will make 
the connection between renormalizability of quantum field theories and fixed points clear. In 
Appendix |C] we review the theory of angular momenta in order to fix its notations and conven- 
tions. It is used in calculations of the dipole-dipole interaction and we also list some convenient 
formulae. In Appendices [D] we show some formula used in this thesis. In Appendix [D] we 
calculate some properties of the effective four point coupling in the vacuum, and they are used 
for the analysis in chap|4] In Appendix |E} we can see details of calculations of the interaction 
matrix of the dipole-dipole interaction which is used in chap [5] 
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Chapter 2 


Functional Renormalization Group 


In this chapter, we will review and explain the functional renormalization group (FRG), which 
will be applied to many-body fermionic systems in latter chapters. This is one of the concrete 
realizations of the renormalization group introduced by K. G. Wilson (6), which was originally 
invented for clear description of critical phenomena. After about a decade, J. Polchinski (9| gave 
another exact formulation of the renormalization group in simplifying the proof of perturbative 
renormalizability of the scalar øf theory in 4-dimensional spacetime. The starting point of 
FRG is based on exact relations among generating functionals of Green functions, and there 
are several and equivalent versions of FRG depending on the choice of types of Green functions. 
In this chapter, we will see two types of flow equations of FRG: the Wetterich formalism for 
one-particle irreducible (1PI) vertices (sec[2.2), and Polchinski’s equation for the Wilsonian 
effective action (sec 2.3). In latter chapters we will only use the Wetterich formalism of FRG, 
however seeing the relationship with other formulations will give us deeper insights for FRG. 

In sec [2.4] we will give a discussion on the choice of regulators, which is an unphysical 
degree of freedom introduced in FRG in order to divide scales of the physics. In principle, the 
results do not depend on the choice of regulators since we remove it at the end of calculations. 
However, in practice some approximations are necessary in solving the physical systems, and 
those approximations can break down the independence of this unphysical degree of freedom 
and cause a delicate problem for the choice of regulators. 

In sec [2.5] we discuss general properties of exact solutions of the RG flow equation. The 
notion of fixed points becomes important for this purpose, and the rescaling operation is intro- 
duced as a convenient tool. 


2.1 Preliminaries 


At first, we give the definition of notations and review basic properties of generating functionals. 
Let us denote the elementary fields by n, where the label n represents the index (x, a) covering 
the space-time coordinate x and the flavor index a. For a functional F'[¢] of elementary fields 
on, F” represents the left functional derivative: 


ÔL 
F” |o] = —F 2.1.1 
(ol = FEF (2.1.1) 
where the left and right derivatives in terms of fields ¢, are defined by 
iu Flé] _ bk (6) 


Sdn. (2.1.2) 


Fo + 66] — FIO] = 86,48 = 


7 
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For two fields A, and B™, we use Einstein’s convention for the integral 
A,B" = f dx Y` Aala) B° (2). (2.1.3) 


Let S|ġ] be an action which is a bosonic local functional of the elementary fields ¢,, and 
for each n we introduce the source J” which has the same statistics with that of n. When 
we need to distinguish the statistics of fields, œ? and ¢; represent bosonic and fermionic fields, 
respectively, and corresponding sources are written as JE and J}. Therefore, J"¢, = J*¢?+Ji¢} 
is a bosonic quantity and we define the generating functional W|.J] of connected Green functions 
as follows: 


Definition 2.1.1 (Schwinger functional). We define the generating functional for connected 
Green functions W|J] by 


exp (W|[J]) = [Pee (—S|¢] + Jon) - (2.1.4) 


W [J] is called the Schwinger functional, and its moments are called Schwinger functions. Here 
we are not interested in the normalization of the right hand side of Eq.(2.1.4), because the 
change of the normalization only affects to add the field independent term in WJJ]. 


Connected Green functions are given by 


ÔL OL 


SJA nt 


(Dy ite = W icc lel | WJJ]. (2.1.5) 
Now we would like to introduce the one-particle-irreducible (1PI) effective action I'[y] as the 
Legendre transformation of W|J]. However, in order that the definition of the 1PI effective 
action makes sense, the convexity of the generating functional W|.J] is necessary and sufficient, 
and thus we need to check it at first. 

Let us make a more precise and concrete statement. Let us define a functional y,,|J] by 


nl J] = Wald] = (nel. (2.1.6) 


According to a general theory of the superspace, the existence of the inverse function J"[y] is 


ensured if the matrix of derivatives zgr Prl] is invertible. Since 
OL 
ggm Prl] = Winn J] = (Omon)cl J], iy) 


the existence of the inverse function J”[y] is nothing but the existence of the inverse Green 
functions with any external fields J. Since any functions are expanded in a formal power series 
in terms of Grassmannian variables, what we need to do is to check existence of the inverse of 
propagators (@mPn)c|Jo], where J, is a bosonic external source. 

Now let us see details about the matrix (dmon)-|Jo]. At first, we can immediately notice 
that the matrix is block diagonal in the sense that (¢°¢1,). = (¢1¢°,) = 0 since Jı = 0, and 
then 


(GnOm)elJo] = ( ‘Pata el Fe (gi Pi Ld ) (2.1.8) 


2.1. Preliminaries 9 


Clearly, the boson sector (¢°¢2,). is a symmetric matrix and the fermion sector (¢)¢},), is a 
skew-symmetric matrix. 

In the case of a theory only with bosons, we can prove the convexity of the effective action 
with a simple calculation. 


Theorem 2.1.1 (Convexity of the Schwinger functional). Let us consider a theory only 
with bosons |'| described by the classical action S|@°|. Then the eigenvalues of the matrix 
(Pp )elJo] are all positive, which may be infinite. 


Proof. Let p” be a real vector, then one can easily find that 


fl DODE GG ee eee 


PP” (PrPmdelFo| J Doe DGe-W 1+ Ig e—MP Tt Jo$” ; 


(2.1.9) 


and the right hand side of the equation is clearly positive. Since the matrix (%2, [Jo] is 
symmetric, we find that all the eigenvalues are positive. 


Theorem ensures that W[J] is convex and then its Legendre transformation can be 
performed. For the case of theories with fermions, we here simply assume |f| that the inverse 
function of yn|J] exists globally, which is denoted by J”[p]. Then the 1PI effective action is 
defined: 


Definition 2.1.2 (The 1PI effective action). The 1PI effective action I'|y] is defined as 


le] = lelen = WL yl), (2.1.10) 
where J”|y] is the inverse function of Wn[J] = pn. 


As a consequence of properties of the Legendre transformation, the following can be con- 
firmed immediately: 


Theorem 2.1.2. We can get W[.J| as the inverse Legendre transformation of The]: 


WIJ] = J’ pal] = Tiel, (2.1.11) 


5 [ylJ]] = J”. Therefore, we can find the relations 
Pn 


SWI] _ SLPmlI] brde0 ly] _ ôL” [y] 
ôJeôJ™ ÔT” ° EpPnôpm Opn 


Therefore, using Eqs. and the Leibniz rule, we get the two relations 


(2.1.12) 


6?7W|J] LRT [y] _ gm ôLôRT [p] 05 WJ] sr. (2.1.13) 


OJ T! 6y16Ym d Ôpnðpı 6J'6J™ me 


'The Yang-Mills theory is not in the scope of this theorem, since it essentially contains F.P. ghosts. Since 
the Faddeev-Popov determinant can change its sign, the proof given here does not hold. Therefore, in the gauge 
theories, the more careful treatment of the path integral becomes important and this would be related to the 
Gribov problem. 

2To my best knowledge, the proof for general cases does not exit. Since the positivity of two-point Green 
functions is related to well-defined field propagation, the Osterwalder-Schrader reflection positivity may 
solve this problem. 
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Therefore, the quantity 
Pm e = — = Tp] (2.1.14) 


is the inverse of full propagators (éndm)- with the source J[y]}. 

In the perturbative expansion, I'|y] is the generating functional of 1PI vertices. This fact 
can be seen via the tree expansion of connected Feynman diagrams with full propagators and 
1PI vertices. It can be proven in a direct way, however we here prove it following the discussion 
in chap.16 of Weinberg’s textbook [2]. 


Proof. Let us define the generating functional of connected Green functions WIJ] with 1PI 
vertices by 


oxp 5 (Wi) = [re exp - (—I[y] + J” pn). (2.1.15) 


Here we introduce the loop-counting parameter h, and expand the generating functional WIJ | 
in the power series of A: 


WIJ] = 5 WJJ. (2.1.16) 


Now, let us take the limit A — 0, then the steepest decent method applied to the path integral 
in Eq.(2.1.15) becomes precise. Therefore, we find that 


WOL = J?vn[J] — Tle], (2.1.17) 
Ort 
OPn _ 
W[J] = W[J]. Connected Green functions are given by the tree-level expansion with full 


propagators and 1PI vertices. O 


with lp[J]] = J”. This is the inverse Legendre transformation (2.1.11) so we find that 


The following theorem turns out to be useful. 


Theorem 2.1.3. Let us consider an infinitesimal change AW|J] in the generating functional. 
Then the corresponding infinitesimal change of the 1PI effective action AT |y] satisfies 


AT vy] = —AW [J y]]. (2.1.18) 
Proof. We define the modified external field J’[y] by 


ôL(W + AW) 
oJ” 


Then the Legendre transformation is given by 
(P+ AP)[y] = J” [pln — (W + AW) Il]. 
Then, setting AJ = J’|y] — J[y], we get 


LW 
ðJ” 


AT[y] = AJ” pn — (ar 


m + awie) = -AW ][J[ġ]], 


which is the desired result. J 
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2.2 Wetterich formalism 


Our purpose in this section is to derive the flow equation for the one-particle-irreducible (1PI) 
effective action 13|, which is called Wetterich’s formulation of FRG. Here we use the path 
integral formalism developed in sec[2.1]to describe Euclidean QFT. 


2.2.1 Derivation of the Wetterich equation 


Let S|] be a classical action. Now let us consider the action with the infrared regulator Rx: 


1 1 
Slo] + 50° Reo = SO] + Ek” bnm. Ga 


Here the label k of the infrared regulator R; describes scale of infrared regularization. We 
assume that the regulator Rọ be bosonic, so bosons and fermions do not mix through this 
term. In the boson sector Rg” = RẸ”, and in the fermion sector R?” = —R". 

Let us define the Schwinger functional Wẹ of the action S with the regulator Rp by 


exp (W;[J )) = | Does (-s l$] — sO Rk- @4+ J: 6). (2.2.2) 


In the limit Rẹ — 0, we can find that W;|.J] converges to the Schwinger functional W|J] of the 
original theory. 


Definition 2.2.1 (The flowing 1PI effective action). We define the flowing 1PI effective 
action [;,[y| as the Legendre transform of W;,{J] in the sense that 


1 
Dele] + eRe = lelen — Wille, (2.2.3) 
where J|y] is determined by the condition 
ôr Wk 
=| = Gp. (2.2.4) 
oF” lat 


In the definition of the 1PI effective action , we have subtracted the trivial contribu- 
tion for the two point vertex function from the ane In the limit of removing the regulator 
Rp — 0, this trivial contribution vanishes and W,|J] > W|J], and then I’,,[y] becomes the 1PI 
effective action T [p] of the classical action S[¢] in this limit. Furthermore, it should be noticed 
that the quantity [',,[y~] has a nice property also in the limit Ry, — oo: 

Theorem 2.2.1. Let us consider the limit Ry —> co. Then the 1PI effective action Tilẹ] 
converges to the classical action S|]. 

OrT 

Proof. Notice that the inverse function J”[y] is given as J"[y] = A 
Pn 


definitions and , we can find that 
1 
exp (lg) = | Desexp ( -S19 ~ FORKA + JMIeI(On = vn) + FeR ) 


dr 1 
= | Doexp (-si6+ 41+ eleg, — Zero). 
Pn 


In the limit Rẹ — oo, the saddle point approximation becomes accurate and then Tyly] > 
Sly]. 


+ YmR;,”. Using 
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Now let us derive the Wetterich equation which is the flow equation for the 1PI effective 
action [,[y]. Let us consider infinitesimal change dR, of the regulator Rg, then it induces 
change of the Schwinger functional W;|.J] given by 


AR ARIE fo (-333 oe z) em (-sia — “0 acts: 6) 


20S od 
16; br 


Therefore, we can find the expression for infinitesimal change dW;,,[{J] in terms of the Schwinger 
functional W;|[.J]: 


(2.2.6) 


1 2 
SWJ] = —SdRR" (on Ma ) 


Je Jm ` §JrôJm 
Theorem 2.2.2. Under infinitesimal change ôR, of the regulator, the 1PI effective action Tyle] 
changes by 
1 —1 
Prle] = 581 br, Ge [ye] + Rx) | l (2.2.7) 


The super-trace represents the ordinary trace of the boson sector subtracted by the ordinary 
trace of the fermion sector. Here we have defined that 


T 
p m ] _ ntl (2.2.8) 


Proof. Using the property (2.1.18) of the Legendre transformation, we have 
1 nm 


Then Eq. (2.2.4) and Eq.(2.2.6) give 


ô? Wi J] 
6J") J™ 


1 1 
Telo] = SORE” = ŻSTr oR (EY lel] + Ra). 


2 


Jie] 


Here, we have used Eq.(2.1.13) and the fact that the factor (—1) appears in the fermion loop 
= 


due to statistics. J 


We express this equation diagrammatically as 
Ô Rk 


m=i Cy , (2.2.9) 


TË + Ry) 
where the shaded blob represents the two-point vertex for infinitesimal change of the IR 
regulator ôR, and the line represents the field dependent propagator (ro ly] + Ry)7*. If we 
assume that we can expand the 1PI effective action Ffy] in terms of the field variables p, we 
can obtain the coupled ordinary differential equation for the vertex functions. 
From now on, we have discussed change of the 1PI effective action [’,[y] or the Schwinger 


functional W;,|J] under general infinitesimal change of the regulator Ry. Now, let us consider 
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the change of the regulator parametrized by infrared scale k € [0, A,], where A, represents the 
ultraviolet (UV) cutoff |*| of the theory described by the classical action S|é]. This produces 
the l-parameter flow of the 1PI effective action I,|y]. In order for desirable properties of the 
flow, we require the following properties of the regulator: In the Fourier space, let us write the 
regulating term as 


where f= f d’q 
i (2r)? 


1. Rz(q) > 0 in the limit q?/k? > 0. (This ensures that R,behaves as an IR regulator.) 


1 1 
a = g | aR Dbaa (2.2.10) 


and Qag = | Pesos, then 


2. Rk(q) — 0 in the limit q?/k? —> œo. (This avoids a new UV renormalization.) 
Pl. R, is an L?-function, and satisfies R, > 0 as k? > 0. 
3. Ry > œ in the limit k > A,. (The theory becomes classical at the UV scale.) 


Let us see details of these requirements. The property [I| guarantees an infrared regularization: 
as an example, R;(q) ~ k? in the region q? < k? for the boson sector satisfies the property [1] 
and it generates a mass k? in boson propagators. 

The property |2| is required to avoid introducing a new UV renormalization. For example 
R,(q) = k? (1 — k2/A2) ~ satisfies the properties |1| and |3| however it violates the property 
and it causes new ultraviolet divergences to the theory because it behaves a new mass term. 
The property [2] ensures the existence of the well-defined generating functional W;|J] and 1PI 
effective action I|] as long as the original theory, i.e., W[J] or I'[y], is well-defined. Therefore, 
the rapidity of the convergence R;(q) — 0 as q?/k? — 0 is important, and we can make a more 
refined statement P}. 

The requirement R — 0 as k? — 0 in the property [2] or [2} is also important, because [;,{y] 
becomes the ordinary 1PI effective action I'|y] at the endpoint of the flow k = 0. The property 
[3] ensures that the flow starts from the classical action P,—,{y] = S|y]. Thus, in principle, we 
can evaluate the effective action I'[y] from the classical action S|] by solving the flow for rfg] 
parametrized by k € [0, A]. 

As a corollary of Theorem [2.2.2| we can derive the flow equation for T|]: 


Theorem 2.2.3 (The Wetterich equation). The flowing 1PI effective action Tpi] satisfies 
the differential equation, called the Wetterich equation: 


(2.2.11) 


2 


1 
oT rle] = 3 STr 


oRk 
Rk 


rØ 


where ð; = kO,. The flow starts from the classical action: Tk-aly| = S|p] and reaches the 1PI 
effective action of the full theory: Ty=o|y] = Ty]. 


Let us assume that we can expand the 1PI effective action Tyf] in terms of field variables 
y and then we can draw Feynman diagrams for the Wetterich equation (2.2.11). For simplicity, 


3For the renormalizable theory, the ultraviolet cutoff scale A, is taken to be infinite. 
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we assume that the theory contains only bosonic fields, and that the Zo-symmetry y => —y 
holds. These assumptions allow us to expand the 1PI effective action as 


l — nm 1 N1,125-++5 N2i—1,N2i 
Pely] = 5¢n(Go* — Dr) pm + ` one RM a. bine Ou, (2219) 


i>2 


Substituting (2.2.12) into the Wetterich equation (2.2.11), we can find its diagrammatic rep- 
resentations by comparing coefficients (See Appendix |A| for the detailed derivation). For 1PI 
vertices with some small number of legs, we obtain that 


(Al x2 = ) — Xk = pe À (2.2.13) 
Tras tar + SOW l (2.2.14) 


and so on. In these graphs, square vertices represent the negative of 1PI vertices —T 4 9;, internal 
lines represent the full propagator (GZ! — Xp + Rp) `}, and the shaded blob represents the two 
point vertex Rp. The negative signature of square vertices reflects the negative signature in 
front of the action in the Boltzmann factor exp(—S|¢]) of the path integral measure. 


2.2.2 The flow of correlation functions 


In terms of FRG, it is quite natural to define composite operators as infinitesimal changes of the 
Wilsonian effective action (241/25). In this review, we have focused on the Schwinger functional 
W[J] and the 1PI effective action ['|y] instead of the Wilsonian effective action, and thus we 
can define the flow for correlation functions Jin a systematic way (26). 

Let I[J,@] be a bosonic functional of fields ¢ and sources J, then the expectation value 
I[J] = (I[J, @]) is defined as 


IJ] = ge WUE Ç z] eW U] = J Dold, o] exp (— S|] Ey 4 o) 


fDoexp (—S[¢] + J - e) (2.2.15) 


This is the quantity defined only for the full theory, i.e., at the end of the flow. In order to 
discuss the flow also for the functional J|.J], we extend this as 


I[J] = eV exp (i2 - Rg: z) (1 [4 Z] omit) (2.2.16) 


One should notice that the endpoint of the flow is clearly given by I,2o[J] = I[J] and that 
I[J] = 0 means that I[J] = 0 at any k € [0, A]. This fact will also be derived in another way 
later. 

Now we can readily derive the flow equation. Taking the derivative of Eq.(2.2.16) with 
respect to k directly will give the flow equation, however we here take another approach. 


4If we assume that I[J,¢] is independent of J, I,[¢] is a composite operator, and then this is a natural 
extension of the flow of composite operators. 
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Theorem 2.2.4 (The flow of the functional J[,|.J]). Let us consider infinitesimal changes 
OR; of the regulator, then the corresponding infinitesimal change of the functional I,|J] is given 
by 

(2.2.17) 


SI J] = oR" (oe Orta J] | 1 Tl] ) 


ol” osm 2 JrðJ™ 
Therefore, the right hand side is linear in I| J| and it implies that the solution of the flow is 
unique P} Especially, if Ip-o|J] = 0 then I,[J] = 0 for any k. 


Proof. Let us rewrite Eq.(2.2.16)) as 


exp (W;[J] + elx|J]) = exp (-3 Or OL 


d Bis z) exp (W[J] + el[J]), 


which is true up to the linear order of e. Here € can be regarded as an infinitesimal real 
parameter. Therefore, the flow equation (2.2.6) is true for W;[J]+e1,[J] up to the order O(e). 


The flow equation (2.2.6) for W;,[.J] is also true, so we find Eq. (2.2.17). 


In the diagrammatic interpretation, the first term of Eq.(2.2.17) is coming from the tree 
structure, and thus even when the bare quantity J|J,¢] has a 1PI structure the quantity [;[J] 
loses 1PI property. Next, we consider the correlation function J;,[y] in the 1PI formalism. Let 
us define 


Lly] = LJ (ell, (2.2.18) 


where J[y] is the inverse function of Wy n|J] = Yn. We would like to derive the flow equation 
for I; k|J]. As we have discussed, the flow equation for [;,{J] can be derived by recognizing the 
functional ¢J;,[J] as an infinitesimal change of the Schwinger functional. Since rafy] and Wz[J] 
is related by the Legendre transformation, this infinitesimal change induces the infinitesimal 
change to the 1PI effective action so that it becomes Ty[y] — el, [y]. 


Theorem 2.2.5 (The flow of correlation functions I,[¢])- The flow equation for bosonic 
correlation functions Ip] is 


= 1 w 1 
Tyly] = —=STr IP ig] OR. (2.2.19) 
2 PP w+ Re TP +R 
=(2)nm LÔRİ; 
Here we have defined that 7?) [yp] = dr6ntale) 
0YnOYm 


Proof. Since Taly] — el,[y] is the Legendre transformation of W;,[.J] + elg[J], it satisfies the 
Wetterich equation: 


1 
Ty] — ef lo] + Ri 


a, (Tale -— eill) = SST ƏR (2.2.20) 


By comparing linear terms in £ of the both sides, we reach the desired result. 


5When expanding the both sides of the equation in the proof with respect to the parameter £, we have to 
assume the absence of phase transitions associated with the perturbation ¢J;,[¢]. This stability of the vacuum is 
related to the cluster-decomposition property (see chap. 4 of |1| and chap. 19 of [2)). For solving degeneracies 
of the vacua, we need to take care about boundary conditions of the infinite volume limit or, equivalently, to 
add the sufficiently many terms to the classical action (see chap. 5 and 16 of B). 
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The flow equation (2.2.19) for I; plp] has exact 1-loop property. Thus, especially when bare 
quantities Ip. |p] are 1PI correlation functions, the 1PI structure of correlation functions I; [y] 
holds for arbitrary k € [0, A]. Diagrammatically, we can express this as 


ae 1 = 
rle] = -3 OR, TP Ty), (2.2.21) 


where the shaded blob ay he the scale derivative of the regulator Rp, the empty blob 
repre enia the vertex for E |p], and the line represents the field dependent propagator (I; Ne] + 
Ri, 

From now on, we have restricted ourselves to the case where bare functionals I[J, ¢] are 
bosonic in order to deal those as infinitesimal changes of the action. Now, let us allow for the 
case where functionals I[J, ġ] are fermionic quantities. This case, however, can be reduced to 
the bosonic case. Let 7 be some fermionic variable, then 7J{.J,¢] becomes a bosonic quantity. 
Therefore, all the discussions given above are valid for 7J{J, | when I[J, ¢] is fermionic. 

Recall that the regulator Rg” is a bosonic quantity, so in Eq. (2.2.17) the fermionic parameter 
7 can be factored out since the labels m and n in the regulator can only run over the fields 
with the same statistics: 


Wl nhl] , 18l] 
I Z nm 
nô lJ] OR, ( 5 J” SJ 2 JJM 
Me (CEJ SJ 28S" 


Therefore, Eq.(2.2.17) is valid also for the fermionic quantity. 
Next, let us check for Eq.(2.2.19) in the case of fermionic functionals. Denoting Ginn|y] = 
We,mn[J[y]], we get 


~ 1 ÔL ÔR + 
noatslel = -30RE" Cmte ($4 E nide] ) Grall 
l 


Therefore, the extra minus sign appears only when m and n are fermionic labels in factoring 
out 7 from the left. Therefore, we get the flow equation for the fermionic functional Iyọ]: 


zs / 1 1 ~ 
alo] = -3R G mi HP" (elGuale] = -3T Pra 
k k 


5 OR 


TP [el + Ry 
(222) 
In other words, since the vertex becomes fermionic, the extra minus sign for the fermionic loop 
does not appear, and thus the super-trace should be replaced by the standard trace operation. 
In the following, we sometimes denote the flow equation as Eq. (2.2.19) instead of Eq. (2.2.22) 
even for the fermionic functional T, for the simplicity of the notation. : that case, we should 
interpret the super-trace as the trace operation which gives the extra minus sign only when 
the fermionic loop is closed. In that interpretation, the fermionic loop is not closed in the flow 
equation of the fermionic vertex [, and therefore the super-trace in this definition coincides 
with the standard trace operation. 
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Dyson-Schwinger equation 


As an example of applications of Eq. (2.2.19), we will find a relation of FRG to the Dyson- 
Schwinger (DS) equation. The Dyson-Schwinger equation is the equation of motion 
in quantum field theories. In the path integral formalism we can understand that the DS 
equation is nothing but the identity, which represents the invariance of the path integral. The 
most standard example of the DS equation can be derived from the invariance of the flat path 
integral measure D¢ under translational shifts ¢ 4 ¢ + ed, where e is an infinitesimal real 
parameter and @’s are arbitrary. Then, we find that 


foo Gi — = o) exp (—S|¢] + J- ¢) =0, (2.2.23) 


which is nothing but the DS equation. 
Let us rewrite the DS equation by using a correlation function. 


Definition 2.2.2 (The DS operator). The Dyson-Schwinger (DS) operator Jpg is defined as 


ôL ÔRS | ôL 
Is | J, =| = J” - ails 2.2.24 
With the DS operator, the DS equation can be written as [fg [J, 24] exp (W[J]) = 0. 


Since we have constructed the correlation function which denotes the DS equation, we can 
apply the flow equation by extending it with Eq. (2.2.16). Before going into this topic, let us 
become more familiar with the DS equation. Since we have obtained the DS equation for full 
Green functions, we can easily rewrite it as the DS equation for connected Green functions: 


n [SRS [SWI] , on] N _ 
r- (=| T Z|) =o. (2.2.25) 


Following the notation given in Eq. (2.2.15), the left-hand side of Eq. (2.2.25) represents the 
normalized expectation value I},|J]. Let us rewrite Eq. (2.2.25) as the DS equation for 1PI 
vertices, which can be achieved by the substitution of J = J|y]. In this change of variables, 
the derivative with respect to J” can be written as 


ÔL 
Jig] Pm 


Ar 
ÔPm 


ÒL ÔLPm 


ðJ? êJ” 


= W nm[J [l] 


Let us denote Grm[y] = W nml], then we can rewrite the DS equation as 


dal [y] _ ORS 
0Pn Pn 


[e + Gly] - A i) =0. (2.2.26) 


This represents TZ [p] = 0. 
Now let us apply the flow equation to the functional Ips. Since the subscript becomes too 
messy if we simply use the notation given in Eq.(2.2.16), we define another notation as follows: 


1 dz ÔL 


IBs[J; Ry] = exp (—W,[J]) exp (-334 Rie: z) (1 Ç z] exp avtu))) 2. G77 
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As we have done in Eq.(2.2.18), we define the correlation function for the DS equation in the 
1PI formalism as - 
Ipsly; Re] = IbslJ ly]; Re, (2.2.28) 


where W;,[J[y]] = Yn. Since J-dependence of the functional I},[J, o] is very easy, we can 
explicitly calculate Ips[J; Ry] and Ips[y; Rx]. By the definition of Eq. (2.2.27), we get 


S Do (I = bmp” — FA) exp (—S[g] - 10- Reo + J-4) 


fee Rg] = 
dsl; Ral f Doexp (—S[d] — $6: Ri 6+ J-¢) 
ORS ô 
= Jp ie (£ eu + i) (2.2.29) 
Substituting J = J|] into Eq. (2.2.29), we obtain that 
oe dal kle] ORS Or 
Islo; = — ii mi=—|1). 2.2. 
DS ly; Rx] bn don p F Wh, Sy ( 30) 


Of course, in the limit k > 0, or Rg — 0, we get Tos lp, Rk > 0] > Tps |. One can find that 
Tos lp; Rg] is a very natural extension of Tosi] by comparing Eqs. (2.2.26) and (2.2.30). 

In the context_of FRG, the Dyson-Schwinger equation_can be regarded as the result of the 
flow equation for Ips|y, Rx], as pointed out in the review given by J. M. Pawlowski. In the 
classical limit k + A, or Ry — oo, we can find that Trake] > Sy] and that Wimnl[J[y]] > 0, 
and thus Tps lp, Rk] —> 0 from the definition of Eq. a Since the functional Tps[y; Ry] 
satisfies the flow equation (2.2.19), the existence and the uniqueness of the solution imply that 


Tos; Ri] = 0 (2.2.31) 
at any k € [0, A]. Especially, this means that the DS equation Ing lyp] = 0. 


2.2.3 Renormalization group flows 


So far, we have considered flows accompanying to infinitesimal changes of an IR scale k, which 
is a parameter of the regulator Rẹ. Since the scale parameter k does not appear in the full 
theory (i.e. in the limit k — 0), the Schwinger functional W|./] of the full theory is independent 
of k: kO,W [J] = 0. From this viewpoint, we can generalize situations in the following way [26]. 

Let us consider a theory with fundamental couplings g;. We would like to study the response 
of the theory under infinitesimal changes of some scale s. For example, we can put s = k, which 
is an IR scale of FRG, and we can also put s = upr, which is a renormalization scale of the full 
theory. Let us regard that the couplings g; and external sources J” can be s-dependent A We 
are interested in a differential operator 


o >. 0 ôl 
Dss dh i a S IPR 
"Os + Wg 0g; +m ga” ( ) 
which satisfies 
D,W[{J] = 0. (2.2.33) 


°For example, let the couplings g; be defined at some renormalization scale upg. Then, g; should depend on 
s if we set the scale s = wr. Furthermore, in this situation the external fields J should also depend on s if J’s 
couple to renormalized fields. 
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In other words, we consider situations in which we can cancel the change of the Schwinger 
functional under the suitable rescaling of the couplings g; and the fields J”. For simplicity, we 
restrict ourselves to the case where anomalous dimensions yz do not changes the statistics and 
are independent of external sources J. 

We would like to calculate the flow D,W;[J]. By the assumption (2.2.33), W[J] does not 
change under the transformation s +> (1+ €)s, gi œ gi + EV 9! gj, and J” œ J” + eyz J™. We 
should recall that 


exp (Wal) = exp ( -5 Re" EEE ) exp WUD. 


Under infinitesimal changes s ++ (1 + ¢)s, the regulator term changes as 


ON E E T Ar hes 

Rem Rm DR E “L — - 

gtk” Sym gym a>" + DSR (2 ewer (a magy 
1 nm nm n plm OL ÔL 


Let us denote as |(D, — y7)R,|"" = D, R}™ — 273” R”. Then, we can derive the flow equation 
of the Schwinger functional W} [J] by replacing OR, by (Ds — yz) Rx in Eq.(2.2.6): 


1 (a i ÔLWg| J] ôL W|I] 


D, e 
Weld] 2 \ SJJ” Je bm 


) (D= ys) Re] - (2.2.35) 


On the other hand, for the flow of correlation functions 7[J], the corresponding operator I[J, 24] 


may depend on s, which makes the equation a little bit complicated. Using the definition 
(2.2.16), we can easily find the RG flow of J;[J]: 

ôr Irl] brWalJ] Or til J] 
AJ J™ ôT»  ôJ™ 
Here, the left hand side of Eq. (2.2.36) is nothing but the linearization of Eq. (2.2.35) and the 
right hand side represents the contribution from the scale dependence of T[J, ¢]. 


Now, let us go back to the 1PI formalism. We would like to define anomalous dimensions 
Ye Of fields y so that 


Dahl) + 5 ( ) (Ds = IDR = (Date. (2236) 


DJ \e| = 0, (2.2.37) 
where 
s= 57. Gin. n>. 
Os Vo59 0g; Vem? Pm 
We can achieve this by imposing 7, = —y, under the assumption of Eq.(2.2.33). 


Proof. Let us calculate the left hand side of Eq. (2.2.37). Recall that T [y] is given by Eq.(2.1.10), 
then we can get 


(2.2.38) 


Dele] = (DI e)n + Yom” Pn — DWI]. 
The last term can be written in the following way: 


oW aJ” e] LW ~ (OW 8J” ôW SLJ! LW 
s- | T8 -mi eae | Pye a ST 
ðs ljg 3s OI J NOG; OG: OI" | i, Opm OF" | 5, 


] 
= (DW)[Jle]] — Vm” lelen + (DsI"[9]) en: 
Therefore, we have found that DTi] = (YI + Ye)mJ” Yn by using DsW [J] = 0. 


DW{J[y]] = 
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We can readily obtain the RG flow for I’,,[y~] under the total scale change of s: 


Dralo] = 581 (Dy +74) Rel (2.2.39) 


r?) ly] + Rk 


We can see that this is a natural extension of Eq. (2.2.11). By substituting J = J[p] into 
Eq. (2.2.36), we can find the RG flow of the correlation functions in the 1PI formalism: 


D, ile] + 581 [asle] (D + o)Ra] - Gell - 7] = (Da Mele, (2.2.40) 


where Gily] = (Ty fe] + Re)? 

The start point of the above discussion is the invariance of the Schwinger functional or the 
1PI effective action in the full theory. Therefore, we can straightforwardly lift Eqs. 
to equations including general variations of the regulator R;. We can achieve this simply by 
replacing 
6s(R, dR) x 


5 
—— as yo 2.2.41 
aR | T (R) (PA) 


The operator Dp represents the total derivative with respect to the regulator Rx. 


D, > DR= far eta 


2.3 Polchinski’s formulation of FRG 


In this section we shall explicitly give relationship between the Wetterich formalism and the 
another effective action formalism given by J. Polchinski (9), which is essentially same with the 
Wilsonian effective action. This is important not only because their relationship is interesting 
but also because it helps us to interpret the physical meaning of the regulator in Wetterich’s 
approach intuitively. 


2.3.1 Relation between the Wilsonian and 1PI effective actions 


Let S[¢] be a bare action of a continuum theory. Let us decompose it as 


=5 | % -paG P hpg + Sinto], (2.3.1) 


where [- [om Im)? Ọpa = y d‘r¢,(x)e””, and G,(p) is the free propagator. The idea of 


the Wilsonian ie action is very simple: by decomposing the fields into low-momentum 
parts and high-momentum parts and by integrating out high-momentum fields we would like 
to obtain the low-energy effective action which is called the Wilsonian effective action. This 
motivates us to write the free propagator in the following way: 


G.(p) = Go(p)K (p) + Go(p)[1 — K(p)], 


where K is some cut-off function which is monotonically decreasing and satisfies K(0) = 1, 
K(co) = 0. Its outline is shown in Fig 2.1} and the equation itself does not depend on details 
of the choice of cutoff functions. If we set K(p) = 1 — O(p?/k? — 1) with some cutoff scale k, 


2.3. Polchinski’s formulation of FRG 21 


the first term in the right-hand side represents the propagator of the fields with p? < k? and 
the second term is the propagator of the fields with p? > k?. Indeed we can easily show that 


Dé exp Psat ae Sine 9] (2.3.2) 
[Pex (-5 | = 


= | DoDo" exp (-4 [tha 5 | nage Rte- Sale + an) 


by performing the change of path integration variables. A to this fact, we define the 
interaction part of the Wilsonian effective in the following way: 


1 G5(p)? 
exp (—Sint,«[¢ p= fpe exp (-4 Ja P a — SintlO + an) i (2.3.3) 
then we get the relation 
Í em 
— = [Dee (5 [eno Es- Sml) ; (2.3.4) 


From Eq.(2 , we can find that the Wilsonian effective action is essentially the Schwinger 
functional ie a ear Let us see this explicitly. By shifting the path integral variable 


o + o" — ọ in Eq. (2 , we find that 
1 Cc. aß 
exp (Sm rlól) = f Do exp (-; J sa Rte Slt" 


+ fonat R Mi ‘the 5 | Oma Pr - na): (2.3.5) 


R = G7! (=; z i) (2.3.6) 


Let us define the regulator R by 


then we obtain that 


EE f Do” exp (-sio" — So" RO" + b(Gz" + Rio" — 50(G5 + n)a) 


= exp (mi . (G71 + R) — sO: (GZ! +R) 6) . (2.3.7) 
k Kalp) 
Kapy 
1.0 p/A 


Figure 2.1: Outlines of the cutoff function K4 and its derivative Ka (p) = AO, Ka (p) 


22 Chapter 2. Functional Renormalization Group 


We have shown that the Wilsonian effective action is essentially the generating functional for 
connected Green functions whose external legs are amputated. 

By using Eq. (2.3.6), we can find the regulator corresponding to the sharp cut-off in the 
Wilsonian effective action. In the sharp cut-off limit K(p) = 1 — O(p?/k? — 1) we obtain 


Rsnarp,k(P) = G3 (p) (soap - i) (2.3.8) 


In this way, the regulator R in the 1PI formalism can be related to the corresponding cut- 
off function of the Wilsonian effective action formalism. Furthermore, the relation Eq. (2.3.6) 
between R and K shows that the appropriately chosen regulators Rg in Wetterich’s formulation 
divide scales in physics and that it is closely related to coarse graining, or decimation of modes, 
in the Wilsonian renormalization group procedure. 

Finally, let us calculate the explicit relation between T|] and Sint,«[@]. Those two variables 
@ and ọ in the effective actions are related by 


_ SLWi[J] —OrWele- (G3 + R)] 


= (G+ R)! l (2.3.9) 
OF | 726(G51+R) og 

Since W, and Sint,« are related by Eq.(2.3.7), we obtain 

p=¢- (G7 +R) el (2.3.10) 
0) 
By using the definition (2.2.3), we can find that 
1 

Tile] +9: R- y= ¢- (G7! + R) -p — Wild: (G7! + R)]. (2.3.11) 


2 


Hence, we obtain that 
1 1 
Pale] = Sin,xl§] — 56° (Got +R) +o (G HR) p5p Ryo (2.3.12) 


= Sim xlo] + Ly -G7 - Zo — y)- (Go + R)-(¢-¢). (2.3.13) 


In the classical limit K — 1, the Wilsonian effective action Sint,x=1 is trivially the classical one 


Sint[@], which can be easily seen from Eq.(2.3.3) or Eq. (2.3.13). 


2.3.2 The Polchinski equation 


Now, let us derive the Polchinski equation (9), which is the flow equation for the Wilsonian 
effective action Sin x|]. Assume that the cut-off function K is a smooth function. Consider 
infinitesimal changes ôK, then they induce changes of the regulator R by 


ôK 
ooa 98 at KG. (Co! 
OR = Gy! am = (Ga +R) 6KG,-(Gz" +B) 
According to Eq. (2.3.7), we get 
ôL W;[J] 1 
Ô Sin = —dW,|J — ġo- R- -o -R-¢. 2.3.14 
szlo] el PEE 2 OF” hean T 3? ? ( ) 
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Here, the first term in the right hand side comes from the dependence of W;,[J] on R, or on K, 
with fixed arguments J, and the second term comes from the change of the variable 6J = ġ- ôR. 


Using Eq. (2.2.6), we find that 


1 ôr Wg ôL Wk ô Wk -o | 
= -§Rr™ — 2bn in 2.3.15 
5S [6] = 55 (SEE + Snape e E 838) 
Here, we should notice that 
| = by — (Got + Rj Sklo] (2.3.16) 
OJ” [j-e +R) Ôm 
Then the second derivative of W|.J] can be represented as 
ô? Wp [J] =i = ô? Z Sint Klel 
2 = (Go + R)im — (G3 + R) (G3 + R)m AA. 2.3.17 
I" I™ | 5 4.(G514R) i. aal mi 01,0 Ptp 


Up to a field independent term, we can obtain the following equation which represents the flow 
of Sint, K! 


Theorem 2.3.1 (Polchinski equation). The interaction part of the Wilsonian effective action 
Sint xlo] satisfies the RG differential equation 


1 ò Sint ô Sint ô? Sint 
aang JOK Gm ( i = “te oe aw 


under infinitesimal changes of K. This is called the Polchinski equation. 


(2.3.18) 


We can see that in the diagrammatic expansion the first term in the right hand side comes 
from the tree contribution and the second term comes from the 1-loop contribution. Thanks 
to this structure of the Polchinski equation, the proof of perturbative renormalizability for 
renormalizable theories becomes straightforward and does not require the analysis on detailed 
structures of each Feynman diagram (9. We need not care about subdivergences and then 
the complication due to Zimmermann forests is absent even in the discussion on all order 
perturbation theory. For details of the proof, see Appendix 

In order to obtain a diagrammatic expression for the Polchinski equation, let us assume 
that the interaction part of the Wilsonian effective action Sint, can be expanded in terms of 
field variables ¢: 


OO 


1 N1,N2,...,2m 
Sim rlo] = > Vien bm Png“ Pn: (2.3.19) 


m=1 


Just for simplicity we assume that fields ¢, are bosonic, so each interaction vertex Vg” is 
totally symmetric in terms of its indices nj,.. . The Polchinski equation (2.3.18) can be 
written as 


Tilya; > ea y% No(p)liq 7!2:.Me (p41) ++%o(m) lila,nay sam 
ôVg 5(0KG,) Jibs ( p)! Vicptl VE apd — Vind ’ 


oESm p= 0 P 


24 Chapter 2. Functional Renormalization Group 


and the diagrammatic representation is given by 


Ni, Nim 
Nis Ninse 
partitions ny n 
m 
Ni Nips Nhe eek Mm-1 


where the empty vertex with @ lines represents the interaction vertex Veo , and the line with 
the blob represents the single scale propagator KG, = (Gj>'+R)7!-6R-(G5'+ R)+. In the 
first diagram in the right hand side, the summation over partitions means that we should sum 
over all the partitions of the labels {1,...,m} into two subsets {71,...,2)} and {ipy1,-.-,¢m}, 
which satisfies 7) < --- < i, and ip4; < > < îm. In these graphs, all the external lines are 
amputated, which represents the fact that the Wilsonian effective action Sint, represents the 
generating functional of the amputated, connected vertices. 


2.3.3 Wick ordered formulation of the FRG 


In this subsection, we shall consider the Wick ordering, which can remove some divergences 
in the continuum limit. Although this procedure does not cancel all ultraviolet divergences, 
the Wick ordering of operators would be natural because renormalization procedures should 
derive not only the low-energy effective action but also observables in the effective theory. 
That is, we also need composite operator renormalization and the improvement with the Wick 
ordering gives better convergence in the continuum limit. The renormalization group differential 
equation in the Wick-ordered form was first introduced by Wieczerkowski [29]. 

At first, let us define the Wick ordering and consider its properties. For simplicity, we here 
again assume that fields @ are bosonic scalars. Suppose that a positive-definite covariance Č is 
given. The free field theory with the covariance C is defined by the following path integration: 


1 
(oles) Aleno = f Does (=56 C8) oler)---0(en) (2.3.22) 
The Wick-ordered monomials of n degree are given by 
166 
OQc|O(21) --- O(an)] = exp (-335035) (21) +++ P(Ln). (2.3.23) 


Clearly, for any permutations o € Sn they satisfy 


Qe [b(21) ++ Azn] = Q0 [b(ory) lo); (2.3.24) 


and the set of the Wick ordered monomials spans the algebra of totally symmetric polynomials 
since the highest degree in Qo(¢(x1)--- @(@n)) ism. An important property of the Wick ordered 
monomials is orthonormality in the sense that 


(Qe [P(a1) ++ P(@n)] Me lln) Pm) = Fam || CEO, yl). (2.3.25) 
The easiest way to show this equality is to use the generating functional. Let us define 


2f] =O = f dea aad o ea h (2.3.26) 
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then Eq.(2.3.25) can be written as 
(Qo [e$] Qc [e*0)] dy = exp(fCg). (2.3.27) 
Let us show Eq.(2.3.27). Since we can easily obtain that 


“1/16 ð 21/1 p 
2S J (-3 3 gC ; =) ey q (-50en) el) = e™CI)/2+4U) (2.3.28) 
n=0 n=0 


(Qe [er] Qc [e%] Ta — e7 FOF)/2-(9C9)/2 (elf +9) o — e FCA)/2—-(9C9)/2 olf +9) +9) (2.3.29) 


which implies the result. Especially, we have proven (2.3.25). Furthermore, we can readily get 


ô 
QolO(@1) +++ O(@n)]QclO(Y1) +++ PYm)] = exp (0; 56, 3o) Melote 1) ++ b(En)O(M1) +++ Alum), 
(2.3.30) 
where = is a functional derivative operator acting only on $(21)--- (£n), and A similarly 


acts only on (y1): -- Plym). 

We can now obtain another diagrammatic representation for the Polchinski equation (2.3.18) 
by expanding the Wilsonian effective action Sint«|@] in terms of the Wick ordered monomials. 
We consider the Wick ordering in terms of the low-energy free propagator C'a, where 


dĉp et? p(x—y) 


page NY (2.3.31) 


Calz, y) = (KaGo)(zx, y) = ls ( 


Under infinitesimal changes of the effective cutoff scale A œ> Ae~, the effective action Dit Fa 
changes as 


1 ÔL Sint Klo] ÒL Sint K lo] Ô? Sint K [o] 
Sin = orem ice ae ee) 2.3.32 
S181 = F.C ee ee e oe 
where 0; = —AO,. Let us expand the Wilsonian effective action Sint, so that 
Sint, [9] -yt ae mr [bn na = Png |: (2.3.33) 
j= =i 7 
Let us calculate the left hand side of the Polchinski equation (2.3.32). Since 
TA E EE Oop aed (2.3.34) 
t84Cn [YP n1 nj] T 2 Anm 5b ddan CFn nj Jo os 
the Polchinski equation for Wick-ordered vertices Va becomes 
2 Tag N TE 130r, von Sint aea LO OSs, rea lAl (2.3.35) 


j= iJ On ÒPm 


The expansion of the Wilsonian effective actions in the right hand side of Eq. (2.3.35) produces 
the products of the Wick-ordered monomials, and then we would like to rewrite them in terms 
of the Wick-ordered monomials by using Eq. (2.3.30). Substitution and comparison of the 
coefficients give us the following result: 
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Theorem 2.3.2. The Wick-ordered vertices Vy satisfy the differential equations 


: (Chii e Chipi) 


YUM =, 1 = 
OV j S2 Ge m)! 


k=1 m=0 0€S; 
T7No(1),.. No (m), tly otk Snol (m+1) J: eono (J) st est, 
x Vimtk Va (j-m)+k . (2.3.36) 


Let us consider the case where Zə- E ot —¢ exists, then diagrammatic expressions 
of Eq. (2.3.36 (2.3.36) ) for small number of legs (j = 2 and 4) are given as 


OVi2 = —@---O@— + © + + Eb + Q 2, (0337) 
ans ood + KW + Oe + QE +. esa 


Here gray vertices represents Wick-ordered vertices, the usual lines represents the free propa- 
gator Ca = KaG, and lines with a blob represent the single scale propagator 0,K,G,. Again, 
all the external lines are amputated. 


2.4 Optimization 


In principle, any formalisms of FRG give the same result and the result does not depend 
on a choice of regulators Rp, which suppress the fluctuation of low-energy modes in the flow 
equations. However, approximations which we will apply breaks this equivalence and the choice 
of regulators R can become important. Let us discuss optimization of the RG flow, which is one 
of the most important problem in practical applications of FRG to realistic systems. Usually, 
we have to adopt some approximation schemes to reveal interesting physics, and we will rely 
on the approximation if we have a systematic procedure of that approximation scheme. Here, 
we should notice that even if we have a systematic approximation, which can be applied at 
any order, we have to truncate the expansion at some finite order in practical calculations. 
Therefore, in approximate calculations physical quantities can depend on unphysical degrees of 
freedom, such as the unit of mass upr in the perturbation theory, or a choice of regulators Rx 
in FRG, whose dependence should vanish for observables when approximations are absent. 

Although it has not been formulated in a mathematically rigorous way, it is very important 
to seek the optimal choice of the regulator Rę in a given truncation scheme. In the context 
of FRG, such optimization procedures are proposed by D. Litim 30}[81], and he applied his 
procedure to the local potential approximation (LPA) and find the optimized regulator R;(q) = 
(k? — q*)O(k? — q?) for boson propagators and R;(q) = ¢ (y: Eo 1) O(k? — q?) for fermion 
propagators in a relativistic system. On the other hand, the agile of minimal sensitivity 
(PMS) is one of the acceptable criteria for optimization of an approximation scheme, which 
was first introduced by P. M. Stevenson in the context of the perturbation theory [32]. In the 
review [26], J. M. Pawlowski indicates that the optimization criterion given by D. Litim can be 
understood in the PMS. 


2.4. Optimization af 


2.4.1 Optimization criterion 


We define the optimization criterion which tells us how to choose the regulator R which is called 
the optimized regulator. In the next section, we will see properties of the optimized regulator. 

In this discussion, we assume that regulators Rọ are monotonically decreasing as k decreases: 
Ry < Ry for k < k’. At first, we define effective cut-off scales keg, which represent mass gaps 
of IR regularized theories. 


Definition 2.4.1 (The effective cut-off scale). Choose and fix a regulator Rpase,x which 
represents a specific flow. The full propagator with a regulator Rẹ is given by 


1 
G[¢, Rx] = m > (2.4.1) 
Tle] + Re 
and we define the norm ||G[Rx]||sup by 
[|G[ Re] | [sup = sup ||Gly, Rill z2- (2.4.2) 
p 


Let dimG be the mass dimension of the full propagator, then for each k € [0, A] we define the 
effective cut-off scale keg(k) by [1] 


kmc (k) = | |G Rbpase,x]||sup- (2.4.3) 


Let us consider a concrete example and become familiar with it. Here we take a scalar field 
theory and adopt the LPA, in which we put an ansatz for flowing effective actions so that 


1 


redy = 3 | P-pP° Pp + Vly]. (2.4.4) 


In order to define keg(k), we have to fix a regulator Rpasex- Here we use the sharp cut-off 
regulator (2.3.8), which gives us an explicit physical interpretation for cutoff scales k: 


1 
ase = LUshar =p —— ~ — 1). 2.4. 
Ry (p) R h pk (p) P (some _ 1) ) ( 5) 
Then we can find that 


Gly, R] = (P/O — 1) + VPI) 


Therefore, we get 
-2 sarda 2, y2) \7 
Keg (k) = sup (k a Vie (el) = (k + a y (2.4.6) 
p 
From Eq.(2.4.6), we can say that keg(k) gives the mass gap of the regularized theory with the 
IR regulator Rsnarp,k, and indeed at k = 0, ke#(0) = Vass is the mass gap of the full theory. 


"If we consider situations, where the renormalization of the full theory depends on the regulator R, we have 
to divide the left hand side of Eq. (2.4.3) by the field renormalization Z. Here we fix the renormalization scheme 
of the full theory and use the same renormalization procedure for theories with IR regularizations. 
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Now that we have another parametrization of regulators R;, by the effective cut-off scale 
keg, we can consider hyper-surfaces {R,}x., whose elements regularizes the theory at the same 
physical cut-off scale keg. Restricting variations of regulators into this hyper-surfaces, we can 
consider the variation of regulators in the same keg. Then we can require that the regulator Rk 
should behave as an IR regularization. 

In order to define the optimization criterion as a stability condition, we need to study 
functional dependence of correlation functions [;,[y] on a set of regulators {Ry}. The flow 
of Talo] is given by Eq. (2.2.19), and since the operator is well under controlled at the starting 
point k = A we find that 


= > ^ dk'1 l =(2)nm 
intel =Talel-+ f F5 Clo Ruhu Glo, Redon ite) (2.4.7) 


Therefore, we can change the point of view: Instead of regarding J;,[y] as a function of a 
regulator R we can consider A as a functional depending on flows {Ry}, and on the position 
of the flow k € [0, A]. To express this in an explicit way, we write the correlation function as 
Talo, R], where the argument R represents dependence on the set of regulators {Ry }xejo,a] and 
the subscript represents the position of the flow. At k = 0, the functional dependence hlp, R] 
on R disappears without any truncations. 


Definition 2.4.2 (The optimization criterion). Let keg(k) be an effective cut-off scale. The 
hyper-surface {R, }k„ with the same physical cut-off scale is defined as 


{Ri hte = {RelVk € [0, A] ||GLRa]||sup = keg" (A) }- (2.4.8) 
The optimization criterion in a closed form is given by 


Dr Ñl, Fl] = 0, (2.4.9) 
—=tstab 


~ ô 
where Ip|p, R] are correlation functions. Here Dp, = f dk’ Lk 5 Beam represents the total 
k! 
derivative in hyper-surface directions. Since the starting point and the end point of the flow is 
fixed, OR, w = 0 at k’ = 0, A. This condition is already non-trivial even without any truncations 


and Rstab is called the optimized regulator. 


By changing the point of view via Eq.(2.4.7), we have reached the idea of stability under 
small variations of a “flow” but not of a regulator. In order to realize this idea, we have to 
construct the space of the flows which can be achieved by parametrization of regulators with 
respect to the effective cutoff scale. This is done in Eq. (2.4.8) so that flows give the same mass 
gap at the same scale parameter k. Thus, Eq. (2.4.9) is the requirement of the stability of the 
flow. 

In this definition we require that dR o = 0. However, we should notice that Tyly, R] only 
depends on Ry with k’ > k according to Eq. (2.4.7). Therefore, we can remove the condition 
Ro = 0 from the optimization criterion. In other words, that condition is already included 
when we have restricted regulators Ry as elements of {R1 keg. 
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2.4.2 Properties of the optimized regulator 


At first, let us rewrite definition (2.4.9) of the optimized regulator Rstab in a suitable form for 
practical calculations to construct the optimized regulator. This can be easily done after the 
change of viewpoint for the flow equation of correlation functions. 


Theorem 2.4.1. The optimized regulator Rstab satisfies 


(Gly, RJ ReGly, Rel) yow 
Rem k 


on ly, R] = 0. (2.4.10) 


R=Retab 


ô 
Here we have defined the notation Dr, = RY” 3Ram and by omitting the subscript k we 
emphasize that it represents a variation of the flow. 


Proof. We should notice that the derivative with E a parameter 0, and Dr, commute 
with each other. Therefore, by applying 0; to Eq. (2 we can find that 


fe (ssn onic te. Ri] ) =0. 


R=Rstab 
By taking functional derivatives of Eq.(2 , we can get Dr, I, Te Me, R] 


can obtain Eq. (2.4.10). 


Remark. This theorem indicates that if we apply the optimization criterion ( (2.4.9) ) for arbi- 


= 0, and then we 


Rstab 


trary correlation functions T the variation of GRG in the hyper-surface {R1 }k,ẹ should vanish 
pointwise, i.e., for all k and for all indices. This would be an over-constraint and in practical 
applications we restrict the set of J used in optimization. 

_ As we have pointed out in the remark, we have to restrict the set of the correlation functions 
Ik, to which is applied Eq. (2.4.10). Then, which J; should be taken for the optimization 
procedure? Here we should realize that in any approximations we believe that there are finite 
number of the relevant operators {Ie} and that they plays the dominant role in physics. 
Therefore, we should take the correlation functions {iat which is considered relevant in the 
approximation scheme we use. 


Theorem 2.4.2 (Principle of minimal sensitivity). Let us assume that a flow satisfying 
Eq. (2.4.10) exists. Then, Eq. (2.4.10) implies the PMS condition on the correlation functions 
{Ia} used in the optimization criterion: for any T € {Ia} we can find that 


Toly, R 


Rstab 


Proof. We should notice that the starting point of the flow Lää is R-independent. There- 
fore, taking the variation of Eq.(2.4.7) with respect to the regulator R, around the optimized 


regulator we obtain 
dk’ 61" (0, R] 
=; sf —(GRG)nm OR, 2 


Rstab Rstab 


JACA 


SRY” S Ram 


Since this is a linear equation for the correlation function Talo, R] and the left-hand side is zero 
at k = A, we can say that the left-hand side is identically zero at any k. Especially at k = 0, 
we have obtained the result. 
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2.4.3 Example: the Litim regulator in LPA 


Let us consider a scalar field theory in the local potential approximation (LPA): 


Vly] = = f P-pP’Pp + Vily]. (2.4.12) 


2 
In this approximation, D.Litim suggests that R,(p) = p?(k?/p? — 1)O(k?/p? — 1) is the 
optimal choice, which is often called the Litim regulator. Here, following the review given 
by J.M.Pawlowski, we derive the Litim regulator from the optimization criterion (2.4.10). In 
the following, we should note that all the quantities in the LPA are evaluated with constant 


fields @. 
Due to the momentum conservation law, we can put 
TP [p (p, a) = Tlp, p) (21)°8 p — q), (2.4.13) 
and in the LPA we can set 
= 1 (27)*8 (p — q) 
Gp, R] = -aa p ODs a = (2.4.14) 
Di [P] + Re P? +V [el + Ri (p?) 


Then, Eq.(2.4.10) tells us that 


dp dĉq n ô IRP) = 
Jae | ae a : + RP) +V T mee 


By the dimensional argument, we can put Rz(q) = qr(q?/k?), and then we introduce the 
variable x = q?/k?. At first, we rewrite the q-integration into the x-integration. 


ae ec Lil, 2)OARK(T) °° gwidre Te(G,2)(—207 der (x) 

kd / q p q =. / 9 (x ae + V2) fale 

=—-O,q i rY tdr Talp aort) : 
i (1+ r(x) + Vy" [pl/2)? 


=0. (2.4.15) 


Rst ab 


(27)° (q? + Rag?) + Ve)? 


(24 is the solid angle of the d-dimensional space. For the integration by parts, let us consider 


ù ”  1+r4V@/xe Jo (l+r4+V@/x)? = 


dart? a ln Z; E ` 
+f O END PEDRA N a a E 


The left hand side becomes —6y2T;(@, 0) since r(x) and V®) /x diverges as x + +0. Therefore, 
we should search extrema of the quantity 


o0 d r+V® /x V® far 
ba2Tx(@, 0 Wee | het pn, 7 
aoLk(P, +f dra k (G ) +00 a el lt+r+V@/x a) 


(2.4.17) 


under variations of r in the r, direction with Z;, fixed. 
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Then, we have to calculate the r-derivative of the parenthesis in the integrand of Eq. (2.4.17) 
and we get (2.4.17) 
[(d/2 —1) + zô, nZ,)(1+7r+V@/x) +V? /x 


(l+r+V©@) /x)8 oe) 
Therefore, we can rewrite Eq.(2.4.15) as 
x dezt? trt? -2I ór Kaam 1) + 2d, MLL +r + VO /2) cl VO fx =0. (2.4.19) 
0 eee (l+r+V0)/x)3 _ 


Tstab 


Let as assume that the quantity of is positive for the most relevant operator J, such as 
T,=a,|[~] = °. This is not a trivial assumption because Vp is not necessarily convex for k 4 0 
and for a non-convex bare potential V,,. Furthermore, we cannot ensure that the inequality 
r+V)/x > 0 holds for any x and G due to the approximation. Here we just expect such a subtle 
region only has a minor contribution to the result for optimized r, and that this assumption 
can be justified in some way | Then, Eq. (2.4.19) says that we should minimize or maximize 
the regulator r in the space {r1},., so that dr, = 0 at the optimal regulator r = rgtap. For the 
stability of the flow, we should minimize the regulator r since it indicates the minimum of the 
quantity and then the minimal flows. On the other hand, maximizing the regulator r 

corresponds to the most instable flow. 
In above calculations we have not taken into account that r belongs to {r1 }k,ẹ. From now 
on, we have to use it in order to proceed the calculation further. By definition (2.4.8), we get 
Va € [0,00] x + eri (x) + V| >14V2 


min — min? 


(2.4.20) 


Eq. (2.4.20), we find that z + xr,(x) > 1. Therefore, in the region x > 1, we can set rı = 0 
and this is consistent with the condition of Eq.(2.4.19), which indicates that the minimal value 
of r is favored. And, we can saturate the ne at 3g with r(x) = 1/x — 1 when z < 1, 
which implies that 


where we have used keg given by Eq. (2.4.6). By subtracting VČ) from the both sides of 


Tstab(Z) = (1/x — 1)@(1/z — 1), (2.4.21) 
which is nothing but the Litim regulator: 
Rinad) =g ra r = (k? =g =): (2.4.22) 


Then we have obtained the optimized regulator for the scalar field theory in the local potential 
approximation. 


Remark. It is very important for the practical applications that the Litim regulator 
has a very simple form. Its derivative 0,R;(q) = 2kO(k?/q? — 1) also does not introduce any 
singularities in the loop calculations, and then the low-momentum region would give a main 
contribution which is a desirable property for the LPA. Furthermore, its simplicity often enables 
us to calculate analytic expressions for the Wetterich equations, and then the numerical effort 
can be reduced significantly. 


8We can see that this statement is glossy. In the double well potential, V?) is most negative at Y = 0, 
which defines the scale of the physics. Therefore, |V? | should be much smaller than the UV cutoff A2 and then 
r+V©)/x is positive for x <1. For large momenta x > 1, this may not be true but such regions should be 
suppressed in the flow equation due to the effective UV cutoff. In Eq. (2.4.19), this effect should be taken into 
account through the factor dr_. 
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2.5 Properties of flow equations 


In this section, we discuss general properties of the exact solutions of the RG flow. These 
discussions were originally developed in the study of critical phenomena [6,{33}]36}, and now it 
becomes an important tool also for quantum field theories. 

The Wilsonian renormalization group consists of two steps: coarse graining (or decimation 
of modes) and a rescaling. In this thesis, we have only considered a coarse-graining in deriving 
flow equations since the rescaling is not necessary for calculations of Green functions. However, 
we will see that the rescaling is necessary to find fixed points of the RG flow and it reveals 
many important properties of exact solutions. 

For simplicity, the discussion given below is restricted to the d-dimensional scalar field theory 
within the Wetterich formalism; arguments in a more general context can be found in (35). The 
description given here is strongly influenced also by the papers [87H40]. 


2.5.1 Rescaling 


With decimation of modes or coarse-graining, an effective cutoff k is reduced to another effective 
cutoff k’ = e~tk. In order to consider fixed points of the RG flow, the lowered cutoff scale should 
be brought up to the original one and then we need the rescaling. Let us consider the dilation 
transformation 

p p = ep, z x =e x, d(x) 6 d(x’). (2.5.1) 


Here, dy = (d — 2 + 7) is the scaling dimension, which is to be determined in searching for 
a fixed point of the RG flow. 7 represents the deviation of the scaling dimension from the 
canonical dimension of the field, and it is called the anomalous dimension. In the momentum 
space, the field transforms as 

Pp > ead) gy, (2.5.2) 


Therefore, under the dilation, changes of fields @ are given by 
Spp = C4) 6 sey — bp = Ôt ((d — dy) + pYOpu) dp. (2.5.3) 
We can find that the dilation changes the functional I'[¢] by 
ôT = 6tGal [¢], (2.5.4) 


where Gai is the operator given by 


7 Ob» 6 5 
Gail = ( J pae + (d dy J hz) Tig. (2.5.5) 


The first term in the parenthesis represents the effect of the scale transformation on momenta, 
and the second one represents the effect of the scale transformation of fields. 

In some literatures, other conventions to represent the dilation transformation Gg; are 
adopted. Here I show two other forms which are frequently taken: 


e By integrating by parts in the first term in Eq.(2.5.5), we can rewrite it as 


ô ô 
g= ( J PO se + dy i tz) r. (2.5.6) 
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e Let us introduce the new derivative ô., which acts only on coefficient vertices, i.e., does 


not act on the momentum-conserving delta-function. With this derivative, we can obtain 


that ; 
Gail = (« - | era A > -ds f doze) F; (2.5.7) 
p 


This follows due to p- 0,54(p) = —dd“(p). 


Here we have denoted a traditional way of the rescaling. We can achieve the same conclusion 
in a neater way by making all the dimensional quantities dimensionless with the cutoff k. 


2.5.2 Fixed points 


We consider the renormalization sronp equation with rescaling. Let us denote the total deriva- 
tive with respect to t = — ln k/ko by 4 = oe to represent Wilson’s RG procedures, and we 
denote the effect only from the first step by 2 3: We have 


dr, — Oly 
dt ðt 
in which the first term is given by the Wetterich equation (2.2.11), which ay tet the effect 


of coarse-graining, and the TE s term represents the effect of rescaling (2.5.1). Let us denote 
the right hand side of Eq.(2.5.8) simply as G(7,T',) to get 


d 

—l,= Tg): 2.5.9 
Here we make the dependence on a choice of the anomalous dimension 7 explicit. It is important 
to notice that G acts on T in a nonlinear way. 


A fixed point r, with the anomalous dimension 7, is given by 


G(T.) = 0 (2.5.10) 


+ Gail’, (2.5.8) 


d 
so that it satisfies a = 0. That is, the fixed point represents the scale invariant theory. In 
the vicinity of T, we have an expansion in terms of perturbations AT (t) so that 
d 
mri 
dt 


where we have separated the right hand side into a linear part £, a quadratic part Q, etc. 


Tr, + AT(t)) = LAT(t) + QAT(t) +---, (2.5.11) 


Definition 2.5.1 (Scaling exponents and scaling operators). Consider the eigenvalue 
equation: 


LOH XO>. (2.5.12) 


The eigenvalues A; are called scaling exponents and the eigenoperators O; associated with A; 
are called scaling operators. 


Let us assume that the set of the scaling operators O; spans the full vector space. Under 
this assumption, we can write any actions [’,, near the fixed point T, as 


Te = D, + >) ai(t)Oi. (2.5.13) 
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a;(t) are nothing but effective coupling constants, sometimes called scaling fields [84]. In the 
linear approximation, they satisfy 


—a;(t) = à;a;(t), (2.5.14) 


and then we find that 
Te =T, +X aie. (2.5.15) 


Recall that the limit k — 0 corresponds to t + oo. The following classification makes sense: 


Definition 2.5.2 (Classification of the scaling operators). We classify scaling operators 
into three different types as follows: 


e If å; > 0, the associated scaling operator O; is called relevant, which emanates I, from 
the fixed point I, in the limit t > oo. 


e If \; < 0, the associated scaling operator O; is called irrelevant, which drives I’, into the 
fixed point I, in the limit t > oo. 


e If A; = 0, the associated scaling operator O; is called marginal and I, + a;O; is also a 
fixed point for infinitesimal a;. 


Remark. In order to describe critical phenomena, we should require that the effective action I, 
approaches to the fixed point I’, asymptotically, and then we have to fine tune parameters so 
that the relevant directions vanish. On the other hand, in order to define the continuum limit 
of the quantum field theory only the relevant operators dominate low-energy physics and then 
the irrelevant operators should be put to be zero (see Appendix [B). 


Especially when a marginal direction appears in leading order consideration, it is very 
important to take into account higher order corrections. Let Omar be a marginal operator, and 
we consider a perturbation Fly] = rly] + Pily] with 


Pill = xt) Omar ly] + 2 r(t) Oly], (2.5.16) 


where scaling fields u;(t) should be considered as quadratically small quantities. In this case, 
we consider up to the quadratic order of the perturbation P;. One possibility is that at the 
quadratic order of the perturbation x(t) shows the t-dependence, and in this case Omar is called 
a marginally relevant /irrelevant operator. In this case, scaling fields j1;(t) should be chosen so 
that the perturbation P; is an eigenperturbation in this order. 

Let us suppose that the flow in the Omar direction vanishes even in the second order per- 
turbation. In this case, we should ask whether the eigenperturbation P, itself vanishes in the 
quadratic order or not. Extracting quadratic order terms of the flow into the O; direction, we 


find that d 
ge = Api (t) + bix(t)’, (2.5.17) 


where the first term comes from L(u;(t)O;) and the second term comes from Q(x(t)Omar). If 
the feedback b; is absent, the flow is eliminated simply by taking u; = 0. Even if the feedback 
exists, the flow of u; vanishes by choosing ju;(t) = —b;x?(t)/A; as long as A; # 0. Therefore, if 
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the mixing with another marginal operator does not exist, the flow of the perturbation can be 
eliminated up to the second order perturbation. If we can kill the flow along the Omar direction 
at any order in this way, the operator Omar is called exactly marginal. An exactly marginal 
operator generates a line of fixed-points. 


Remark. From now on, we have discussed the notion of fixed points in the Wetterich formalism, 
but all the statements given above is of course true for other versions of FRG. However, we 
should make some comments on the relation between the Wetterich and Wilsonian effective 
actions. 

Originally, F. J. Wegner discussed properties of critical fixed points based on the invariance 
of the partition function under RG transformations [33135]. Therefore, his arguments can be 
directly applicable to the Wilsonian effective action, and then fixed points in the Polchinski 
version of FRG have well-known properties. That is, if it is a critical fixed point then a line 
of equivalent fixed points exists, and the line is generated by a marginal, redundant operator. 
This property is well studied by O. J. Rosten [40]. 

On the other hand, the connection between the Wetterich formalism and the Polchinski 
formalism is not still clear when the rescaling operation accompanies coarse-graining in RG 
transformations. Recently, a clear connection between them is given just above critical fixed 
points by O. J. Rosten [39], and it is explicitly shown that a line of equivalent fixed points also 
exists in the Wetterich formalism. 
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Chapter 3 


FRG for many-body fermions 


In this chapter, we will consider interacting many fermions using FRG. In deriving the low 
energy effective theory with the renormalization group method, identification of low-energy 
excitations is most important but one of the most difficult problems. Indeed, there does not 
still exist a systematic way to derive low energy excitations if it were composite particles, 
topological solitons, etc. In the case of fermionic systems, the identification is still difficult even 
if we assume that low energy living particles are fermions given by elementary fermionic fields. 
In the Fermi liquid picture, such excitations are given by pseudo-particles in the vicinity of the 
Fermi surface but the calculation of the Fermi energy from the microscopic model is not trivial 
at all. 

In this chapter, we at first define the setup of physical systems to be considered in a general 
way. In later chapters, we will treat specific models, e.g., BCS-BEC crossover and dipolar 
fermions, which are contained in the general discussion given in sec[3.1] In sec[3.2| we will 
derive the RG flow equation up to some lower orders in terms of the field expansion of the 1PI 
effective action. 

In order to get a physical intuitions for the RG flow in fermionic systems, we will review 
some previous studies on the Landau Fermi liquid theory in terms of the FRG approach in 
sec|3.3] After that, we consider the validity of such approach and pose some problems in 
treating the Fermi liquid from the microscopic theory in sec[3.4] Especially it will be turned 
out that the self-energy correction is important for calculations of the Fermi sphere for the 
interacting fermions, because we need to identify gapless excitations to separate high-energy 
and low-energy scales of physics quantitatively. 


3.1 Notations and setup 


In this section, we describe notations and the set up for many-body fermionic systems which 
we will consider. 

Let Wa(x) be fermionic fields with the flavor index |'}a € {f,4} and with the spacetime 
label x = (7,a). Here 7 € |—8/2, 8/2], with the inverse temperature 6 = 1/T, represents the 
imaginary time and æ € R? runs over the spatial directions. We require fermionic fields to be 
antiperiodic in the imaginary time direction: 


Yo(B/2, £) = —Ya(—B/2, 2). (3.1.1) 


'In this section we treat this label as a spin index, but sometimes it represents another degrees of freedom. 
For example, in the case of the BCS-BEC crossover this label represents hyperfine states of an alkaline atom. 
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For the simplicity of the notation, we sometimes extend the fields y(x) as an antiperiodic field 
on R+! with the period £ in the imaginary time direction, that is, Y(T + 6,2) = —W(r, £). 
The conjugate field (2) is also introduced and satisfies the antiperiodic boundary condition 
Y(T Ep; z) = =r, a). 

We assume the homogeneity and thus the momentum space representation is quite useful 
in calculations of physical quantities. For that purpose, let us define the Fourier transform of 
fermionic fields as 


Upa = dre tals), Ypa = f doe =p (a), (31.2) 


where p = (wn, p) and pr = wpT +p: x. Here the temporal component p° is given by fermionic 
Matsubara frequencies wn = 7(2n+1)T with n € Z. The difference of signatures e*’?* in 
can be understood by regarding w as a linear functional and yw as an antilinear functional: 
wef) = (f) and y(cef) = c*(f) for a test function f and a complex number c. The inverse 
Fourier transform of them is given by 


(T). o Do 
TE / ey), B(x) = i, cm (3.1.3) 
p p 


T) 1 d¢p 
where == ` —— and the sum is the fermionic Matsubara sum. 
p BÈ] QnA 


Let us dictate transformation properties of fermionic fields = and %. We will require the 
symmetry of the action under the following transformations: 


e Global U(1) symmetry: Let Yy + ew and Y +> ew with a constant a. The chemical 
potential u defined below is the Lagrange multiplier corresponding of this symmetry. 


e Rotational symmetry |} Let us consider infinitesimal spatial rotations w+ æ + O(n x a) 
around a unit vector n, then fields transform as 


wry! =exp(—-ién- S)v, pH Y = wexp(idn - S), (3.1.4) 


where S = Ș is the spin matrix and ø are the Pauli matrices. Under this transformation, 
w' Sw transforms as a vector field. 


e Euclidean time reflection [} Let us consider the time reflection 0 : (T, £) + (—7,a), and 
we define the antilinear operation on the field algebra defined as 


Op(T, 2) = Yal =T, x), Ob, (7,2) = Va(—T, £), (3.1.5) 
and 0(AB) = (0 B)(0A) for functionals A, B. 


e Parity: The discrete symmetry given by Y(T, Œ) œ> Y(T, —x) and by yt (T, x) = yt (T, =x) 
is called the parity. 


2In the case of the BCS-BEC crossover, the labels a do not represent real spins, and we further require that 
the action is invariant under the spin SU(2) rotation. 
3The invariance under the Euclidean time reflection is strongly related to the Hermitian property of the 


Hamiltonian ; 
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To calculate the thermodynamic quantities, the imaginary time formalism will be used and 
the grand partition function Z is given by 


2(8,n) = Te [exp(—A(H - nN))] = | DEDWexp(-SH, vi), (3.1.6) 


where H is the Hamiltonian, N is the number operator, 8 is the inverse temperature, and LL is 
the chemical potential. S is the corresponding microscopic action, which is a functional of Y, a, 
and the path integral in Eq. (3.1.6) is given by the Berezin integration. As usual, we decompose 
the action © into the free part So and the interaction part Sint. The free quadratic part is 


(T) 
S.(0,0] = f AO AE E (), (3.1.7) 


where the kinetic energy is €o(p) = p”/2m. This free action satisfies the symmetries above, 
and furthermore it is also symmetric under the spin-SU (2) rotation. 


3.1.1 Possible forms of the effective action 


Let us write down the interaction part Sint. Although we require later that Sint is quartic at 
the classical level, it is important for us to consider about the most general form of actions 
since such terms will appear in the quantum effective action. The global U (1) symmetry and 
the momentum conservation suggest that 


X% i (T) 
Sm = zr | (pi ++ +++ Pn — G1 — +++ — Gn) 
n=1 (n!) Pl yeesPn5Q1 5-05 dn 
xD erg (pı, see > Pn; fns- +) G1 Pex a on an Van Ba i Pbs (3.1.8) 


where 5“*!(p — q) = 86,0,q0(27)76"(p — q). The Fermi statistics of fields gives totally antisym- 
metric properties to vertices: for any permutations o,o’ € Sn 


Ao (1)++-Oa(n) 


Q1...An 
k Bot (ny-Bat (1) 


Bn. BY (pı, -<+ Pni fns- -> 1) = sgn(o ` o'I (Do(1)s +++) Poln); Jol(n)s ++ +5 do'(1)) (3.1.9) 


Using the parity invariance, we get the equality 
Daran (Dr... , Pri gns- -3 Q1) = T320 (Pp, ---; PPni Pans-+-; Pa); (3.1.10) 


where P is the parity operator: P(p°, p) = (p?, —p). The invariance under the Euclidean time 
reflection implies that 


Daage (Prs Pri Gns 1301) = Darien Oats + + s Oqni Dns «+ 0p1)]*, (3.1.11) 
where * represents the complex conjugation and 6(p°, p) = (—p°, p). 
Especially, we shall pay attention to the self-energy correction and the four-point vertex. 
Let us specify useful properties of them concluded from these symmetries. For the self-energy 
Dsp) = —G(p; p), the rotational symmetry suggests that 


=3(p) = X(p°, p?)6g + Us(po, p?)p - S3. 
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However, the second term does not satisfy the parity invariance and then it vanishes as long as 
the parity symmetry is not broken. We conclude that the self-energy takes the form 


x3 (p) = X(p", p?)dg, (3.1.12) 


where }(p°, p?) = X=*(—p®, p°) due to Eq. (3.1.11). 
Next, let us consider the properties of the four-point vertex. Let us write down the momen- 
tum dependence of the four point vertex as 


Vea (Pq, d) = T38 (P/2 +q, P/2—9¢;P/2—q',P/2+¢), (3.1.13) 


where P represents the total momenta and q + q’ represent the momentum transfer. Let us 
decompose it with respect to the total spin S. Components with S = 0 consist of two terms 
aQ dd a1 5a ss aQ a 

Tan (Png) haa = F°°(P3 4, ¢/) 626g? + T*(P; ¢, ')(Sgi - $32), (3.1.14) 
where those coefficient functions are scalars under spatial rotations due to the parity invariance. 
Since those two terms are mixed under interchange of particles, the consequence of Fermi 
statistics is unclear with this expression. In order to make it obvious, we rewrite (3.1.14) using 
Fierz identities |*]so that 


Teer Prge] 


BoB =P Pa el eae (3.1.15) 


S=0 88281 a B2 By 


where Pegg; = O Bo + Òg, Pe x (a;e1)2(€0;) 4,4, is the projection onto the spin-symmetric 
Q12 


component and Pagg, = (53°63? — 0a Oar) œx (e1)™%e,,4, is the projection onto the spin- 
antisymmetric component. From the Fermi statistics, we can conclude that 


M(P3q,q/) = -T4(P; —q, q') = -T4 (P; q, —q'), (3.1.16) 


Mirne 9,4) Ea: (3.1.17) 
The Euclidean time reflection requires that 
r4/S(P;q, d’) = [T^ ((-P°, P); (—4°, q'), (9°, 0)" (3.1.18) 


When we can require the invariance under spin SU (2) rotations, the four point vertex only 
contains the S = 0 component and following terms do not appear. 
The S = 1 component consists of three terms 


Tea (Pia d) = TPP, Sigh og? + TH (P; ad)o Sa + TPF (P; a, dS X Se? des 

(3.1.19) 
where the parity invariance requires that those coefficients are pseudovectors. If we would like 
to find the consequence of the Fermi statistics, the Fierz transformation to the particle-particle 
channel again gives clear descriptions: 


ro 


Fag gay = Va e T maa E e a TS x ES) (3.1.20) 


4TIn order to make the rotational symmetry manifest in the particle-particle basis, we use £,£ø as a basis of 
2 x 2 matrices and e', oe! as a dual basis. 
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where the momentum arguments are omitted. Now the Fermi statistics suggests that 


r*(P;q,q) =-=], (3.1.21) 
etc. The consequence from the Euclidean time reflection also becomes clear in this basis: 
[*(P;q,q’) = [I (OP; 0q', 0q)]*, T*(P; 4,4’) = —[E* (OP; 60, Oq)]*. (3.1.22) 
The S = 2 component consists only of the tensor coupling 
Tee? (P30,9)| 5. = rT Pan Ba? (3.1.23) 
where T;; is the tensor operator, given by 
aja 1 a a aiga Oi Oy a 
(Toni = 5 (S299 + SSi) — SECS SR) (3.1.24) 


We will find in chap[5] that the dipole-dipole interaction contains this term at the classical 
level and then the spatial and spin rotations cannot be treated separately. We will postpone 
studying details of this term until chap}} 

Finally, let us comment on forms of possible interactions in the vacuum. In the nonrelativis- 
tic system, we can further require the Galilean invariance in the vacuum. In this case, vertices 
depend on the center of mass momentum P only through the Galilean invariant combination 
iP? + P?/4m, and the phenomenological description of effective interactions becomes much 
simpler. 


3.2 RG flow equation 


In this section, we will write down the Wetterich equation of fermionic systems for the self- 
energy and the four-point vertex. We introduce diagrammatic expressions for convenience: 


Eralp) =—Tea(p) = Q tl (3.2.1) 
Pipi, ph p2 

-Ires (Pid) = -Trap (P1, P25 Popi) = ; (3.2.2) 
Qı Pı P2 Q2 


where pı = P/2 + q, p2 = P/2 — q, p = P/2+ q', and p = P/2 — g. 
Let us write down the Wetterich equation for these vertices. For derivation, we use the 
diagrammatic method described in Appendix[A]} For the self-energy correction, the flow equation 


is LB 


g; (o =m a’) = Op () (3.2.3) 
p p a al’ 
_ P P 
where we have introduced the new derivative 0;, which acts only on the k-dependence of the 


infrared regulator Rẹ in the loop propagators. The analytic expression of the flow equation for 
the self-energy is given by 


” Oore ae (p—1)/2) 
Ope (P) = -3 | e “GH = MO m Rad) (3.2.4) 
(T) e—ilo pR; 
— -f Te32,(p + l; (p — 1)/2, (p — 1)/2) CO sh UTF (3.2.5) 
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Therefore, we get 


= P -1 e P08, R,,(1) 
(G 


T) 
dk Xr(p) = -f [BTE +r] (> +h ; LE, + Re)? (3.2.6) 


The other terms in the four point coupling drop out since they do not satisfy the restrictions 
coming from symmetries. 
For the four-point vertex, the diagrammatic expression is 


— AP ne22(P 34,4) -å ("ROO + a. & FY) (3.2.7) 


The first and the second diagrams represent particle-hole fluctuations, and their sum satisfies 
the restriction coming from the Fermi statistics. We call them ZS and ZS’ diagrams, respec- 
tively. These two diagrams are important in the description of the Fermi liquid theory since 
their singularities for small momentum transfers q +q’ ~ 0 are related to the existence of “zero 
sound” as brought up by L. D. Landau [41/42]. The third diagram represents particle-particle 
fluctuations, which plays an essential role in the formation of Cooper pairs in the BCS theory, 
and then we will call the diagram BCS. The last one comes from feedback from the six point 
vertex. 

Let us calculate the analytic expression of these diagrams. In the following, we assume the 
spin SU(2) symmetry for simplicity. At first, we consider the ZS diagram. Using the Feynman 
rule, we find that 


v A T Pees {P, nd hE rag alts {P. 4, 7}) 
1 2 Ji ([G3!-— Er + Ri]()(G5+ — Er + Rs)(g — g' +1) 


where 
a wap ( P gal Saag =g=1 
mea {P,a d) = rit (g tarni 2 z 5 ; (3.2.8) 
P 1 P 1 
—vya a P —5 +2q- q +l —5> +g +l 
Tsao {P.q,q}) = Dearg ( > q ! l; 2 2 , 2 2 j (3.2.9) 


Q201 


It satisfies relations Dios (l — (q — q); {P,q,q' }) = Pea (0 {P, —q, —q'}) and Piara o == 
q); {P,q4,q}) = Praia, (b {P, —q,—q'}). We can calculate explicitly the summation over spin 
indices. Then we get 


ey 1 A Te ZA ng ! aja 
R T kogoli {P,a d }) = 5 (ores + (CATE + PETA) + PETE) (G (Fg, HPS (3.2.10) 


1 exe ae pe ee 
+5 (areri + 3(T4T5 + rr) — rér¥) iiPad)e 
The contribution of the ZS diagram can be written as 
a ~ p (STT + (DALE + EPA) + rire) G {P,ag 
—ôkr i (P;q, d = -ð / ; 
E(P Dios =—% a O t Ada T) 


(3.2.11) 
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ay (3rATÉ + a(PAV§ + EPA) = PEPE) (G {P0 

-ATPa =- fae bi 2 
p 2(G31 = Xr + Ril) (G3! — Xr + Rilla- g +1) 

The Doi of p ZS’ a is very similar to that of the ZS diagram. 


kaos RIP- a iP, —q'}) (3.2.13) 
a CA See Bet Rilla+ a+ Oy 0 


The summation over spin indices also gives a similar result: 


(3.2.12) 


Goh at, 


1 
raa Tego (GLP, a,—a'}) = 5 [BAD + ATE + LETS) + TEL §] (0 {P,a -0 pP 
1 ~ ~ ~ ~ 
—5 [arava + 34T + Tera) - ring | 0 {P,a -d DP 
As a result, the contribution of the ZS’ diagram is given by 
(514TA + (CATE + DEPA) + PEPE) (G {P,a -0 


-Lr (P: qg, lan =O f , (3.2.14 
EPa a R at ar ES 
~ poo aren + (4DE + Pgs) - Perf) (G {P,e} 
—0,05(P; 4, „=ô . (3.2.15 
TPD EA RONG SE Maller E 
Next, let us obtain the analytic expression for the BCS diagram. 
/ / 
ba aD ARP gl- EPH (Ps 2,q) 
= 3.2.16 
key of (Gt 3s RIGS RP an 
Since ~B6 
Tye Cea i =or4rap satan + ee, (3.2.17) 
the analytic expression of the BCS contribution is 
T Tiel Di erase 
-TÉ (P; q, 9) lacs =a, f| =I a z) a z4) , (3.2.18) 
1 (G31 — Er + Rel) UG! — Ue + Rel(P — 1) 
= [m Teea SPs) 
—ô, r$ (P; q, q' -3 f ae ee a 3.2.19 
TEP a Docs = J TOAS + RI e+ RP Hy 8O 


3.3 Fermi liquid theory from the RG viewpoint 


In this section, we will review the Landau Fermi theory as a low energy effective theory 
arising from the RG method. This point of view is developed in 1990s and reveals a universal 
nature of interacting fermionic systems ae We can find some mathematical rigorous 
treatments of this problem in the textbook |46]. 

Let us consider the low energy effective action of interacting fermions. For this purpose 
we must identify the correct degrees of freedom in low energy physics, but this is a difficult 
problem and depends on details of dynamics. Here we just guess that they are also a spin-5 
fermion and also denoted by w. The quadratic part is again given by Eq. (3.1.7). The ground 
state of this theory is the Fermi sea, with all states e(k) < yu filled and all states e(k) > u 
empty. Therefore, the low lying excitation are obtained by adding a particle just above the 
Fermi surface, or removing one just below. 
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3.3.1 Tree-level RG analysis 


Suppose that we already have integrated out fields with momenta outside the shell of thickness 
2A, around the Fermi surface. For simplicity, we assume that the Fermi surface has rotational 
invariancd?| and we also require the spin SU(2) symmetry. Let k be a momentum in the shell 
then we can decompose it as 

k=kp +2, (3.3.1) 


where kp and £ are parallel to k and kp lies on the Fermi surface. Since A, < kp, the 


No 
momentum integration becomes f dk ~ hir! / dé f dQ,, and then the free action in 


Eq.(3 becomes 
kī 1 
So = Tonya xf. at f aod) [—iw + upl]va(k), (3.3.2) 


where vp is the Fermi velocity. Now we can readily find that Eq.(3.3.2) is invariant under the 
RG transformation: 


1. Decimation: we integrate out modes with momenta lying between sA, < |l| < A, with 
0<s<l. 


2. Rescaling: we recover the original momentum region by changing variables w + sw 
( B — s718), l — sé, associated with the rescaling Y => s~*/2y, Y > s—3/2y. 


Now we have determined the scaling property of fields, then we can judge naturalness of the 
theory by studying the scaling behavior of interaction terms. 

Apart from the momentum delta function, we can easily find the scaling property of quartic 
terms. By multiplying scaling factors of each f, , of a factor 0. +w2,03+w4, and of fields (yip), 
we find 

(s*)* x s7! x (s79/?)4 = s. 


Therefore, if we have totally neglected the scaling behavior of the momentum delta function 
04(k, + ky — ki, — ki), the quartic interaction is irrelevant and then all the interactions are 
irrelevant. This result would be wrong, since we already know that gapless fermions are unstable 
under, for example, the spontaneous symmetry breaking. Indeed we should take into account 
an important subtlety in the kinematics. Momentum conservation requires that 


k = kr(nı +N — ni) + (£ + by — £), 


where n; is a unit vector in the direction of k;. Since the momentum k4 must lie inside the thin 
shell around the Fermi sphere, possible configurations of directions n; are restricted. Especially 
in the limit A,/K yr — 0, the restriction becomes 


lini tng- ni| =1. (3.3.3) 


To make the point clearer, we restrict ourselves to consider 2-dimensional cases in the 
following. The possible solutions of Eq. (3.3.3) can be classified into two types: 


>This is only true for the case where the system is invariant under the spatial rotation. Therefore, the Fermi 
sea of the electrons in conductors does not satisfy this condition because crystals break this symmetry. 
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1. If nı # —Ng, the possible solutions are n; = ni or ng = n}. 
2. If ni = —nz, the direction ng is not restricted. 


These are a restriction on directions and give (d — 1) constraints in the d-dimensional space. 
In this situation, there remains a delta function which only contains @;’s in its argument and it 
scales as s7}. In these special configurations, the interaction becomes marginal and produces 
various interesting physics such as superconductivity. 

The interaction channel which corresponds to nı = n) is called the Landau interaction 
channel and the corresponding effective action is 


1 -D : 
Kl= ioe 2 u dé; Jaf” re (kL, k2; -Q) 


i=1,27—A 


x Waki, cae W2)Wy(k2 — Q, we — OQ)ps(ki + Q, w1 +9). (3.3.4) 


Here we have introduced the following notation for the nearly forward scattering vertex 
Taos (ki, kz; Q) ='a55 (ki, ka; k2 — Q, kı +Q) (3.3.5) 


with the transfer vertex Q = ki — kı = (Q,Q). Since the frequency dependence and the 
dependence on the deviations from the Fermi sphere £ of the vertex are irrelevant, we simply 
put kı = (nT, kL) and ky = (nT, k% ). Nevertheless, we will retain the Q-dependence of the 
vertex since we will find in the 1-loop analysis that the expansion of the vertex around Q = 0 
is not well defined when the cutoff A becomes of order of the temperature T/vp. Therefore, we 
can define the two distinct limits of the vertex functions 


aß , 7 

rR (kr kp) = dim, [r3 eh, k}; Q, D l (3.3.6) 
ap ; a 

Liy (kr kp) = lim [rsp Ky: at}... : (3.3.7) 


In a rotationally invariant 2-dimensional system, P@/° are functions only of the relative angle 6 
between kı and kz. The only remnant of the Fermi statistics is the constraint reao =0)=0 
and re (0 = 0) = 0. For simplicity of the notation, we do not write the subscript A explicitly 
in the n 

In a 2-dimensional case, functions of a relative angle can be expanded as a Fourier series [f] 


= `. Xe’, X, = SEO (3.3.8) 
T 


l=—oo 


The symmetry property X(@) = X(—0) becomes X; = X_; and the consequence of the Pauli 
principle, P°/°) (0) = 0, becomes >, cee, =Q. 


In the 3-dimensional system, this expansion should be replaced by the expansion with Legendre polynomials 
X(0) = ico Xi Pi (cos 0). Due to the difference of the normalization, we should replace X; — X/(2l + 1) to 
obtain formulae for the 3-dimensional case from ones for the 2-dimensional case. 


46 Chapter 3. FRG for many-body fermions 


3.3.2 1-loop analysis 


Let us calculate 1-loop contributions to effective four-point vertices. There are three types of 
1-loop contributions, two of which mainly contribute to the Landau channel and are called ZS 
and ZS’ and the remaining one contributes to the BCS channel 


Y yy n : +> 
ZS = , Z9 = , BCS = oO. (3.3.9) 
1 2 1 2 top 


The analytic expressions can be found in sec[3.2| 

Let us consider the Landau channel at first. For this purpose we set outgoing momenta 
ki = kı + Q and ki = k — Q, where Q is small enough, i.e. |Q| < A. In this channel we 
simply neglect the contribution from the BCS graph because the configuration space of the 
loop momentum ks is strongly restricted due to the center-of-mass momentum P (see Fig[3.1). 
The fermionic Matsubara sum over ws can be evaluated in a standard manner and gives 


1 nr(e(ks)) — nr(e(ks — Q)) 
= S G(ks)G(ks — = — : 3.3.10 
B ” ( 5) ( 5 Q) iQ + elks) — e(ks — Q) ( ) 
Eq.(3.3.10) gives zero in the Q-limit and then 
—A Te(0 0 3.3.11 
DA 12 aes f (3.3.11) 
On the other hand, the Q-limit of the same product is not zero and gives a factor = ; 
4cosh*(Buprés/2) 
With the sharp cutoff, the flow equations (3.2.11]3.2.12) of the vertex in this limit are given by 
B 
A sUFA  N(er) (5 [.9 47? Ags 1 5s]? 
A == Ir? | re reS 4 = ee | (3.3.12 
A*I cosh? (£ vrA) 2 € l +i, l T 5 l ( ) 
= N(ép) (5 A]? A_os 1 sy? 
-S = re*] +are"re?-= [re] }. 6-313 


Here N(ep)/2 is the density of states at the Fermi surface, and here it is given by kp/2rup. 
The similar analysis on the ZS’ graph shows that unless |ki — kə|/kr < A/kr(« 1) the 
loop integration is suppressed due to the restriction on possible configurations. Therefore, in 
the following analysis we simply neglect that graph P] 
According to Eq. (3.3.11) and neglecting the ZS’ diagram, we find that T? is approximately 


invariant under the renormalization group flow: 


T90, A) = POS)" (0). (3.3.14) 


In order to refer each diagram many times, we have assigned names to them. This nomenclature is based 
on R. Shanker’s paper [44], which is also used in [45]. ZS stands for zero sound and we will find that ZS and 
ZS’ diagrams are important for the Landau Fermi liquid. We will find later that the BCS diagram plays an 
important role for the Cooper instability. 

8Due to this approximation, the consequence on the Pauli exclusion principle on the vertex function T9/® is 
violated, i.e., even if we have put kı = k2 or 012 = 0, the vertex does not vanish. For the rigorous treatment of 
this problem, see [47]. 
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Figure 3.1: In the BCS diagram, the possible configuration of the loop momentum is restricted 
into the gray circle. As long as the total momentum drawn with the blue arrow exist, this 
restriction suppresses the BCS diagram significantly. Here we have drawn the figure for the 
3-dimensional case, but the same is true for other dimensions. 


We introduce the dimensionless Landau functions F'*/* by 


F5 = Nr) ar) dt re), Fe = Mee) jp _ TS, (3.3.15) 


They are suitable for the channel decomposition of particle-hole channels. Next let us try to 


solve Eqs. (3.3.12]3.3.13), and for this purpose we introduce the running dimensionless couplings 
A, B by 


is MEF) (grow +709) B= MEF) (paw _ pasy), (3.3.16) 
For these couplings, the flow equations become 
Bop 
—Ad,A; = -—2 FO ap, (3.3.17) 
cosh [SurA| 
Fup 
—Ad,B, 20°F" p (3.3.18) 


cosh? [ZurA] 
At the UV scale Ag, A) = FF?” and Bo = F?". Then at IR scales we find the fixed point 
Fs* ee 


* * 


es ee es 
Po 1+ Fe l 1+ F* 


From the flow equations, Eqs.(3.3.17]3.3.18), and the expression of the fixed point, Eq.(3.3.19), 


we find the stability condition for the Landau Fermi liquid, which is called the Pomeranchuk 
stability condition |48|: 


(3.3.19) 


pe ss a, Pe Ses (3.3.20) 


If one of the conditions does not hold, description with the Landau Fermi liquid be- 
comes inconsistent and there are many possibilities for physical origins of these instabilities. 
Pomeranchuk originally derived these conditions by studying energy gains under the deforma- 
tion of the Fermi surface in [48]. 

In the rest of this subsection, let us concentrate on the RG flow in the BCS channel. For 
this purpose, put ky = —k, and k3 = —k/. In this case, Q = ki — kı or Q! = ki — kı need not be 
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small and then the ZS and ZS’ graphs are suppressed due to the restriction of configurations of 
loop momenta. Therefore, we neglect these graphs and only take into consideration the effect of 
the BCS graph. This approximation respects the symmetry of the vertex, and thus it makes the 
analysis much simpler than that of the Landau channel. Since we have set kı + k2 = 0, analytic 
expressions for the BCS diagram Eqs. (3.2.18]3.2.19) are simplified. The 1-loop contribution of 
the BCS graph to this channel is given by the flow equation 


—AQT? = -tanh (ved) Ner) paz (I : odd) (3.3.21) 
—AQT? = —tanh (Gerd) NE) (T?)?. (1: even) (3.3.22) 


By solving the flow equations (3.3.21} |3.3.22), we find the relation between the renormalized 
: A/S* . A/S(0), 
couplings T; and the bare couplings T; : 


A/S(0 
past _ re) 


; = f 
1+ Aapa ma tå tanh (Sur A) 


(3.3.23) 


The Landau Fermi liquid remains stable in the BCS channel at any temperature when the 
interaction in the BCS channel is marginally irrelevant, i.e., bare couplings satisfy 


TAA So. (3.3.24) 


On the other hand, if one of the conditions is violated, the coupling becomes marginally 
relevant and a Landau pole appears at some low temperature. To evaluate the critical temper- 
ature, let us consider the low temperature limit. In the low temperature limit, we can calculate 
the integration over A quite accurately: 


A 
e dA B 2eEvp Ao 
— tanh | —vgA | ~ In | ————_ ). 3.3.25 
| o (Se a( rT ) l ) 


Therefore the critical temperature is 


2eE 


1 
Ne ——— ], 3.3.26 
ao (ee) een) 


Lee 


where ©; is the maximally negative bare coupling constant which violates the condition given 
in Eq. (8.3.24). 

We should emphasize the difference between singularities in the BCS channel and Pomer- 
anchuk instabilities. In the case of Pomeranchuk instabilities, the Landau Fermi liquid theory 
itself contains self-contradictions and some of the assumptions on low energy excitations must 
be wrong. However, in the case of singularities in the BCS channel, the Landau Fermi liquid 
theory is consistent as long as temperatures are higher than the critical temperature given in 
Eq. (3.3.26). When T > T., we can find no singularities from these analyses and then it is 
natural to conclude that we can rely on the Fermi liquid description at T > T, but some phase 
transition occur at T = T, so that it breaks down below T.. 
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3.4 Open problems to be considered 


In this section, we pose some open problems based on the lessens obtained in the previous 
sections. 

As we have seen in sec[3.3] as long as we assume that low energy excitations are given by 
fermionic quasiparticles and that the effective cutoff is sufficiently lowered to justify Eq. (3.3.2), 
the Landau Fermi liquid theory is justified and many important properties of many-body physics 
can be described in a unified way. However, if we have satisfied at this stage, then we cannot 
predict many-body physics quantitatively from data of few body physics. In order to use results 
in sec 3.3] we should already know parameters of the low energy effective action. We will see 
in sec[4.5] what happens when we calculate the critical temperature using the renormalized 
coupling of scattering physics. 

In sec [3.3] we have also seen that the Pomeranchuk stability condition (3.3.20) is derived 
from the RG flow of particle-hole channels. We would also like to ask what if the coupling is 
so strong that the Pomeranchuk instability occurs. The answer would depend on situations. 
In sec{4.5] we will see that weakly interacting bosonic bound states condense at sufficiently low 
temperatures if its density is very low but the Pomeranchuk stability condition is violated due 
to the sufficiently strong attractive interaction. We should emphasize that deriving this result 
would become difficult if we have performed the Hubbard-Stratonovich transformation in the 
particle-hole channel, although the bosonization in this channel may seem to be plausible from 
observation of the Pomeranchuk instability. 

In order to clarify physical origins of these instability, it is much better to use computational 
methods which are less biased. From this viewpoint, we will use FRG for fermionic systems 
and do not rely on a bosonization such as the Hubbard-Stratonovich transformation, although 
it is often a quite powerful tool when we have already identified the origin of instabilities. 

In chap|4] we apply FRG to the BCS-BEC crossover, which is a familiar system and we 
already know its properties to some extent. However, physics of the BCS-BEC crossover cannot 
be described with the perturbation theory and we will encounter many nontrivial things as we 
shall see. 

In chap{5} we will consider dipolar fermionic systems with FRG. This system is not yet well 
studied and it is fruitful to analyze what kinds of instabilities can occur even at qualitative 
levels. 
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Chapter 4 


BCS-BEC crossover 


Let us discuss the BCS-BEC crossover in our framework. In the weakly attracting two- 
component fermionic system, the Cooper instability appears due to the fluctuations around 
the Fermi surface at low temperatures. By A. Leggett [49], and by P. Nozieres and S. Schmitt- 
Rink [50], it was suggested in 1980s that the BCS superfluidity smoothly continues to the 
BEC of bosonic molecules as the attraction increases, which is called the BCS-BEC crossover. 
Although it had been considered as a theoretical toy model for a long time, the BCS-BEC 
crossover now becomes one of the interesting and realistic physical problems after successes of 
the superfluidity of fermionic pairs in the cold atom experiments in 2004. 

In sec|4.1] we give a brief and qualitative description of the BCS-BEC crossover, and we also 
review some experimental successes of this field. Since the experimental techniques in the cold 
atoms are developing, we can expect that this provides a severe testing ground for theoretical 
physicists. 

We will use the FRG formalism developed in previous chapters to analyze the BCS-BEC 
crossover. For that purpose, we need to renormalize the bare coupling. In sec[4.2} we calculate 
scattering processes and renormalize the bare coupling ga, using the s-wave scattering length 
ag, which provides the unique length scale in the vacuum. Furthermore, we also consider the 
atom-dimer scattering problem using FRG, and explicitly show the equivalence between our 
new method and conventional ones. 

For the BCS-BEC crossover, the relevant extensive quantities are the number density n and 
the superfluid gap A. In sec[4.3] we will provide several ways to calculate the number density 
n: one is a direct calculation from the full 1PI effective action [,~9 and the another is the use 
of the flow equation of the scale dependent number density ng. Those ways are equivalent in 
full calculations, but in practical calculations the different methods may give different results 
due to approximations. 

In sec[4.4| we derive the condition for the appearance of the off-diagonal long range or- 
der (ODLRO), which gives the Thouless criterion (51). At first, we consider the consequence 
obtained by symmetries in the Wetterich formalism of the FRG in a general way using the 
Ward-Takahashi identity. We will apply it to spontaneous symmetry breaking to derive the 
Thouless criterion for superconductivity. 

In sec{4.5| we reproduce the NSR theory in the FRG formalism from a purely fermionic 
theory. We discuss to what extent we can justify that theory from the RG viewpoint, and we 
estimate the possible correction which arises in a natural way. Especially, we discuss Gorkov 
and Melik-Barkhudarov correction in the BCS regime from the viewpoint of scale separations 
and clarify its physical origin. We will see that similar ideas also give free bosonic picture in 
the BEC regime without introducing any bosonic auxiliary fields. 
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4.1 Introduction to the BCS-BEC crossover 


In this section, we give a brief and qualitative introduction to the BCS-BEC crossover, which 
is a target of our applications of the FRG method. 

We consider the two-component fermionic system with the contact interaction. That is, it 
is described by the classical action 


s= fate (Be) (@- FT -n) vw) + WHat]. 4.4) 


Here Y, w describes the two-component fermionic field. In order to circumvent possible misun- 
derstandings, we use 1, 2 as the flavor labels, which should be replaced by f, J in the comparison 
with notations in sec[3.1| Here we consider the attractive interaction (g < 0) between the dif- 
ferent fermionic species. We describe the nature of the BCS-BEC crossover before considering 
how we can realize such a system. 

At first, let us suppose that there exists the attraction but it is weak enough. In this case, 
if the temperature is not so low the Fermi liquid should be a good picture for the model, as we 
have seen in sec[3.3| However, the Fermi liquid picture is not valid at T = 0, since there is a 
phase transition to the BCS superfluid at the critical temperature T = Te. This is suggested 
by the existence of a Landau pole at temperatures lower than the critical temperature T, given 
in Eq. (3.3.26). 

Next, let us consider what happens if the attraction is very strong. It is easy to imagine 
that for such strong attractive interaction two fermions 1 and 2 form a bound state, which is 
often called a dimer. Since the contact interaction only has the s-wave component in the partial 
wave expansion, that is, its orbital angular momentum L is zero, the bound state would also 
be an s-wave pairing. When the system is cooled down, the state will show the Bose-Einstein 
condensation (BEC) of bosonic dimers 

In the seminal paper by A. Leggett, it was suggested that these two different superfluid 
states connect to each other smoothly as the coupling changes. He pointed out that the BCS 
ground state wave function 


gcs = J [0x + URC at pilva) (4.1.2) 
k 


is applicable not only in the weak coupling region originally suggested by Bardeen, Cooper, and 
Schrieffer (BCS) but also in the strong coupling region if the chemical potential u is treated 
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Figure 4.1: The schematic view of the BCS-BEC crossover. The horizontal axis represents the 
inverse of the scattering length ag, which controls the interaction between different fermion 
species. 
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Figure 4.2: A schematic of the Feshbach reso- Figure 4.3: Characteristic behavior of the scat- 
nance for ®Li. tering length ag of °Li. 


self-consistently. The theory was extended to finite-temperature treatments by Nozieres and 
Schmitt-Rink (NSR) [50], which is called the NSR theory. This smooth connection between 
the BCS state of fermionic atoms and the BEC state of bosonic molecules is the BCS-BEC 
crossover (see Fig/4-1). 

Now we discuss the experimental access of the BCS-BEC crossover. From now on, the 
unique possibility of the crossover is provided by ultracold atomic gas experiments of fermionic 
atoms. By cooling down the system (typically about 100nK), energies of particles become 
sufficiently low so that the s-wave scattering becomes dominant. Effectively, such interactions 
can be treated as contact interactions, and due to the Pauli exclusion principle the scattering 
between the same species is forbidden. Remarkably, in these systems the interaction between 
different species can be tuned with the external magnetic field via the Feshbach resonances. 
We will see some details below. 

Let us explain about the Feshbach resonance [52]. To be specific, we consider the case 
of Li, which is an alkaline atom, and Figf.2] gives a schematic view of the inter-particle 
potential. Roughly speaking, each curve is described by the Lennard-Jones potential. With the 
magnetic field in the broad Feshbach resonance region, we may concentrate on the two states 
characterized by |F,mpr). Here F represents the total spin and mp represents its projection 
along the magnetic field. The two states are given by |1) := |5,$) and |2) := |3,—3). If 
we alternatively use the electronic and nuclear spins as labels, the states |1) and |2} have the 
electronic spin -4 and have the nuclear spin 1 and 0, respectively. 

S= lin Figf4.2] represents an electronic spin-triplet state, and it is called an open channel. 
Similarly, S = 0 represents an electronic spin-singlet state. Each curve represent the interaction 
for the corresponding channel. The electronic Zeeman coupling to the triplet state enables us 
to tune the relative position of the curves, described by 6(B) in Fig/4.2| with magnetic field. As 
a result, the energy level of the singlet bound state can be tuned around that of the scattering 
state of the open channel continuously. When the detuning v = 0, the scattering length 
ag diverges and this is the Feshbach resonance. Since the ultracold atomic clouds are dilute 
enough, the density, the temperature, and the scattering length are the only scale parameters. 
Therefore, after the normalization of the thermodynamic quantities with the density, they only 
depends on (krag)~' and T/ep. At the resonance ag = 00, any scale parameters except for the 
temperatures does not affect the thermodynamic properties, and this region is called unitarity. 
The characteristic behavior of the scattering length ag can be found in Fig/4.3| 

Let us comment on experimental developments of this area. In 1999, JILA group 
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Figure 4.4: “°K Cooper pair condensate [53]. 

These images are taken after the projection of Figure 4.5: °Li Cooper pair condensate |54|. 
the fermionic system to a molecular gas with The initial temperatures of the atomic clouds 
the detuning AB = 0.12, 0.25, and 0.55G (left are T/Tr = 0.2, 0.1, and 0.05 (left to right) in 
to right) on BCS side. the BCS side (B = 900G). 


succeeded to achieve quantum degeneracy of the Fermi gas (K). In 2003, some groups (JILA, 
MIT, etc. ) succeeded to realize the BEC of dimers using different elements (°Li in 
MIT, K in JILA). The observations of the BCS superfluid of these gases are done in 2004 by 
the same groups [53154]. In Figs/4.4| and we show the density profile of these condensate 
of the Cooper pairs in the momentum space. 


4.2 Scattering problems in the vacuum 


Our final goal is to describe many-body physics from the knowledge of few body physics. For 
that purpose, we need to calculate the flow equation in the vacuum to calculate scattering 
properties in the formulation of FRG. In the followings, we often use the unit 2m = 1 for 
simplicity of calculations. 

In order to calculate the scattering physics in the vacuum, we should take the zero density 
and the zero temperature limit. This limit should be taken with some care. We should take 
two limits n + 0 and T —> 0 satisfying that the temperature T is always above the critical 
temperature Tę. For that purpose, recall that the partition function Z is given by 


Z=Tr [exp —8 (4 — pn)! ; (4.2.1) 


where H is the Hamiltonian and N is the number operator. As long as the ground state energy 

per particle is bounded from below, we can make it by taking sufficiently large negative values 

of the chemical potential u and T — 0. If there is no bound state we can simply take u = 0. If 

the bound state appears, the binding energy Ep of the bound state can produce a pole in the 

vertex function to give a finite density and we should let 2u < — F, to circumvent this problem. 
Let us consider the self-energy correction. We obtain that 


ox [fl fa T(p+de#O 
kdr (p) = ôr J EAE J E + Ral) 2-2) 


We can close the contour of the 1° integration encircling the upper-half plane, and then the 1° 
integration vanishes since the integrand has poles only in the lower-half plane. As a result we 
find that 

z(p) = 0 (4.2.3) 


in the vacuum. 
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4.2.1 'Two-body scattering 


Next, let us consider the renormalization of the four-point vertex. Since the feedback from 
the higher order vertex and the particle-hole loop vanish in a similar way, the only possible 
quantum correction come from the particle-particle loop. The loop correction of four-point 
vertices is given by 

1’ 2’ 


_ Tris (P: g i= Pie (py = È q’) 
9.04/5(P: 9, q) = — a i is f BN ae agd 
meee key PSA) (G+ RIG + RIP — oy OP") 


At the bare scale, iy = 0 and ry. = ga,- Therefore, at any scale the anti-symmetric part 
vanishes; T4(P; q, q’) = 0. Even for the symmetric part T$, the relative momentum dependence 
does not appear, since at the bare scale rX = gẹ is a constant. Therefore, we may put ['?(P) = 
['?(P;¢,q). Now we have found that the vertex functions T? in the right hand side of Eq. (4.2.4) 
can be put outside of the integration, we can readily get 


1 -a i , 
OTE) ~ pee i | (Got + Ra] (G7! + Bal(P — I) (4.2.5) 


Since the k-dependence does appear only through R, in the right hand side, On coincides with 
ôk. Then we obtain that 


1 1 al 1 
T'e(P) p Ino 7 | (2Q7)4 (aes + RJD) (Gs! + Ryl(P —D) (4.2.6) 


1 
~ (G51 + Ra ]OMIGs? + Ral = y) 


Let us choose the regulator given by R,(l) = (k? — 1?)0(k? — I’), which is a nonrelativistic 
version of the Litim regulator discussed in sec [2.4] Recall that we use the unit 2m = 1, then 
the integration over 1° gives 


1 1 al 1 
TR (P) 7 GN. 7 / (27)? (= — 2p +P + R0) + (P -1 + R(P —) (4.2.7) 


1 
{po — 2u +P + Ry (l) + (P—-1?4+ Ry, (P — 5) 
In order to fix the renormalization condition, let us put P = 0 and iP® = 2p in Eq. (4.2.7). It 
is easy to calculate the integration analytically, and then we get the result 


1 1 1 
— = Ao — k). 4.2.8 
T$ (iP? = 2u, P = 0) Oro 572 | ) ( ) 


At k = 0, T'2_,(P° = i(2q? — 2), P = 0) represents the T-matrix for the scattering between 
two different fermion species in the center of mass system. Therefore, it should be related to 


the scattering length ag in the case P? = —2iu and P => 0. 


2 4 
The T-matrix and the scattering length a, is related by lim T(q’,q) = BOE e Gog 
|q|0 m m 


r 


the particles with different species |*| Therefore, we should fine tune the bare coupling ga, of 


'In general, we can show that in the center of mass system 
m 

f(q, q’) = =a E q) 

T 
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the contact interaction so that 
ee ee doa 
ga,  8ras 6r? = 


By the dimensional analysis, we can recover the formula in the original unit if necessary: 
1 m Ao 


Jho ~ 4nag 6m? 
Let us explicitly calculate the P-dependence of the 1PI four-point vertex ['?_,(P) of the 
singlet channel. In the limit A, — oo we get 


3 1 1 
Ce f N n ae F (4.2.10) 
Peg(P). 8ras (27)? | iP? — 2u + (lye Ge Hl W 
We can explicitly calculate this integration without any divergence to find that 
S(p0 p2 87 
rS(P?, P?) = (4.2.11) 


1/as — (iP? — 2u + P2/2)/2 


We should notice that this expression satisfies the Galilean invariance since the center of mass 
momentum appears only through the combination iP? + P?/4m. We can find calculations with 
general cutoff scales k in Appendix[D| There we also discuss how to remove the UV cutoff scale 
Nos 

This is obtained in the imaginary time formalism, and then we need to perform the Wick 
rotation P? +2iu = iE +0 in order to get a result in the real time formalism. In this scattering 
process, we need to input the total energy EF = 22 = q’/m in the center of mass system 
P =0, and then the scattering amplitude is given by 


f __ mMm 8T = 1 
fla q) E AT 2m/as = /—(2m)3(q?2/m) /2 ilq| = l/ag (4.2.12) 


This is nothing but the effective range expansion of the contact interaction. Since we have 
taken the limit A, — oo, the effective range becomes zero. 

When the scattering length is negative, ag < 0, then zeros do not appear in the real axis of 
the energy E. Then, there is no bound state and we can take the limit «7 0 as a result. On 
the other hand, if ag > 0, then the pole of T’ appears with the dispersion relation 


P 2 P ıl 
2 a% 4m mag 


(4.2.13) 


This implies that there exists a bound state in the spin singlet channel, and its binding energy 
EF, is given by 1/maz. This is a well-known result for the two-component fermionic system 
with the contact interaction. We should restrict the chemical potential u < —E,/2 = —1/a% 


for the zero density condition. 


for the non-relativistic elastic scattering. Here q’ = |q|r and my is the reduced mass. We can find the proof 
of this relation in some textbooks on the scattering theory, e.g., in chap.10 of the book [58]. We can find the 
result for the general situations in chap.3 of the textbook [1]. 
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4.2.2 Three-body sector: atom-dimer scattering 


From previous calculations, in the case ag — —0 the notion of weakly attracting fermions is 
a good starting point and the Landau Fermi liquid theory gives a good description of many 
body physics. On the other hand, in the limit ag — +0 a bound state, or a dimer, appears and 
fermions are gapped by half of the binding energy E,/2 = 1/a%. Therefore, we have to derive 
the bosonic picture from our fermionic theory in order to describe low-energy properties in this 
regime. For this purpose, the calculation for the dimer-dimer scattering is important. 

As suggested by S. Weinberg in his work on composite particles and on multi- 
particle scattering problems, the vertex for the dimer-dimer scattering should be represented 
with a vertex for the atom-dimer scattering. We can easily see this fact from the viewpoint of 
the Wetterich equation, and we at first consider the scattering problem between an atom and 
a dimer. 

In the case ag < 0, the effective action up to the quartic terms is given by 


T = f Pp(ip? + p° — utp (4.2.14) 


ah T 8T 
+f (f Prerna) 1/as — VUP? — Ou + P?/2)/2 (/ Yupapa | ; 


The pole of the four-point vertex is given by the dispersion relation —iP° = P? /2 —2/az — 2p, 
and then we should take u < —1/a2 in order for the zero density condition. We can take the 
limit u Z —1/az to get 


i 16r /as (1 + 1+ uP + P?/2)) 
T= |y, (iP +p +> -]/ S p. (4.2.1 


Here we have introduced the following short-hand notation for the simplicity of the expressions: 


bp = f E EA Dee / TEEN ae (4.2.16) 
q q 


We should emphasize that Eq. (4.2.16) does not represent any transformations such as the 
Hubbard-Stratonovich transformation and that it is just a definition of the convenient notation 
for the composite operators. From Eq. (4.2.15), we can identify the boson propagator by iden- 
tifying the pole which is necessary for the Lehmann-Symanzik-Zimmermann (LSZ) reduction 
formula in calculating scattering processes. That is, we can identify the fermion four-point 


vertex as 
. = xK l (4.2.17) 


where in the left hand side the arrowed double line represents the dimer propagator D(P) = 
—(1+ y1 + È (į p2 + F*))/(2iP° + P?) and the vertex represents \/327/ag = \/8m/m?as. 
Now let us consider the scattering between an atom Î and a spin-singlet bound state. We 


would like to determine the scattering length aaa of the atom-dimer scattering. For the atom- 
dimer scattering, contributing Feynman diagrams consist of tree and 1PI diagrams. The tree 
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diagram for the atom-dimer scattering T-matrix is represented as 


4/32 z 8 
Tae = = l . an) a) z = lOnas = ll (4.2.18) 
i(—t/a%) + 1/az m 


Here we should notice that momenta of external legs are given by p = p' = (i/a, 0) for external 
legs of a fermion and P = P’ = 0 for external legs of a dimer. Thus, the transfer momentum 
of the internal fermion line is given by q = (—i/a%,0). Since the T-matrix is related to the 
atom-dimer scattering length aaqa through the formula T°! = 27aqq/m, = 3Taaa/M, we get 


8 
ac? Jas = 35 2.67 (4.2.19) 


as a tree diagram contribution. 
Next, let us evaluate contribution of 1PI diagrams to the atom-dimer scattering. Diagram- 
matically, the flow equation can be represented as follows: 


on = dp + Se + X + P (4.2.20) 


Here we should regard that external legs of each diagram are antisymmetrized. The first 
diagram may seem to vanish at first sight since it represents a feedback from higher vertices, 
but it can contribute due to non-locality of effective vertices} Analytic expression of Eq. (4.2.20) 
is given by 


-3 Tk (p, P2, P3; P3, Pa P1) 


= f De Pi, p2, l, p3; P3, l, Ph, p1) Pe pupas L, ps, P,P) 
[Got + Rel) Go! + Rel) 


(4.2.21) 


45, f T$ (p243) (p1, l, poss — l; P3, Po, P1) 4 T$ (pi43) Pelt (1, p2, Piss — l; p3, P2, P1) 
t JA [Gs FRIOG + Rilos) | [GF Rel [Gs + Rips — 1 
Pin (D1, Pa, Psi Pha — LA PAE Oa) Taed Pigs — Pa W 
[Got + Ral (Q[Go* + Rel (Po43 — 9) [Got + Ral (Q[Go* + Rel (piss — 1) 


2 1 
4s | eID 


x (- De (p2+3)TR (Phs) R (P1 +1) a De (Pisa) R (Po43) ER (P2 + 1) 
[Got + Ril + piv )[Go' + Re|(po43— 1) [Got + Ri](l + po-v)[Got + Re] (pi4s — !) 
a: De (p2+3) E (Pas) Pe (p: +1) 7 De (pita UE (Piya) Pe (p2 +1) ) 
[Got + Rall + pi-a [Got + Ral (pers — 1) [G71 + Re](0+ p22) [Got + Re] (Piss — l) 


where we have introduced the abbreviation p213 = p2 + p3, Pi_v = Pi — p4, etc. 


?Let us emphasize that this fact is consistent with a general property in nonrelativistic scattering physics: 
We can solve problems for n-body scattering without knowing about (n + 1)-body scattering problems. As we 
will see later, the feedback from 1PI eight-point vertices can be written with vertices in three-body sector. Since 
this point is often missed, this feedback is neglected in many previous studies. 
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We would like to decompose it corresponding to the diagram 


K = K (4.2.22) 


Therefore, we define a new vertex function g motivated by Eq.(4.2.22) so that 


Penia P's Pa, Ph) (4.2.23) 
= —T? (po+s) 96" (p1, P2+3; P243: PLIE (Po43) T T? (piss) 98" (p2, P1433 Pasg: POLE (Po+3) 
+r? (pois) 9¢"(p1, P2+3; Pijs w Pus) = TÈ (piss) 90" (pe, P14+3; P43) paT iss): 


Indeed, we can show that this decomposition is possible directly from a formal diagrammatic 
expansion of six-point vertex in terms of the coupling ag. With this decomposition, we can 
relatively easily represent the three-body sector in eight-point vertices, which can contribute 
the first diagram in the right hand side of Eq. (4.2.20). In eight-point vertices, the three-body 
sector is composed of the following diagrams: 


Now we can immediately obtain an expression for the feedback term in the flow equation of 
vertex functions gł. Amputate the bosonic external lines in each graph, and close two fermion 
lines with the two-point function 0;R, so as not to make particle-hole loops. This procedure 
is clearly possible in each graph in Eq. (4.2.24), and then we here explicitly show that feedback 
from higher-order vertices is possible even in the vacuum, but we also find that higher-body 
sectors do not come in for nonrelativistic scatterings. 

In order to obtain other part of the flow equation of 0,g%4, let us substitute Eq. (4 into 


Eq. (4.2.21) and use of the flow equation for [?. Besides the above eight terms coming Ai the 
feedback, we have another three terms for the flow equation of g4: 


3-body 


Onge (p1, P243; Pie, p,) = 8 feedback terms 


15, [Aes — Doh, p1+2+3 — l; poss P1) | 9R (P1, D243; Pryoes — LOTR (Pit243 — D) 
IA [Got + Ri] (D[G + Ri] (pots — 1) Got + Rel(O(GS* + Rel(P243 — 1) 


r? (pı +1) 
"Te PRMIG hae A Ga = 5) . (4.2.25) 


Diagrammatically, Eq.(4.2.25) can be represented as 


~ Ý 
On 4 = ók S + a, + yí + 8 feedback terms, (4.2.26) 
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where “eight feedback terms” are immediately obtainable from Eq. (4.2.24) as we have men- 

tioned. Precisely speaking, we should multiply a factor —327/ag to the both sides of Eq. 
to get Eq. (4.2.26). The function gz is related to the 1PI contribution to the atom-dimer scat- 
tering T-matrix, and the LSZ reduction formula gives 


6 327 a 
TH = Ge Io ((i/a%, 0), 0; 0, (i/a3,0)) . (4.2.27) 


Let us see some properties of this new vertex function g%%. The restriction coming from 
Fermi statistics is automatically respected in Eq.(4.2.23), and then it does not restrict the form 
of g@7. Using the Euclidean time reflection invariance, we find that 


I (Pr, P2433 P3243 pi) = at" (Opi, 015433 Op2+3, Op1)]". (4.2.28) 


Next, let us reduce the number of variables in gł. Due to the momentum conservation, 
P := pı + poss = p + Po43- Furthermore, we can easily find that we may put P = (i/a%,0) in 


Eq. (4.2.25) to extract the quantity in Eq. (4.2.27). Introducing 
lnd) = 98(P +9,-G-d,P +4), (4.2.29) 


with P = (i/a2,0), we can rewrite Eq.(4.2.25), or Eq.(4.2.26), as 


ai Ao re (1) a a) g(a.) 
ostlai | CAES AET (a ER sPo0 2s CoP oe) 
1 
+ [G51 + R,|(—l — P — q)[G7} + R,|(-l- P — =) + feedback. (4.2.30) 


In Appendix|D| we solve this equation with an approximation in which we neglect the feedback 
effect from eight-point vertices. In the following, instead of solving this flow equation we discuss 
its properties and find formal solutions using Dyson-Schwinger equations. 


Formal solutions 


Let us discuss a formal solution of Pen , or głf. As we will see below using Dyson-Schwinger 
equations, the vertex function g?7(q, q’) satisfies the integral equation 


a e rg(—0) IL) 
Ik (%4) = le + Ry](l+ P) (a Fir — P= @) 


1 
JC + R,)(—l — P — 9 [G51 + R,|(-1- P- 5) (4.2.31) 


and similarly satisfies 


eh ré(-) g (a,l) 
IGT) = he + Ril(0+ P) (a + Ral(-1- P- q') 


1 
Ñ [G51 + R,](-1 — P — q)[Gs1+R,](-1- P — z) . (4.2.32) 


At k = 0, these integral equations reduce to the integral equation for the atom-dimer scattering 
T-matrix, which was derived by G. Skornyakov and K. Ter-Martirosyan [63]. 
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Using the Dyson-Schwinger equation (2.2.26), or oh (2.2.302.2.31), we can derive inte- 


gral equations given above. Let us recall feet EIA suggests that the Dyson-Schwinger 
equation must be true even for IR regularized E e “ho the six-point vertex [k P in our 
model, we have two Dyson-Schwinger equations in the vacuum. Diagrammatically, they can be 


represented as 


where black dots represent the bare contact interaction, and as 


KERK = 


In deriving the second Dyson-Schwinger equation, Eq.(4.2.34), for the six-point vertex, we 
use the Dyson-Schwinger equation for the four-point a ae the first Dyson-Schwinger 
equation, Eq. (4.2.33), for the six-point vertex. This second Sen (4.2.34) of the Dyson- 
Schwinger equation is equivalent to the integral equation for three-body scattering problems 
derived by S. Weinberg in (61). 

Combining these two Dyson-Schwinger equations and applying the decomposition (4.2.23) 
of the six-point vertex, we can obtain an integral equation. Its diagrammatic expression is 


given by 
kd pf nits 


This is nothing but the integral equation (4.2.31). One can find direct ooh ae of this integral 

equation by considering all possible Feynman diagrams in the paper We here invented 
a new way to derive the integral equation for the three-body co ering T-matrix by 
using FRG, and explicitly proved that those two formalism can give the same answer. 

Let us explicitly show that the solution of this integral equation satisfies the flow equation 
(4.2.30). For this purpose, it is better to write down the integral equation (4.2.35) symbolically 
as 

g@=G-K-g%4+G-K-G, (4.2.36) 


me ost «= t (4.2.37) 


Here we made their k-dependence just for simplicity of notations. We can formally write down 


the solution of Eq.(4.2.36): 


g@=(1-G-K)'-G-K-G=G-K-G4+G-K-G-K-G+---. (4.2.38) 


where 


In order to derive the flow equation, let us take a derivative with respect to k of the both sides 
of Eq. (4.2.38). Leibniz rule and formal resummation gives 
0.9% = 0,(G-K-G)+O(G-K)-(G-K-G+---)+(G-K-G+---)-0,(K-G) (42.39) 
+(G+G-K-G+.--)-K-OG-K-(G+G-K-G+---) 
+(G-K-G+G-K-G-K.-G+.::-)-0kK-(G-K-G+G-K-G-K-G+--:-). 
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Combining Eqs. (4.2.38]4.2.39), we obtain 


kg = O,(G-K-G)+0,(G-K)-9%4+9%7-0,(K-G)+(G+9%)-K-O.G-K-(G+g%)+9%4-O,K «9%. 

(4.2.40) 
Now it is clear that Eqs. (4.2.24]4.2.26) and Eq. (4.2.40) are equivalent as it should be since 
Dyson-Schwinger equation must hold at any values of k. 

It would be very important to comment that, in the flow equation of gt, quadratic terms 
of the vertex g?" appears, and that it originates from the feedback from eight-point vertices in 
the flow of six-point vertices. This structure must become important when we discuss Efimov 
physics without introducing auxiliary composite fields, since limit cycle structure can occur 
under existence of this quadratic term. 


Remark. We can expect that a formulation for the dimer-dimer scattering problems is also 
given in a similar way, but this task is still open in our formalism. Even though we can solve 
such problems by deriving integral equations directly from Dyson-Schwinger equation, as we 
have done here for the atom-dimer scattering, formulating that problem in the context of FRG 
is still important. Using FRG, we can study many-body physics by reflecting knowledge of 
few-body physics thanks to ideas of scale separations. Therefore, formulating few-body physics 
in FRG and solving it, we can immediately find some knowledge of many-body physics. 


4.3 Number equation 


In specifying thermodynamic states, it is natural to fix the number of fermions Np rather than 
the chemical potential u especially in the strong coupling region. Of course, those two quantities 
are related via the Legendre transformation and then they have equivalent informations on a 
system in principle. However, in the strong coupling region we often need a crude approximation 
for calculations and lose the correct relation between them. In order to obtain sensible results 
even in such situations, we should treat relevant operators themselves, and in this section we 
will explicitly show why the NSR theory gives a correct picture of the BCS-BEC crossover from 
the viewpoint of the composite operator flow in sec [2.2 
Let us define the number of fermions ng in IR regularized theories by 


In stead of using this formula for the number density n = ng=0, we will try to use the composite 
operator flow equation. The number density operator is given by 


nip, y] = Y(T +0, x) (7, £), (4.3.1) 


and then we can extend it as a correlation function n,[v), Y] according to Eq. (2.2.18). It obeys 
the flow equation (2.2.19). Let us denote the expansion of the correlation function n,|q, Y] in 
terms of fields Y, Y as 


At k = Ao, na, = 0 and ng” = ( since they are absent in the classical limit, and n” = 
O30" (x — y) in the coordinate basis. At k = 0, nk=o should give the number density of the 
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full theory n. Let us introduce a diagrammatic expression for the expansion of Eq. (4.3.2) as 
follows: Y w 


niy, y]= A + y eae y+ Ww free, (4.3.3) 


Y Y 
We can obtain the flow equation for each coefficient. Let us write it down using Eq. (2.2.19) 
and using the technique described in Appendix [A] The flow equations for the constant and 
two-point parts are given by 


a a Cy E Cy asa 
ðk | y >A y| = (R 0) (4.3.5) 


Triangle vertices represent operators in Eq. (4.3.3 , square ones represent the negative of 1PI 
vertices, and the blob represent the derivative a the IR regulator 0, Ry. Lines with an arrow 
represent propagators (Gj! — ©; + R,)~!, and we should regard that all external lines are 
amputated in the both sides of equations. 

On the other hand, we can calculate each coefficient in ngi, y] directly using formulae 

(A.2.3]A.2.5). Let us denote the insertion of the two-point vertex y(x) with the crossed 


vertex. Then, we get 
Je-il?o 
4.3.6 
-f ica ESTA ( ) 


y =A y = (2) + . (4.3.7) 


For later convenience, we give the analytic expression of Eq. (4.3.6). In the numerator, the 
negative sign appears due to the fermion loop, the factor 2 appears due to the number of spin 
degrees of freedom, and the factor e™#°? appears due to the infinitesimal difference of imaginary 
times of field arguments in the sal} of Eq.(4.3.1 (4.3.1). 

Substituting the result of Eq.(4.3.7) into the right hand side of Eq. (4.3.4), we can see 
consistency between the flow oo and its solution. For this purpose, let us substitute 
Eq. (4.3.7) in to the flow equation of ny to get 


Recall that the flow equation for the self-energy is given by Eq.(3.2.6), which can be diagram- 
matically represented as 1. B 


A, D42(p) = de , © mi (4.3.9) 


p p 
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Combining Eqs. ((4.3.8]4.3.9), we obtain that 


(T) —2e7 i0 
OkNk = f (en Zig AOA eRe) + kkl) = On C) . (4.3.10) 


As a result, we find the consistency between the flow equation of the correlation function 
ng, y] and its analytical expression with 1PI vertices. 


As a by-product of this consideration, we have now obtained two n ways in order for 
evaluations of the number density. One way is to use the formula (4 at k = 0, then we can 
obtain the number density directly: 


T) _9 —il?0 
n= Neo = f oa a (4.3.11) 
l 


where X = X}=0. As long as we have solved the flow equation of the self-energy, we can obtain 
the number density using it. 

Another way to get the number density is to use the flow equation (4.3.8). The analytic 
expression is given by 


Da _ f Je #8, Ri, (1) P T De ilo +1'°)0 ‘(ary r$] (l Ll : n EE, LL) ra) k RCl”) 
aO d e O — Ee + RO)? 
(4.3.12) 
If we solve the flow equation containing up to four-point vertices, we can immediately evaluate 
kng. Therefore, we can also evaluate n = no using this flow equation, which does not require 
to solve another renormalization group equation. 

Of course, Eqs. (4.3.11]4.3.12) are equivalent, and they should give the same answer. How- 
ever, in practical calculations, we have to adopt some approximation schemes in solving the 
Wetterich equation. Approximations break this equivalence and they will give different results. 
For example, in the NSR theory the self-energy correction is totally neglected and then if we 
use Eq. (4.3.11) 4.3.11) the number og EST is ee that of free fermions. That must be totally wrong at 
least in the BEC regime. In Eq. (4.3.12), on the other hand, the second term in the right hand 
side contains the effect of the Peer of dimers through the four-point vertex, and it would 
give a better solution when X% is put to be 0 due to approximations. 


4.4 Thouless criterion 


In order to realize good approximations to the Wetterich equation, we need to have a physical 
insight to the system which we are studying. Especially if the system has a symmetry then we 
can use it as a guiding principle and we would like to respect the symmetry in the flow of the 
effective action. In this section, we would like to study the effect of symmetries in the FRG 
formalism. 

We also consider a second order phase transition, and especially we will derive the Thouless 
criterion for the superfluid phase transition, which plays an essential role in the determination 
of the critical temperature. 
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4.4.1 Ward-Takahashi identity 


At first, let us briefly derive the Ward-Takahashi (WT) identity for the full theory and 
by finding the correlation functional Jwr|J, ¢]. Let us assume that the theory is invariant under 
an infinitesimal transformation of fields ¢, + ¢, +¢®,[@], i.e., the product of the measure Do 
and the Boltzmann factor exp (—S[@]) does not change under the transformation. Here @,,[¢] 
is a local functional of ¢. Change path integral variables Øn 4 Øn + EPn in Eq. (2.1.4), then we 
obtain 


0= f Dé [exp tsis reta en (—Sl4] + a] 
= eJ" | Dé $,[6]exp(-S|6] + 7"). 


Therefore, this is the WT identity and we can rewrite it as 


ô 
I"P,, H exp (W[J]) = 0. (4.4.1) 
Performing the Legendre transformation and substituting SE = Yn, we can find that 
ôrT [e] ÔL 
o G, Gk- 1 \=0, 4.4.2 
Ce (aa e+ Gla) Se (4.4.2) 


where G[y] = W)[J|y]]. Especially when we are considering a symmetry generated by a linear 
transformation, we can simply put (Baly + G - (6, /dy)|1) = nly], and then the 1PI effective 
action turns out to have the same symmetry with the classical action. 

An immediate application of the above consideration shows the following: Suppose that the 
classical action S/¢] and the flat path integral measure D@ are invariant under the symmetry 
transformation independently. Then, as long as one uses the regulator Rẹ which preserves the 
same global symmetries, the effective action I% at all scales k also satisfies those symmetries. 

It is usually a difficult problem to find a symmetric regulator for a nonlinear symmetry. 
Investigating the situation where the symmetry is explicitly broken by the regulator is thus 
important, and indeed we will encounter such situations when the Nambu-Goldstone bosons are 
treated as elementary fields and when the non-Abelian gauge theories are considered. According 
to Eq. (4.4.1), the correlation function representing the WT identity is 


IwrlJ, ¢] = J” Balo]. (4.4.3) 


In order to apply Wetterich’s formulation, we define the flowing correlation function Iwr|J; Ry] 
by 


16 ô ô 
W Ro fJ: Rp] = Ma a 7 reed N, mas 4.4.4 
e wrlJ; Rx] exp ( 25J Rk z) ( WT J, 5J exp(W|[J]) ‘ ( io. ) 
Theorem suggests that Iwr[J; Rx] = 0 at any k, which is already trivial from the expression 
of Iwr[J, Ry] in Eq.(4.4.4) and from Eq. (4.4.1). Explicit calculations show that 


_ [Dd (I"Pnld] = on RE" Pm g]) exp (—Sl9] - 50: Re:e + J -9) 


Fel J; Re] = J Doexp (—S[6] — 46- Ry- 6+ J- 4) 
= Mil (ro, il SAE n, 541) Wil (4.4.5) 


= (J™ — vnl J] RE") (on fet + z] i) =R" (on et + z] i) 
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Let us perform the Legendre transformation in the sense of Eq. (2 , and then we obtain the 
flowing correlation function 


Iwely; Ri] = Iwe[J[y]: Rel (4.4.6) 
=: (on [e outl: f])) mrs (hen [pte J1) 


As we have seen in the case of the Dyson-Schwinger equation, we can regard the WT identity 
Iwr = 0 as a consequence of the integrated flow. Let us consider the classical limit k > Ag, 
then the regulator becomes heavy Rg —> oo. Since Gle] = r? lp] + Rk)™t, this limit implies 
that Rg -Gk —> 1 and Gk — 0. We have obtained that 


ôrS [e] 


i ò 
jim Iwrly; Re] = om Palp] — ST (= 


Pm 


20) . (4.4.7) 


In Eq. (4.4.7), the first term in the right hand side gives the infinitesimal change of the classical 
action S|{¢], and the second term represents the Jacobian of the path integral measure D¢. 
Indeed, under the transformation ¢, 4 n + €®,|¢] the path integral measure becomes 


Do = | [den | [dom (æ + 7+ onlol) = DdSDet (æ + eia, a): (4.4.8) 


Up to the first order of the infinitesimal parameter £, the Jacobian becomes 


ÔL 
SDet (æ + eg nll) = exp €STr (he nll). (4.4.9) 
and the second term in Eq.(4.4.7) appears. Since the product Do exp(—S[¢]) is invariant under 
the symmetry transformation, the right hand side of Eq.(4.4.7) should vanish by hypothesis. 
Hence, we get the unique solution of the flow equation (2.2.19) so that Iwr[{y; Ry] = 0 for any 
k, especially for k = 0. 


Remark. With the help of the Batalin-Vilkovisky formalism, or the anti-field formalism, one 
tries to control the explicit violation of gauge symmetries due to IR regulators. The studies 
in this direction is reviewed in detail in the review (67, which is based on the Polchinski 
equation not on the Wetterich formalism which we use here. One = also studied FRG in the 
background field gauge, and for example we can find the papers (26}/68| 69|. In these studies, 
the modified version of the Slavnov-Taylor identity, which extends Eq. (4.4.6) so that it contains 
the BRST transformed fields, often plays an important role. On the other hand, there are other 
approaches to invent a manifestly gauge invariant formalism of FRG. For example, see 
for the approach to construct the manifestly gauge invariant effective action using the Wilson 
lines. Any programs in these studies are not completed, and further investigations are still 
necessary. 


4.4.2 Spontaneous symmetry breaking 


Let us discuss the spontaneous symmetry breaking (SSB) in the context of FRG. For our 
purpose, it is enough to discuss SSB associated with the second order phase transition. 
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Figure 4.6: A space of theories with a second order phase transition. The red curve represents 
a line of theories with different temperatures. We indicate two RG flows, one of which shows 
the flow for a symmetric theory and another one represents the flow for a symmetry broken 
theory. 


To solve degeneracies of vacua, we add a symmetry breaking term O[¢] to the bare action 
S|] at the UV scale k = A,. In this case, we also define the flowing Schwinger functional 
W,,| J. j A Nol by 


explo) = f Doexp (-stel — Ay,O[6] — 50° Ry O+ T- 6) | (4.4.10) 


where Aq, is a new coupling constant associated with the symmetry breaking term O. The flow- 
ing 1PI effective action T;,[y, Aa,] is again defined by the Legendre transformation of Wz[J, Aa,| 
in the sense of Eq. (2.2.3). Consider the infinitesimal transformation bn > n + Enl] which 
preserves Dé exp(—.S|@]), then we find that 


Dé [(J —d- Ry) - ldllexp (—S[d] — An, Old] — t0- Re G+ I-¢ 
2 


=o | po (FIE. oio] exp (Sie) - AsO -30 Redo). ALn 


a) 
The left hand side of Eq.(4.4.11) is nothing but the WT operator Iwr|J, Aq,; Rg] defined in 
Eq. (4.4.4). Let us denote that dO[¢] = we - b|¢], then the right hand side of Eq. (4.4.11) 


becomes gdO;|J, g]. Therefore, the WT identity should be modified as 
Fwrly, An; Re] = g5Ouly, An). (4.4.12) 


Let k = 0, then the discussion reduces to the one for the original theory. If the symmetry is 
not broken, all physical quantities are analytic in terms of the coupling g and the WT identity 
holds in the limit Aa, — 0. On the other hand, if the theory is in the symmetry broken 
phase, the right hand side of Eq. (4.4.12) does not vanish at k = 0 and in the limit Ay, > 0 
and the WT identity does not hold. 

Since we assume the second order phase transition, the critical fixed point exists. To be 
specific, let us consider a line of theories parametrized by temperatures T which is drawn as a 


68 Chapter 4. BCS-BEC crossover 


red curve in Fig/4.6] At T = T,, the line crosses the critical surface and the second order phase 
transition occurs. Let us consider the RG flow in a symmetry broken phase. Here we assume 
that the symmetry is not explicitly broken through the regulating term. At k = Aj, the theory 
is classical and the symmetry is not broken and then there exists a critical scale k = ke, at 
which the IR regularized theory lies on the critical surface [P] For k > ke the RG flow coincides 
with the symmetric one in the limit Ay, — 0, and for k < ke flowing 1PI effective action Tẹ 
does not satisfy the symmetry. 

In order to get the Thouless criterion for the BCS-BEC crossover, let us consider fermionic 
theories with the singlet pairing. Since our regulator YR manifestly satisfies the global U(1) 


symmetry, Eq. (4.4.12) becomes 
STA, Y] = gO, 4. (4.4.13) 


Our purpose is not to derive the gap equation but to get the critical temperature, and thus we 
can neglect anomalous vertices with legs > 4 because they also vanish in the limit k > ke or 
T > T, with An, > 0. 

As a symmetry breaking term, we take a fermion bilinear 


O= f des (W(x)e(x) — B(a)eB(a)) . (4.4.14) 


Within the approximation neglecting higher order anomalous vertices, we find that Eq. (4.4.13) 
becomes 


A; := >| = Ay, | >e + ; (4.4.15) 


If the original theory shows the spontaneous symmetry breaking, scale dependent gaps A; does 

not vanish in the limit Aa, — 0 for k < ke. Therefore, at k = ke the four point vertex in the 
singlet channel r3 (P) diverges if we put P = 0. Especially at the critical temperature T = Te, 
the critical scale ke = 0, and then we obtain the Thouless criterion for superconductivity (51): 
(PP) > œ as P 0. 


Remark. We should emphasize that Eq. (4.4.15) is a result of approximations, which are valid 
only when Ay, > 0 and k Z ke. In order to calculate the correct gap in the symmetry 
broken phase, we have to invent an approximation which is compatible with the WT identity 
Eq. (4.4.12) or Eq. (4.4.13). Even though our model only has a simple symmetry, i.e., the global 
U(1) symmetry, it is still nontrivial how to keep the constraint in the RG flow of symmetry 
broken phases. In this thesis we do not attack this problem, but since we have derived the 
Thouless criterion for the superfluid phase transition we can calculate the critical temperature 
of our model. 


4.5 Nozieres-Schmitt-Rink theory 


In this section, we try to reproduce the NSR theory from the viewpoint of FRG. Since it is 
valid only in the deep BCS and in the deep BEC regimes, we consider those two limits and 


3Here we implicitly assume that any exotic behaviors of the RG flow is absent: the RG flow crosses the 
critical surface only once. In our application of FRG to the BCS-BEC crossover, this assumption is plausible. 
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simply connect the results. At the same time, we also consider some possible corrections to the 
NSR theory using the RG flow equation. We will invent the better approximations motivated 
by those analyses. 

Again we use the unit 2m = 1 for simplicity of calculations. 


4.5.1 Deep BCS regime 


Let us consider the asymptotic behavior in the limit (kpag)! —> —oo. We call this the deep 
BCS regime. At least in the vacuum, fermions interact via the weak attraction. At first, assume 
that an effective cutoff scale k? is much larger than the temperature T and the chemical potential 
u. Then the flow equation does not change so much from that in the vacuum, and as long as 
k? >> T, u we may approximate the effective action T% as 


(T) (T) 
E= f BO OW +o f orv ari PrP intette (4.5.1) 
p pp" ,q' 4 
1 1 k 
— TS(P=0) Sras 6r? 
should break down if the cutoff scale becomes of the order of physical scales T, u. 

At first, let us evaluate an amount of the self-energy correction. It would be natural to 
input the external momentum of the order of y/u ~ kp, then for evaluation of the self-energy 
we should use the effective coupling at some scale k ~ A ~ kp since the external momentum 
behaves as an IR cutoff for the coupling in the self-energy diagram effectively. At the 1-loop 


level, we find that 
eee Jkr e—il20 
Q. il? +P — u 


7 a 
Te J Ons FC $1 


_ 4,/pias/7 = : x? 
> i | UG eee) 


Here G51(p) = ip? + p? — u and gz" Clearly, this approximation 


Therefore, it roughly gives the correction of the Fermi momentum kp from a naive estimation 
u = kz, since kê = u+ X in our approximations here. 


= k? Ea aia LES ipa (4.5.3) 
MEE an Fos} ~ Om, 3n OS Di 


Since the chemical potential can be regarded as an energy which is necessary to add a particle 
to the system, this indicates that due to the attractive interaction u is smaller than the Fermi 
energy £p = k? /2m. Here we have derived the leading term in terms of kpas and the momentum 
integration of l is evaluated at T = 0. Anyway, we have a that in the small coupling region 
the self-energy just shift the chemical potential u in Eq. (4 to the Fermi energy €p = k?./2m. 

Since we have found that the Fermi surface is a aa notion in the deep BCS regime 
(kpas)~' < —1, we can use the Landau Fermi liquid theory described in sec [3.3] A consequence 
of those analyses is that if we have integrated out the modes except for those in the vicinity 
of the Fermi surface then the scaling ansatz works well. In this scaling regime, only the BCS 
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diagram can contribute to the superfluid phase transition and other diagrams, ZS and ZS’, drop 
out. 

As a first step, let us totally neglect the effect of the ZS and ZS’ diagram since those diagrams 
cannot contribute in the both high- and low-energy regions (Fig 4.7). In this approximation, 
we can evaluate the flow of the coupling. Again only the spin-singlet part of the four-point 
vertex survives in this approximation, and the relative momentum dependence of it does not 
exist (see SecA. 2). Then, the flow equation becomes 


Pe Lm ~ lt) rsp) 
OLTE (P) ~ —ók a = -a, | (G1() — Ey (l) + Ra)? (4.5.4) 


Since in the deep BCS regime we have already found that the self-energy correction is quite 
small, let us neglect X% in the flow of the coupling. In this approximation, we get 


1 bo f| 1 E 1 | (4.5.5) 
Ti(P) gas Ji L(G + Re)? (G30) + Ra) o 
We find that in the limit A, > œ and at k = 0 
2 2 
1 1 +f dêl 1-nr ((2-1) -u) -nr (E41 - 1) LI (466) 
T,-0(P) 8ras (27)? iP? — + (2 - i P(E D 21? B 
1 
The Thouless criterion gives T5(P =0) = 0, which leads us to the condition 
1 2 f% tanh ((e—u)) 1 
-— =- d = —~—|. 4.5.7 
ag a vade | 2(e — u) 2e ( ) 


This is the gap equation at T = T, used in the NSR theory, which just reproduces the result of 
the BCS theory as we have seen in sec[8.3| Indeed, in the low temperature limit and putting 
u = kz, we find that the critical temperature T, is given by 


T 8eF : T 
c= ex 
mez? i 2kras 


) ~ 0.613... X Ep exp (- a ) (4.5.8) 


2krlas| 


where ep = k? /2m is the Fermi energy and yp = 0.5772... is the Euler constant. 

From now on, we have totally neglected the contributions of the ZS and ZS’ diagrams 
although we have no reason to neglect them in an intermediate energy region characterized by 
k? ~ T, u. Let us evaluate the effect of those diagrams 


and then we will reproduce the leading correction to the BCS theory originally suggested by 
Gorkov and Melik-Barkhudarov [73]. Similarly, the self-energy correction also emerges from 
the same energy scale as we have seen in Eq. (4.5.2), but we neglect its effect to the RG flow of 
the coupling constant for simplicity. Recall that low energy excitations occur around the Fermi 
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Energy scale 
A 
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RG flow is described by the flow in vacuum. 
Only the BCS diagram renormalizes effective 
couplings between fermions. 


| Matter starts to affect the RG flow. 
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Figure 4.7: Separation of energy scales in the deep BCS regime. 


surface. Therefore, a typical scale of the center of mass momentum incoming to the effective 
coupling in the ZS or ZS’ diagrams is again of the order of \/H ~ kp. Furthermore, the transfer 
momenta in the ZS or ZS’ diagram is also of the order of kp. Let us assume that we can use the 
1-loop approximation for these diagrams in the weak coupling region, then their contribution 
el to the effective coupling becomes 


ZS+ZS' T io (4.5.9) 
fi 1 GOG + D 


where |g°| ~ 27T and |q| ~ kr. Here we have not restricted the possible configurations of the 
loop momenta l, but this expression automatically satisfies desirable properties. To see this, 
let us perform the Matsubara sum, which gives 


g@St28! = g? f d’l npl) -nr((l+q)’) 
ir ke J (27) ig’ +(L+qP?—-P ` 


(4.5.10) 


In the high-energy region specified by 1°, (L+ q)? > k}, the Fermi distribution function np 
suppresses the integrand exponentially and this is consistent with the fact that the ZS and 
ZS’ diagrams vanish in the vacuum. On the other hand, in the low-energy region specified by 
(L+ q)? ~ ki, the possible configurations of the loop momentum I is too restrictive and 
the quantum correction from this region should be small. This is consistent with the analysis 
in sec [3.3] which implies that the ZS and ZS’ diagrams does not contribute to the superfluid 
instability in the scaling regime. Therefore, we can crudely regard that Eq. (4.5.10) represents 
the loop correction coming from the intermediate region and that it roughly gives the correction 
to the effective coupling g;,,, used in the BCS diagram in the scaling region. 

Now, let us evaluate the effect of the ZS and ZS’ diagrams more quantitatively. In Eq. (4.5.10), 
the result depends on the length of transfer momenta |q|. Since we plug the result into the BCS 
diagram evaluated in the lower energy region, it would be natural that we apply the projection 
to the s-wave scattering. Since the transfer momentum is represented as q = p+ p’ with 


incoming and outgoing momenta p, p’, we redefine the effective contribution g as 
gS TZS = g? [ d cos 4 f dl np(l?)-— ne((l + 4(9))*) (4.5.11) 
i wda 2 J aF ig? lte a 
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where |q(@)| = kr./2(1 — cos 0) and @ is the relative angle between p and p’. For simplicity, 
let us evaluate the loop integration in Eq.(4.5.11) at T = 0, then it becomes 


Í dl np((l-— q/2}°) — nr((l+q/2)}°) 

(27)?  (L+q/2}? - (l - q/2)? 

_ T Pdl [ des 6(k2, — (P — lq cos 01 + q?/4)) — 0(k2. — (I + lq cos 0i + q?/4)) 
o 4r? Ja i 


2lq cos 0; 
1  dcos6) q? 
S A yi- qu — cos? 6). 
We now find that Eq.(4.5.11)) becomes 
' k  dcos@ o 1 
ZS+ZS" __ „2 F l 1- 5 a {= 20 
Ikr Ikr 78 D 5 L 5 a cos 0)(1 — cos? 0;) 
14+1n4 
= N(Er) 9p 3 


Here we have introduced the density of state N (ep) = kr /4r? = mkp/2n?. Therefore, at the 
lower energy scales we should use the effective coupling constant replaced as 


1+ln4 1+1n4 2kpras 
— = kp 1 
3 3 T 


 Ikp X 0.795... x (kras). (4.5.12) 


) . (4.5.13) 


This replacement affects the exponent in Eq.(4.5.8) in the leading order to add a constant 
1+1n4 


Ikr ™ Ikr = Gas = kp (1 + N (EF)Jkp 


. Therefore, the critical temperature is reduced by a factor (4e)~!/3 ~ 0.45 and then 


we should replace Eq. (4.5.8) by 


T 
T, = Fie 2 (tne exp (a) (4.5.14) 
in the deep BCS regime (kras)! < —1. Using the idea of the RG flow, we have found that 
the Gorkov and Melik-Barkhudarov correction is reproduced in a quite natural way. What is 
important is that, since we have a scale parameter k, we can flexibly change approximation 
methods at each scale and find physical insights of important diagrams. We expect that these 
insights obtained here will help us in going beyond Gorkov and Melik-Barkhudarov corrections, 
especially when we approach to resonance regions from the BCS side. 


4.5.2 Deep BEC regime 


Next let us consider the opposite limit (krag)! >> 1, which is called the deep BEC regime 
In the vacuum a bosonic dimer is formed as a spin-singlet bound state in this region and its 
binding energy Æp is given by 1/maz according to the analysis in sec|4.2| which is large in the 
deep BEC regime. 

In the vacuum, the effective action is given by Eq. (4.2.15) in which the chemical potential u 
is put —1/a%. In order to clarify the matter effect, let us introduce another chemical potential 
a defined by u = —1/az + pa/2. 


4We use the same notation kp as a length scale given by the number density n via the formula n = k}./37? 
even not in the BCS regime. We also use the notation £p as a unit of energy scale: ep = k?,/2m = (30?n)?/3 /2m. 
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Again, we evaluate the self energy correction with the 1-loop approximation with the full four 
point vertex in vacuum. We should emphasize that this is still a nonperturbative approximation 
motivated from the RG analysis, and it is consistent with a Dyson-Schwinger equation. 


—il?0 


l, p 
(7) : ; 
Dp) = ( ) = Trg- l 
(p) a ; , a! ae ea tolp +1) 


a2 . m 
__ lr fg st + yh + F + ip? + (p +1)?/2 — Ha) (4.5.15) 
as |) MER- pO Tp + AP 


Let us perform the Matsubara sum, which gives 


1+ 1+ 2G + ip? + (p+ 02/2 — pa) 


1 y e7110 
b ogar BOTË + 1 /ag — pa/2]fil + ip? + (p + 1)?/2 — pa] 


a2 . 
T 1+ $u + ip + (p +0)?/2— pa) meee 
Ja miee +1 [-w +P + 1/a% — pa/2][-w + ip? + (p+ D)?/2 — pd] E 
Here C, is a contour enclosing the imaginary axis in the w plane. See Fig[4.8] and we can 
change the contour of the integration as shown there. Since we assume that the binding energy 
is large, 1/a2 > ua/2 and then the branch cut in Fig/4.8| does not intersect the imaginary axis. 
Therefore, we have found that 


1+ 1+ ŽO + ip? + (p+ 0)2/2 — pa) 


1 5 e7 ilo 
B oaar DOP + 1/as — wa/2|[el? + ip? + (p + 1)?/2 — pa] 
2 


=e (P+ 2 — Ha) oT ee a P 1a — pal? 


yu [1+ & (ip + (p + 02/2 — B — 1/02, — pa/2) 
ip? + (p +1)?/2 — l? —1/as — pa/2 


(4.5.17) 


—np(P? + 1/02 — ua/2 


a2. » 
E dw =} 14 yı H (~w + ip? + (p + 1)?/2 — pa) 
ca 2Ti ef” + 1 [—w +l? + 1/03 — pa/2][—w + ip? + (p + 1)?/2 — ua) 


where Cù is the contour wrapping around the branch cut (see Fig/4.8). In calculating this 
result, we should notice that p° is also a fermionic Matsubara frequency. The second and third 
terms in the right hand of Eq. (4.5.17) are exponentially suppressed in the limit 3/a23 —> oo. 
Let us neglect these higher order corrections and simply approximate the self energy as 


1)2 
(p) = =| a 2np (E3 — pa) (4.5.18) 
PT as J CTP ip + (p+ 02/2 -P — 1/03 = pa/2 = 


We can now evaluate the number density n using the formula of Eq.(4.3.11) with the self 
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w = —il° 


Re > 


Figure 4.8: Change of contours in the Matsubara sum. It also reveals the pole positions and 
the branch cut of the integrand. 


energy given by Eq.(4.5.18): 


o (Pp) — E(p) 
(T) —2e7ip0 

— 4.5.19 

l i PP oe = l ) 

l 2 
ft ah ran Coroa) 
z 3 38 J 2 : 
(27)? (27)? B m fip + p24 2 = a fip + oy _ pe x — a 


We expand the full propagator in terms of the self-energy X(p) in this expression, and we 
will discuss what happens when we do not expand it at last of this section. The first term 
in Eq. (4.5.19) gives the Fermi distribution function. Since fermions are gapped by a half of 
the binding energy, this term is exponentially small. Let us evaluate the second term. The 
Matsubara sum gives 


2 
2 = 2? 
serer g [weg] (eo 
(4.5.20) 


where we again neglect exponentially small contributions in terms of 3/a%. We obtain that 


2 
/ Bp dl (320 /ag)np (e2 =f. = _ us) np (22 O pa) 
n= 


(27)3 (27)3 [(p — 1)2/2 + 2/a2]? 
ËP (P? da Btne (Z -(F-aP-2-#) 
7 >| Or P ( 2 pa) J Ga [(2q)2/2 + 2/a2]? oe 


Here we have changed integration variables so that P = p +l and q = (p — l)/2. Let us 
evaluate the integration over the relative momentum q. Since the center of mass momentum P 
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should be of the order of VT due to the Bose distribution function, we may put it to be zero in 
the Fermi distribution function. As a result, since we may approximate the Fermi distribution 
function as ne(—q? — 1/a3) œ 1, the relative momentum integration gives 


1 d°'q  327/as _ 
(2m) (2g? + 2/2)’ 


Therefore, number density becomes 


n=2 I Sm (= = pa) ; (4.5.22) 


We have confirmed that the free bosonic picture naturally emerges in the deep BEC regime. 
In the deep BEC regime, the biding energy and the energy scale of thermal excitations are 
well separated. Therefore, the change of the RG flow from the vacuum one can occur only in 
low energy scales. As a first approximation, let us neglect its effect to the four point vertex. 
Then only the BCS diagram contributes and the Thouless criterion is again given by Eq. (4.5.7). 


Let us substitute u = —1/a2 + Ha/2 into Eq.(4.5.7) to get 
S 


1 Bye tanh (&(e +1/a& — ua/2)) 1 
E. vaa | Ae +- a2) a m 


Since 3/a%, >> 1 in the deep BEC regime, we may put the hyperbolic tangent in the integrand 
to be 1 neglecting exponentially small corrections. In this approximation, the right hand side 
becomes —,/1/az — fig/2 and then the Thouless criterion gives pa = 0 at T = T,. 

Now we can calculate the transition temperature T, substituting ua = 0 into Eq. (4.5.22). 
Since n = k3./37?, we get 


kr qd? P 1 1 jr 
a = TH?) 56(3/2). 4.5.24 
3r? eee c = 560012) (4.5.24) 


Since ep = k% in our unit 2m = 1, we find that 


9 1/3 
T./ér = (sa = 0.218... (4.5.25) 
This is the BEC transition temperature of the free bosonic system. 

Here again, ideas of scale separation were so powerful. Thanks to the fact that RG flows in 
the vacuum and in matters coincide at high energy scales k >> n!/3, we can concentrate on low- 
energy dynamics even if we use purely fermionic formulations in the BEC regime. This enables 
us to extract nonperturbative features in our calculations, and to discuss phase transitions of 
BEC. In order to go beyond this approximations, it is very important to take into account the 
interaction between dimers and for this purpose calculating dimer-dimer scattering in terms of 
the flow equation is necessary. 


4.6 Summary 


In this chapter, we consider the two-component fermionic system with attractive contact inter- 
actions. Our central goal of this chapter is to calculate critical temperatures of superfluid phase 
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transitions from the knowledge of scattering physics without putting artificial assumptions in 
those calculations. The main idea of those analysis is the following: in the high energy region 
important physical process should be approximately same with that in the vacuum, and effects 
of many body physics mainly affects low energy physics only. In order to realize this idea it is 
quite natural to separate energy scales, and by employing FRG we can consider the BCS-BEC 
crossover in a systematic way. 


At first, in order to find properties of physics in the vacuum, we calculate scattering process 
using FRG. Due to the absence of antiparticles or holes, the FRG equations automatically 
close at each sector and are free from infinite hierarchy. Thus, in principle we can calculate 
the scattering process in a precise way. As a result we have seen that the two body sector in 
the vacuum is precisely solved and we have obtained the four-point vertex in the vacuum. We 
also discussed the atom-dimer scattering by studying the three-body sector, and derived the 
flow equation for a new vertex whose dimer legs are amputated. Important lesson learned from 
atom-dimer scattering is that, feedback from higher-order vertices exists in general cases but 
still flow equations decouple from higher-body sectors. This is very important property which 
is missed in previous studies. 

We derived some formula to calculate number densities n using FRG, and also derived 
the Thouless criterion using the Ward-Takahashi identity. Especially, we discussed the Ward- 
Takahashi identity in a general way and showed that it can be understood as an integrated flow 
of the corresponding correlation function as long as phase transitions do not occur in the RG 
flow. We also proposed that for symmetry broken theories there exits a critical cutoff scale k, 
at which the second order phase transition occurs by crossing the critical surface, although it 
is not still a mathematically rigorous statement. 


Finally, we reproduce the NSR theory by considering the deviation of the RG flow of many 
body fermionic systems from that of the vacuum. In the BCS regime, high energy physics 
is controlled by scatterings in vacuum and low energy physics is controlled by the Landau 
Fermi liquid theory. Therefore, it is easy to control resummation of Feynman diagrams based 
on scaling ansatz on those energy scales, and by considering the higher order effect in the 
intermediate regime we have reproduced the Gorkov and Melik-Barkhudarov correction in a 
natural way. 

In the BEC regime elementary excitations are bosonic dimers, and in order to describe the 
physics based on fermionic theories, nonperturbative description is required even at the lowest 
order approximation. Based on knowledge of the RG flow in the vacuum and separation of 
energy scales of many-body physics from that in the vacuum (binding energy), we evaluated 
the fermion self energy correction so as to find a bosonic picture. As as result, we can obtain 
free bosonic picture in our formalism, and then we can consider scattering between dimers as 
a future work. 


As a concluding remark, let us comment on previous studies by other groups. Researches on 
the BCS-BEC crossover using FRG is done in the paper by introducing the Hubbard- 
Stratonovich field in the particle-particle channel. In this case, one can apply the local potential 
approximation to the effective potential of the bosonic field and they gave quantitatively correct 
results in the BCS and BEC regions. In order to circumvent ambiguities due to the Hubbard- 
Stratonovich transformations associated with the Fierz transformation to some extent, they 
invented the rebosonization technique to eliminate dynamically generated fermion quartic cou- 
plings [76]. Their numerical calculation gives almost correct dimer-dimer scattering length, first 
calculated by D. S. Petrov et al. E and asymptotic formula of critical temperatures of 
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interacting Bose gas so called Lee-Yang-Huang correction [481]. 

Without introducing the Hubbard-Stratonovich field, R. Haussmann calculated transition 
temperatures of the BCS-BEC crossover [82||83}. Later, his group established the variational 
method with the 2PI formalism for the BCS-BEC crossover, which is called the self-consistent 
T-matrix method [84]. Since their approximation breaks the global U(1) symmetry the phase 
transition becomes first order, but their analysis provides plausible numerical results in the 
BCS and BEC regions. Especially, by introducing the external fields for two point functions 
the description of the BEC region becomes transparent. These are not in the framework of 
FRG, but combining their analysis to our FRG method would be an interesting problem. 
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Chapter 5 


Application of FRG to dipolar 
fermionic systems 


So far, most ultracold atomic experiments are performed with alkaline atoms, and in such 
systems the dominant interaction can be treated as a contact interaction (see sec/4.1). In 
this chapter, we consider a totally different system: ultracold atoms with strong dipole-dipole 
interactions. In such systems the interaction is long-ranged and anisotropic, and we can expect 
unprecedented phenomena absent in conventional systems with the contact interaction. 

Experimental successes in realizations of ultracold dipolar gases reported in the papers 
encourage us to perform theoretical calculations in detail for such systems. These are 
very recent progresses in this area. Indeed, those papers were published in May 2012. Therefore 
we will explain them in sec. 1] and especially we will see details about the fermionic dipolar 
system realized with dysprosium-161 [86]. 

Our original motivation to study dipolar systems is to reveal physical properties of dense 
nuclear matter. Although such dense systems of nucleons have not yet been realized in experi- 
ments, it is believed to be realized inside of neutron stars. In the dense system, the strong and 
anisotropic nature of the nuclear force become important and it is theoretically predicted that 
many unobserved phenomena such as nucleon superfluid, meson condensation, etc. are real- 
ized [87]. Acquiring new insights for dense nuclear matter from ultracold atomic experiments 
is our purpose of this study. 

In sec{5.2| we will see properties of the dipole-dipole interaction briefly and it will provide 
physical insights for our system. In sec[5.3] we establish the Landau Fermi liquid theory with the 
dipole-dipole interaction using RG methods and find a new kind of superfluidity characteristic 
to a dipolar fermionic system. 


5.1 Ultracold atomic systems with the DDI 


Let us explain experimental current status of dipolar gases in this section. The dipole-dipole 
interaction (DDI) between two magnets separated by r is given by 


2 
a Y Q Q a Zs 
Vaal") = 75 Ee - $55, — 3(S1%, - #)(S28, -#)] , (5.1.1) 
where S; represents the spin vector of a particle i. Here y is the coupling constant for the 
DDI, called the magnetic moment or the gyromagnetic ratio. It is often measured in the unit 


of the Bohr magneton upg = A Since the spin and orbital angular momenta are coupled in 
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Figure 5.1: The dipole-dipole interaction (DDI) between two magnets. 


Eq. (5.1.1) it reveals an anisotropic nature (see Figf.1). We should also notice that it decreases 
as 1/r?® as the distance r increases, and its behavior is much slower than that of the van der 
Waals interaction which decreases as 1/r®. 

The first experimental study of quantum dipolar gases was the realization of the dipolar 
BEC of chromium-52 (ł4Cr) by the Stuttgart group |88|. Chromium is a transition metal 
with a magnetic moment of 6ug. In order to make the DDI between the chromium atoms 
dominant, the researchers used the Feshbach resonance to turn off the contact interaction and 
they observed the d-wave collapse of dipolar condensates [59]. 

In the field of the dipolar BEC, the BEC of erbium-168 (!$$Er) is reported in May 2012 
by the Innsbruck group [85]. Erbium is a lanthanide atom with the magnetic moment Tug. A 
remarkable property of erbium-168 is that it has 6 Feshbach resonances in low magnetic field 
environments. Indeed, they observed the d-wave collapse of the dipolar BEC in erbium about 
1 Gauss, although very high magnetic field lager than 500 Gauss is necessary for the case of 
chromium-52. This may open possibilities to study a purely dipolar bosonic gas thanks to the 
weakness of external fields. 

Next, let us care about fermionic dipolar gases. In May 2012, the achievement of the 
quantum degenerate dipolar fermionic gas was first reported with dysprosium-161 ('@sDy) by 
the Illinois and Stanford group |86|. This had been one of the unrealized experimental challenges 
in ultracold atomic fora group FA Therefore this step is qualitatively new, and the observation 
of the Fermi sea of the dipolar fermions itself is a very important result (see Figf.2). 

Now we will explain the experiment of dysprosium-161. Dysprosium is a lanthanide atom 
and its isotopes '®!Dy and '®*Dy are the most magnetic fermionic atoms with the dipole moment 
104g. Usually, achieving the quantum degeneracy of identical fermions is difficult, because the 
Fermi statistics prohibit s-wave scattering with the short-range interaction and because the 
rethermalization process from the p-wave collision ceases due to the centrifugal barrier below 
some threshold temperature about 10-100uK. 

From this fact, they chose the sympathetically cooling of '*'Dy with its bosonic isotope 
162Dy, Both are spin polarized into high-magnetic-field seeking states: |F = 2, Mp = —2) 
for 1©'Dy and |J = 8,m y = —8) for 1*Dy, where F = I+ J is the total spin, J the nuclear 
spin PI and J = 8 the total electronic angular momentum. By switching off the optical dipole 
trap gradually, spin-polarized !°'Dy are evaporated to a quantum degenerate state. Fig [5.2] 
shows the result, and the population shows the better coincidence to that of the Thomas-Fermi 
distribution, which reflects the Fermi-Dirac distribution, than that of the Gaussian fit, which 
reflects the Boltzmann distribution. The estimated temperature is T/T = 0.21(5). 


lFor other highly magnetic fermionic atoms °°Cr (6g) and !°’Er (7g), the temperature is still not below 
10uK. For the fermionic polar molecule *°K®’Rb (0.57D) the temperature is about the Fermi temperature 
(T/Tr ~ 1.4), but the dissociation process from chemical reactions prevents further evaporative cooling. 

#7 = 0 for "Dy, and T = 5/2 for "Dy 
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Figure 5.2: Quantum degeneracy of the dipolar fermions [86]. (a) Single-shot image of the 
population in the momentum space. (b) Average of six images. (c),(d) Density integration along 
the horizontal, vertical axis, respectively. Green curves represent the fit with the Gaussian, 
while the red curves are fit with the Thomas-Fermi distribution. This reveals that T/Tr ~ 0.2. 


We should also emphasize that they have also achieved the quantum degeneracy with the 
forced evaporative cooling of spin-polarized '©'Dy without the bosonic isotope '°?Dy. Their 
analysis shows that T/Tr = 0.7 with Tr = 500nK. This would suggest the importance of 
dipolar collisions and they might provide rethermalization processes. 

In following sections, we will consider two-component fermionic system with the DDI. Hith- 
erto, such systems have not been realized but the recent progress is opening the new physics 
with dipolar systems. We can expect the realization of similar systems with the model which 
we will consider. 


5.2 Basics of the DDI 


At first, let us remark on basic properties of the DDI. We consider the magnetic dipole inter- 
action (5.1.1) between spin-4 fermions. We should be attentive to following properties of the 
DDI: 


e The interaction is invariant under the combined SU(2) spin rotation and SO(3) spatial 
rotation, however the orbital angular momentum L and the spin S are not conserved 
independently. 

e The DDI has a long-range feature, i.e. it damps as 1/r?. 


Since we will frequently use the momentum-space representation, it is convenient to give the 


Fourier transform of Eq. (5.1.1): 
Q 1 Q ^ ^ Q 
Vaala) = <7 [3(S% ASE â) — Sè- S$ |. (5.2.1) 
The interaction term V in the Hamiltonian is given in the coordinate space as 


gg 5 : PrdyYale) paly) Vzw lE — yhe lya (z). (5.2.2) 
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In the momentum space, Eq.((5.2.2) becomes 
1 — — a 
= J f Tapii awe aa Ak — k'y Pj- Ya, P/24k': (5.2.3) 
P Jk,k' 


We decompose Eq. (5.2.3) into some channels with spherical harmonics. Here we regard that 
the total momentum P = 0 and that all the relative momenta k and k’ lie on the Fermi sphere. 
Therefore, the tensor interaction Vf, (k — k’) is a function only of directions k and k’. We 
can find details of calculations in Appendix [E] to get an explicit form of the interaction matrix, 
but here let us discuss properties of the DDI to get a physical insight of the result of those 
calculations. 

For later use, let us consider decomposition of more general interactions Ve AR: k’) and 
obtain expansion formulae. We define 


(S3') S's! a 
Veo (kk) = >- > Cie gCa a ipi [Vein (k — k’) — Vga (k + k')]. (5.2.4) 
a,B a yc 
The inverse formula of Eq.(5.2.4) is given by 
Vik ke) = =S DO ae P a i mies (k, k’). (5.2.5) 
SS! S.S! 


Since we assume here that the spin of particles is 1/2, the Clebsch-Gordan coefficients only takes 
the form C?°:,,, and S and S’ only take the values 0,1. In the case of higher spin species, we 


2 a," 
can obtain a similar formula. Decomposition with respect to relative angular momenta gives 
(SS’) A A 
Ves '(k, k= Dy » Vins, rms YLM(R)Y py (k'), (5.2.6) 
LM L'M' 


where coefficients are given by 


ge =j dO; f POW ulk) Yum (BV E> (k, k^). (5.2.7) 


We can interpret this quantity as a tree-level transition amplitude from the channel 
labeled by (L', S"; M’, s’) to another channel labeled by (L, S; M, s). Since this is the decom- 
position of the interaction for general cases, the above statements is also applicable to effective 
interactions after a cutoff scale is lowered. 

Let us now restrict possible terms in the effective actions using the symmetry arguments. 


e The Fermi statistics, Y, ade, y= oe. vUa p implies that under the interchange of the 
spin variables a + s — a in Eq.(5.2.5 5.2.5) associated with k + —k should give a factor (—1). 
This suggests that the oe 8", urs does not vanish only when L is odd and 
S = 1, or when L is even and S = 0. This should be true for L’ and S’, too. 


e Next, let us consider the invariance under the parity transformation Y(w, k) > w(w, =k) 
and similar for 7. This implies that the interaction in Eq. (5.2.6) is invariant under the 
transformation k +» —k© and then oa wy does not vanish only when (—)’+"" = 
+1. This means that L+ L’ should be an even integer, i.e. L = L,L+2,L+4,---. The 
combining the above two arguments, S is a good quantum number for our model. 
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e The rotational symmetry implies conservation of the total angular momentum J? = 
(L + S)? and its projection quantum number J, = M + s. Since the total spin S is 
also a good quantum number, we can consider the spin-triplet couplings and the spin- 
singlet couplings independently. For S = 1, we need to require J = J’, which implies that 
L = L',L'+1, L’+2, and then the invariance under the parity suggests that L = L’, L' 2. 
For the spin singlet case, we can conclude that L = L’ and M = M”. 


Since L and S are not good quantum numbers in general, and furthermore M and s are 
not good quantum numbers even in our model, we would like to use the representations which 
diagonalizes the total angular momenta J and its projection J,. For this purpose, let us define 


the coefficients 
vised Jt JJ J, S’) 
Vis. L's! =) X ) CiM,ss.C ae S/S! ae L'M's!s (5.2.8) 
M,Sz M',S! 


or more explicitly we have 


JJa I, — JJz Ors: SS: Sst 
VESS = ` ) ` ) CEM.ss.C L'M',S' S! X D C AO sf! 


M,Sz M’,S! aß a'b! 


x f ao | dOQwY p(B) ow (k' [V74 (Ke —k')—-V. Dk +k’)]. (5.2.9) 
This is explicitly calculated in Appendix [E] 
The DDI in Eq.(5.2.2) is non-local, however we can obtain the local interaction by going 


back to the original action with the magnetic field. This procedure is equivalent to eliminate 
the quartic coupling (5.2.2) with the Hubbard-Stratonovich transformation: 


= f a fee (T00 Fue) +} vas") 


+f ar | ae y (Pv) -(V x A) (5.2.10) 


To see this, let us integrate out the gauge field A from the action, then the effective interaction 
becomes 


exp (-[ drV(r 25 [PAcw|- [area (5074, + re nv) . (5.2.11) 


To see that the interaction V introduced here coincides with the one in Eq. (5.2.3), we complete 
the square with respect to A: 


[ee (5074)? rein S00) As ) 

= J (S Ai-aAig + Yei (PW) (iA) 

7 f T = (miy). 3) (Ano yews (Fv) =) 
E fener Qt (B84) (08%), 
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Therefore, up to an additive constant, the interaction becomes 
g? AA 0i — Oy 
V= fai- su) (00) PE). (5.2.12) 
q —q q 


Let us evaluate the Fourier transform of the spin composite operator: 


TOi 3, —i(k'—k)-æ —ip-x ip -aT 
Ee) =a Te ( ) J P [ie Ypi Yp 
= f Yp4kTiYptk' (5.2.13) 
p 
Now we have found that the traceless part in Eq.(5.2.12) becomes 
2 
Y T A A dij \ = 
Vtraceless = rf , Wor quilp (4a _ s) Wy —q?5 Vp! 
P.p'.q 
1 = = 
= 2 , Va pta VB p -qVob b'w (OQ Up p Va! ps (5.2.14) 
P.P'.q 


which coincides with the one in Eq.(5.2.3) after the change of the variables P = p+p’, 
k’ = (p — p’)/2 and k = q+ k’. On the other hand, the trace part in Eq.(5.2.12) gives a 


contact term: 
Vine = -T S (Eg) 00), 
= — f Pa (Fv) (æ) (Gv) (æ). (5.2.15) 


Since there usually exists a centrifugal barrier, and such a contact term does not play an 
important role. The classical equation of motion implies that 


1q; (—O 
Aig = Yei A (=y) ; (5.2.16) 
q 2 q 


or 


-V° A = 7V x (Y Sy). (5.2.17) 


We should also notice that V- A = 0 from Eq. (5.2.16), and this is a consequence of the fact that 
the action in Eq.(5.2.10) describes the dynamics of the magnetic interaction with the Coulomb 
gauge condition. 


5.3 RG study within the RPA 


We at first study the weak coupling regime using the random phase approximation (RPA). Not 
only just applying the RPA, we also consider the justification of those analyses based on the 
RG flow. The following discussion will be performed in a parallel way with the analysis in 
sec[3.0] 

Basically we rely on the Landau Fermi liquid theory again. We assume the existence of the 
Fermi sphere and suppose that low energy excitations live in the vicinity of the Fermi surface. 
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Due to weakness of the interaction, we neglect renormalization of the coupling in vacuum. 
We also neglect the self-energy correction, since unless the parity invariance is broken the self 
energy correction just shifts the chemical potential u to the Fermi energy €p. 

Then, in order to discuss the stability of the Fermi liquid theory, we need to check the 
absence of singularities in the four-point vertex. In the RG flow equation, there are three 
diagrams contributing to the flow of dipolar couplings: 


4 2/ i 3 {o> 

78 = 78! = _ BCS = oO. 
1 2 1 2 

3? 


We will confirm that ZS and ZS’ diagrams lead to generalizations of the Pomeranchuk insta- 
bility to dipolar systems, but in the weak coupling region such instability is absent. From the 
analysis of particle-particle fluctuations coming from the BCS diagram, the superfluid instabil- 
ity is induced in the low temperature we will find that the new type of superfluidity (the ?P, 
superfluid) appears. 


5.3.1 Study of the Landau channel 


Let us discuss whether the diagrams ZS or ZS’ can affect low energy physics. We discuss the 
ZS diagram at first: 


In order for the diagram not not to be suppressed in terms of the cutoff scale, it is required 
that both of vertices denoted by gray squares be automatically marginal at any loop momenta. 
In other words, it is necessary and sufficient for the ZS diagram to become marginal that the 
momenta kz (or equally k3) is not restricted after we specify the momenta kı and k}. This is 
only possible if kı ~ k{, indeed if so it implies that kə ~ ki and we may freely choose a point 
of the Fermi surface for ky because ki, automatically lie on the Fermi surface too. Similarly, 
the ZS’ diagram is not suppressed only when ky ~ ki and especially when ky ~ ky œ k4 ~ ki 
both of the ZS and ZS’ diagram contribute to the quantum effect, which respects the Pauli 
exclusion principle. 

We call the channel such that kı ~ ki as the Landau channel and let us analyze this 
channel. If kı ~ ka further, not only the ZS graph but also the ZS’ graph contributes in 
the same order of magnitude, however here we do not take into account the ZS’ diagram for 
simplicity of calculations. 

Let us perform the Matsubara sum in the ZS diagram. We neglect the dependence of the 
vertices on the Matsubara frequencies of the external legs, then we should calculate 


1 1 


1 1 
= G(ks)G(ks — = - - ? 5.3.1 
pe i) GM @) Ao ae) ee ( ) 
Here Q = (iQ, Q) represents the transfer momentum. The summation over the fermionic 
Matsubara frequencies can be written as 
1 dp? —1 1 1 
— G(ks)G(ks — Q) = -—— - a: (5:3:2 
z> SOUPER (ar E a aal 
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where C, is the closed curve enclosing the imaginary axis in Fig [5.3] Performing this contour 
integration, we immediately obtain that 


1 yy _ sinh $ (elks) — elks — Q)) —1/2 
B Ds SENG a= iQ + [e(ks) — e(ks — q)] cosh 5 (e(k5) — u) cosh 5 (e(ks — Q) — u) 
(5.3.3) 
In the limit |Q| < |ks|, we can approximate it as 
17 sinh (Sek ks - Q) 1 
B 5 Le kis)G Q) = 2 iQ+vupks-Q cosh? (ls) eae) 


and the naive limit Q — 0 is ambiguous as usual. This singularity reflects to singularities of 
effective couplings and two different low-energy effective coupling constants are defined: 


ap 1 a yl 

DRgwlk k’) = lim ara nk -Qk +Q] (5.3.5) 
Qal AA- J] aß I 

DA gra (k, k’) = lim [Tag(e ck Q+ Oaa . (5.3.6) 


By neglecting the ZS’ diagram, we found that the Q-limit V® is invariant under the RG flow, 
and that the Q-limit V@ obeys the flow equation 


A -A\ 12 
Qos ALAS j kpdls B 1 
NOT palk k ) = Ada (/ +f ) 2T? ( 4) cosh? £ 5UFls 


d?ks oa ' i 
J eR ray (k, Bes) ER y (ks, k’) 


CurA dks Q&I 
=N Te" (k, kpr? (ks, k' 3 
Crna) e iek kk) 687) 


where N (ep) is the density of states of a single spin component at the Fermi surface: 


1 dk k2 
N = — k — S 5.3.8 
(eF) 272 F deg e Onur ( ) 
Im 


iQ + e(ks =q) — 
X 


Figure 5.3: The position of the poles in the frequency plane p° = iws. 
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That is, d®k/(2m)° ~ N(ep)dé. By defining the dimensionless coupling constants 


F= an (5.3.9) 
we get 
AN or QY'B 1 
AA Vy k k’) = re AUN ak, ks) VX (ks, k’), (5.3.10) 


The channel decomposition in the Landau channel is given by 
a Qa . s S'S, f * f, 
Vzal khh =S Dros ul O Ver YE kVa (8.3.11) 
L,Lz S,Sz L' L, S'S} 
where A = (L, L}, S, S) and A’ = (L', L, S', S1). The angular integration in Eq.(5.3.10) 


becomes 


f dR" YO Ve ulk, k") VOLEC", k') 


yy 


QR s-a! WSS sty sus A : sy Q 
fa k ) X ) (— ) O al ) Cay: ;s— YLL (k)Yirru(k Vý an 
vy’ LLzSS:z L” LUS" sy 
= Se s+ S'S, A ¥ A , Q 
x ` ` (=) IO ee a ) Crane ;s— Yrm (k)Yi r, (k Vän ay 
UELS SUSY L'LL, SS, 


S Ss’ S} Aa 
a ae ee a coe. oe PE B; rp YLL (kK) You (k’) 
LLzSSz L' LY S'S, 


Si sy Sg 
ie s—y st+y 
x ) ` ` l A; L” LYS" SY Vit T wd. ie ) Cay; ;s— "a ) Cy; 33-7 


L” LY S's OOS 


s s—a ž s sist * o 
= (-)? we ) Ge. a! {= ) "O TT SVa a (5.3.12) 
A,A’ A” 


Let us introduce the matrix notation V? = Un 4), then it obeys the RG flow 


B 
UFA 2 
NAVE = (-)*#—5— (ve) , 5.3.13 
se (-) cosh? 2 sUFA i ( ) 
and its solution is given by 
V2 
Ve = g l (5.3.14) 


1+ (=) V£ (tanh Furl, — tanh Sur) 


Since A, is considered to be sufficiently large, we may approximate that tanh (GuvpA,/2) ~ 1, 
then we can readily identify the stability conditions for the Landau channels. Let the fermion 
spin s be 1/2. If and only if all the eigenvalues of V,! is less than 1, singularities in the 
renormalization group flow do not arise and the Fermi liquid picture is consistent. 
Especially when we restrict ourselves in the weak coupling region, the stability condition is 
manifestly satisfied and the Fermi liquid picture should be valid. 
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5.3.2 Study of the BCS channel 


Here let us study the BCS diagram to see the BCS instability, and let us concentrate on the 
RG flow in the BCS channel. For this purpose, put ky = —k, and ki, = —k’‘. In this case, 
Q =k, — kı or Q' = k, — kı need not be small and then the ZS and ZS’ graphs are suppressed 
due to the restriction of possible configurations of loop momenta. Therefore, we neglect these 
graphs and only take into consideration the effect of the BCS graph. In these approximations, 
the 1-loop contribution of the BCS graph to this channel is given by the flow equation 


No —A 2 tanh g 2p! 
Að TA, (k, k’) = Adj (-5] (| +f j= k odes an (5 UF 5) F k P (k; a (k", k') 
A —No 


2r? Qupes An 
N(er) B d2k" B 6 
= tanh | —upA ee k, k” ry a(k”, ke’ .3.1 
EP tann (Sued) S SE ree Rg (RH) (5.3.15) 
where N(ér) = ff is the mn EG BID o£ as for one spin projection. With respect to the 


dimensionless couplings given in Eq.(5.3.9), we get 


Ady Vee, (k, k’) oe 


1 p FIr " "yl 
= 5 tanh (2a [ve Ve" (ke, ke \V3% (k jk): (5.3.16) 


Let us calculate the angular integration in Eq.(5.3.16) using the decomposition (5.2.5]5.2.6). 
Substituting the decomposition of the interaction 


yo 1 Ka a’ +B'17(SS") a * an 
Ve k, k’) = >, > Cy aera ip ee a+; L'M' a! rg Yrm(k)Yi w (k ), (5.3.17) 
LSM L'S'M' 
we find that 
i ea Vie (RR Va (k") (5.3.18) 
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LSM L'S'M! i 
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LSM L'S'’M’ 
(S'S) (S”"'S’) 
x Vim a+; L!M" 3! Vi MI" s: Li M' ol 48" 
LYM" s", 
Therefore, the expansion of Eq. (5.3.16) gives the set of flow equations: 
SS’) 3 ) À V SS") S" S’ 
—NAOd,V, een L'M’'s! — T tanh Gay VO athe ” Pat E (5.3.19) 


L!M" S's A 
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Formally, we can solve the differential equations given in Eq.(5.3.19). Notice that it takes 
the form 


1 
AOAVA = 5 tanh (Sued) vVř, (5.3.20) 


where Va is the matrix with coefficients Poy m's at the scale A. The formal solution of this 
matrix differential equation is given by 


1 


va =v (1+v, i tanh (Sena) ) (5.3.21) 


Indeed, this solution satisfies at A = A, V = V, and Eq.(/5.3.20), as desired. The Landau Fermi 
liquid remains stable in the BCS channel at any temperature when the interaction in the BCS 
channel is marginally irrelevant, i.e. the bare couplings satisfy 


V, > 0. (5.3.22) 


On the other hand, if a condition in Eq.(5.3.22) is violated, the coupling becomes marginally 
relevant and a Landau pole appears at some low temperature. To evaluate the critical temper- 
ature, let us consider the low temperature limit. Then, the critical temperature is 
2e7E 1 
Tsc = ——vprA, exp | —— ], 5.3.23 

sc = re pP ( wn) ( ) 

where Vinax is the maximal negative eigenvalue of the matrix V. 
The interaction matrix elements is given by 


= et EPP ay (5.3.24) 
where 
Cie td 3 
Tisis = (=P sb. SVS + 1) (28+ 1) (5.3.25) 
gazos? 
xy (L + 1)(2E' + 1)(28 + 1)(28’ +1) l F s s 
LI J i] 


1 1 2 
x | (Hr + Hr) Cio, 100700 + (— 10D HiCtoao Chono d L L' l } i 


The details of calculations to obtain this result can be found in Appendix|E] In the case of spin 
s = 1/2, only the spin triplet channel S = S’ = 1 remains. We get 


: 3/4 0 7 -3/8 0 = 0.075 0.061 
J=0 _ J=1 _ J=2 — 
na ( 0 0 ) = ( 0 0 ) (ae ( 0.061 0.050 ) fe 


where we have picked up the values for L, L’ = 1,3 in (5.3.26). The maximal negative eigen- 


value of the matrix V is given by Vinax = —gN (er) and the corresponding channel is 3P}. 
T 
Therefore, we can expect the 3P, superfluid below the critical temperature 


2e7E 8 
T= A — . 
Eh e ( a) 
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Let us consider the situation in which we add the contact interaction 
7 g 
ve= $ | atout ouiet). (5.3.27) 


In this case, g > 0 means the repulsive contact interaction, and g < 0 the attractive contact 
interaction as usual. The channel decomposition shows that this remains non-zero only if J = 0 
and L = S = L’ = S' = 0. Therefore, the dipole-dipole interaction and the contact interaction 
do not mix at least in this approximation for the BCS channel, and the transition temperature 
of this channel (tSo) is given by 


IQeVE 


e 1 
vrAexp | -—— ), 5.3.28 
r F p( E Ce 


if g < 0. We can judge the relevant channel by comparing the relevant scale for those channels. 
The border of those two phases are given by 


T= 


3 


y? 
= -—., 5.3.29 
gl = 5a (5.3.29) 


5.4 Summary 


From the RPA analysis of a dipolar fermionic system with spin 1/2, we have found two stability 
condition. The first one corresponds to a generalization of the Pomerunchuk stability condition 
(3.3.20) of the usual Fermi liquid to the dipolar system. This instability is absent when the 
coupling is sufficiently weak. The second one is the superfluid instability, which appears if the 
temperature is sufficiently low even when the coupling is very weak. 

The phase diagram of our analysis is shown in Fig5.4 The vertical axis shows the strength 
of the dipole-dipole interaction and the horizontal axis shows the coupling of the contact inter- 
action. When the attractive contact interaction is dominant, the system shows 'So superfluid, 


Break down of Fermi 
Liquid picture 


yE 2/7 m N 
S > 
"o O 
superfluid | superfluid `~ 
o i 0 1 
N(er)g 


Figure 5.4: Phase diagram of a dipolar fermionic system in sufficiently low temperatures. The 
horizontal axis is strength of dimensionless contact couplings, and the vertical axis shows that 
of the dipole-dipole interaction. The green region shows emergence of !Sọ superfluid due to 
the attractive contact interaction, and the blue region shows that of °P, superfluid due to 
the dipole-dipole interaction. The white region shows break-down of the Landau Fermi liquid 
theory due to the Pomeranchuk instability. 
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and it is shown as the green region in the figure. On the other hand, if the dipole-dipole inter- 
action becomes dominant it shows °P, superfluid as shown with the blue region in the figure. 
This type of superfluidity has not been found in other systems, and it reveals characteristic 
behaviors of the dipole-dipole interaction. 

Before closing this chapter, let us comment on the relation of our studies to previous works 
of other groups. In Fig [5.4] the white region represents break-down of the Landau Fermi liquid 
picture due to the extended Pomeranchuk instabilities of Eq. (3.3.20). These instabilities are 
already shown in the paper by T. Sogo et al. using the RPA analysis at T = 0. They 
introduced bosonic auxiliary fields in particle-hole channels and identified those Pomeranchuk 
instabilities using mean field approximations. However, both of their analysis and ours can 
only be justified in the weak bare coupling region, we need more refined calculations to reveal 
the correct physical origin of those instabilities. 

In strong coupling regions, unbiased analyses have not been done and we can only access 
variational calculations at the zero temperature with some ansatz on ground state wave func- 
tions. In the paper by B. M. Fregoso and E. Fradkin, they revealed that the spin-polarized 
state is energetically stable if the contact repulsion is so strong, and that between the param- 
agnetic and ferromagnetic phases the intermediate state called ferronematic states is a possible 
ground state. In the paper by K. Maeda et al., they propose another ground state, called 
antiferrosmectic-C phase, motivated by the alternating-layer-spin (ALS) structure asso- 
ciated with meson condensations in high-density nuclear matters, and they showed that it is 
energetically stable if the dipole-dipole interaction and the contact repulsion are strong enough. 

It is important to confirm that whether those possible ground states indeed occur. In order 
to give convincing explanations, we should provide more sophisticated calculations and discuss 
the competition of possible ground states from unbiased viewpoints. Since the competition of 
possible instabilities often appears in many body physics, establishing those analysis provides 
understandings not only on dipolar systems but also on many physical open problems. Espe- 
cially, it will give deeper insights on the high-density nuclear matter which are important to 
understand neutron stars. 
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Chapter 6 


Summary and perspectives 


Main purpose of this thesis is revealing properties of many-body fermionic systems based on 
nonperturbative and unbiased method of quantum field theories. In order to realize it, we use 
the functional renormalization group and discuss its approximations in a systematic way to 
clarify physical origins of interesting phenomena. 

Chapter |2| is devoted to review methods of the functional renormalization group, which 
is our main tool to study many-body fermions. We introduced the Wetterich formalism and 
discuss its properties by comparing with other versions of FRG. We also consider optimization 
of the renormalization group flow proposed by D. Litim. 

In chap {3} we reviewed the Landau Fermi liquid theory based on FRG. We have learned that 
the Landau Fermi liquid theory is justified as long as we assume that low energy excitations 
are fermionic quasiparticles. We also learn how possible instabilities of the Fermi liquid appear 
from the RG point of view. The idea established here is a foundation in analyzing the BCS 
regime in the BCS-BEC crossover and dipolar fermions. 

We applied FRG to the BCS-BEC crossover without any bosonization in chap [4] Starting 
from scattering problem of the model in the vacuum, we consider superfluid phase transition 
using the RG flow in normal phases. In scattering problems, we considered the atom-dimer 
scattering in detail, and we obtained the atom-dimer scattering length without solving integral 
equations, which is not rigorous but improves the tree approximation significantly. In order to 
calculate thermodynamic properties with our formalism, we derived some formula to calculate 
number densities n using FRG, and also derived the Thouless criterion. In that argument, we 
discussed the Ward-Takahashi identity in a general way and showed that it can be understood 
as an integrated flow of the corresponding correlation function. 

We considered the BCS and BEC regimes separately, and as a possible connection between 
those two in leading order approximations we reproduced the NSR theory. In the BCS regime, 
we have seen that high energy physics is controlled by vacuum physics and low energy physics 
is described by the Landau Fermi liquid theory. By taking into account possible corrections in 
intermediate energy scales, we have reproduced the Gorkov and Melik-Barkhudarov correction 
in a natural way. In the BEC regime elementary excitations are bosonic dimers, and in order to 
describe the physics based on fermionic theories, nonperturbative description is required even 
at the lowest order approximation. Based on knowledge of the RG flow in the vacuum and 
separation of energy scales of many-body physics from that in the vacuum, we evaluated the 
fermion self energy correction so as to reproduce the bosonic picture. 

In chap[5| we discuss possible instabilities of dipolar fermionic systems using the scaling 
ansatz of the Landau Fermi liquid theory. We will find that the ansatz is justified in weak 
coupling regions and that °P, superfluid is predicted from our calculations. We also reproduced 
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possible instabilities in particle-hole channels suggested by T. Sogo et al. which corresponds 
to Pomeranchuk instabilities in the usual Landau Fermi liquid theory, and gave some discussion 
on relations of our result to previous studies. 

Throughout this thesis, we have seen that FRG provides a very powerful tool to study 
nonperturbative quantum field theories, but also that it is not yet completed and requires 
improvement especially for fermionic field theories. Separation of energy scales provides very 
important ideas to study physics, and establishing FRG is directly related to polishing up 
techniques to realize it in practical computations. Indeed, it is applicable not only many-body 
physics but also few body physics, and it may open a new way to analyze interesting physics 
clearly and intuitively. I believe that FRG will provide more powerful means to study very 
wide range of physics in the future. 


Appendix A 


Quick derivation of the l1-loop RG 
expressions 


In this appendix, we derive the abridged way to derive 1-loop expressions of the Wetterich 
equation and the composite operator flow equation (2.2.19). Of course, we can obtain 
them by expanding the both sides of the flow equation in terms of field variable since we 
already know their nonperturbative expressions, however it becomes very complicated as the 
order of vertices increases. We here would like to obtain the quick derivation of them using the 
knowledge of Feynman diagrams in the perturbation theory. 


A.1 Feynman rules for the flow of vertex functions 


Let us rewrite the Wetterich equation (2.2.11) as 
aT sly] = ds STrIn [rẹ Ny] + Ry. l (A.1.1) 


where Ôp, acts only on the k-dependence of the IR regulator Rẹ. In order to use Feynman 
diagrams, we decompose the effective action into the quadratic part and the others. Let us 
denote the decomposition as 


1 € 
Lely] = zr Gyo + Vilel, (A.1.2) 


where V;[(y] does not contain quadratic terms. The inverse propagator G" is often written as 
G7! — Xp, where G, is the free propagator and ©, is the self-energy correction at the scale k. 
In order to derive Feynman rules, we rewrite 


Est A re? + Ry. = EST In (Ge + Ry) + VO Ip 


1 iyor" 
= 5STrn(G," + Ry) — a ee Ee (org + Re) Vg (el) 


The first term is independent of field variables, and for the most purposes we may simply 
discard that part. Eq.(A.1.1) becomes 


-ôT lp] = SD =i (Gy + Ry VP Tel) (A.1.3) 


n>1 
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We can easily interpret the right hand side of Eq.(A.1.3) diagrammatically. To simplify the 
interpretation, we at first neglect the derivative ô% in the right hand side of Eq. (A.1.3) for a 
moment. Then it is just the summation over all the possible 1-loop diagrams with the following 
rules: 


1. Write down the all possible 1-loop diagrams with the interaction V;,|y] with the given 
external lines. Regard that all the external lines are amputated. 


2. To each vertex, attach the corresponding coupling in —V, |v]. 
3. To each internal lines, attach the propagator (G3! + R)! = (G3! — Ep + Rp). 


4. Take the super-trace. In the momentum space, take the trace of the flavor indices and 
d 


attach the loop integration f 7: For the fermion loop, attach the factor —1. 


d 
(27) 
5. The symmetric factor appears as the same way in the usual Feynman rule. 


After all these procedures, we get some analytic expressions and we get the result by taking 
the derivative of the k-dependence of Rọ. We should not forget the negative sign in the left 


hand side of Eq.(A.1.3). 
A.2 Feynman rules for the flow of composite operators 


In the same way, we can derive the quick rule for the composite operator flow. Again, we denote 


Eq. (2.2.19) as 


F2) 
rOl TR 9 . (A.2.1) 


Substituting the decomposition given in Eq.(A.1.2) into Eq.(A.2.1), we find that 


= ~ 1 
Ont lo] = O15 8Tr 


Yo (Get + Re) EVP tel) (Gat + RPI . (A22) 


n>0 


In this case, we should write down all the possible 1-loop diagrams which contains one of the 
couplings in J,|y] only once. The other rules to get diagrammatic expressions are the same 
with those of vertex functions. 

In this case, the expectation value itself of the composite operator is important for many 
cases. Therefore, we usually should not discard the field independent terms. When the ul- 
traviolet divergence appears in the continuum limit, we have to apply the composite operator 
renormalization. In general, operator mixing appears, that is, not only the multiplicative renor- 
malization but also the subtraction of lower-dimensional operators is necessary to get a finite 
result {10}. 

Finally, let us derive useful formula to get explicit forms of correlation functions. We often 
have to calculate the quantity Wk nm using the 1PI effective action for obtaining the explicit 
form of the correlation functions Talo]. Of course, when we decide that we calculate it using 
the flow equation (2.2.19), the calculation for the explicit form is unnecessary. However, it 
sometimes provides the way to check the consistency of the calculations and also gives the 
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physical insights for the problems. Therefore, let us write down diagrammatic rules to obtain 
correlation functions. 

For simplicity, we assume that the original correlation function I[J, | does not depend 
on sources J, that is, it is a usual composite operator J[¢]. For this case, the formulae 


(2.2.16]2.2.18) give 


a ô 
Tele] = 1 le + Grlo] - z| 1, (A.2.3) 
where 
7 W[J] 
J=R E +e Ri 


Therefore, it suffices to give the diagrammatic rule to calculate G;[y]. Generally, we can say 
that this quantity is the functional inverse of the second derivative of the effective action. Again 
using the decomposition (A.1.2), we get 

_ ÔLÔRIk 


Gly) = jose TT (Gy! + Re) + VP l. (A.2.5) 


Formally, we can take the inverse of the quantity (A.2.5) to get 


Gale] = (Gg + Re) — (Gal + Be) VO le] (Gat + Ba) (A.2.6) 
+ (Gig! + Ba) VKP lel- (Ga! + Ba) VP le] (Get + Re) 


Diagrammatically, Eq.(A.2.6) tells us that we should sum up the unbranched tree diagrams 
with the full propagator (G,* + R)! and with effective couplings —V;,[p] at the scale k. 


98 


Appendix A. Quick derivation of the 1-loop RG expressions 


Appendix B 


Continuum limits 


In sec [2.3] we have derived the Polchinski equation, which represents the renormalization group 
flow of the Wilsonian effective action. By introducing the explicit ultraviolet cutoff A, in the 
original theory, we are free from UV divergences in all calculations. In the relativistic system, it 
is important to control short distance physics without knowing its details and for that purpose 
we require renormalizability for relativistic quantum field theories. In this appendix, we shall 
take the continuum limit A, — oo and explicitly see how perturbative renormalizability is 
ensured in this formalism. For simplicity, we restrict ourselves to consider the ¢*-theory in the 
Adimensional Euclidean spacetime (the ¢j-theory) and prove its perturbative renormalizability. 
The proof with FRG was first given by Polchinski (9| and it was greatly simplified by Keller, 
Kopper, and Salmhofer [94]. 


B.1 Continuum limits of quantum field theories 


The nonperturbative realization of the renormalization of the quantum field theory is sketchily 
summarized in Fig/B.1]| The critical manifold containing the fixed point is the set of the points 
which flow into the fixed point in the RG flow; around the fixed point, the critical manifold is 
spanned by the infinite set of irrelevant and marginally irrelevant operators. 

Let us consider the flow starting from the bare action which is sufficiently close to but not 
in the critical manifold. As shown by the dotted line in Fig/B.1| the effective action flows into 
the vicinity of the fixed point and then runs away along one of the relevant directions. As a 
result, it will reach a trivial theory which represents an infinitely massive theory. 

We need a trick in order to get the continuum limit of the theory with finite masses. At 
first, we should tune the bare action back towards the critical manifold in order to prevent the 
theory from becoming infinitely massive. Simultaneously, we should rewrite physical quantities 


Critical manifold 


UV fixed pointy Renormalized trajectory 


Infinitely massive 
theory 


Figure B.1: Schematic view of the RG flow in the space of bare couplings. 
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in terms of renormalized quantities. After these renormalization procedures, we take the limit 
that the bare action touches the critical manifold, and then the RG flow splits into two parts. 
The first one goes into the fixed point, and the another one emanates from the fixed point along 
a relevant direction, which path is called a renormalized trajectory. In terms of renormalized 
quantities, the far end of the renormalized trajectory obtains a finite limit. In this sense, the 
continuum limit is said to be defined at a UV fixed point. 

Using this idea, we will explain what is expected as a consequence of renormalizability of 
the ¢j-theory following the discussion given by Polchinski (9). We give the bare action, or the 
Wilsonian effective action at the UV cutoff scale A,, in the next form: 


bake lo] = 2 K(p?/A2) P-pPp + Sint, Ka, lo] 


1 [Am > AZ 2 Ào 4 
= : Kag r + f x (Ao + -y (On) = a? ) « (B.1.1) 

We would like to integrate out high-momentum modes of a field, so that we reduce the cutoff 
A= A, in Eq. (B.1.1) to a much lower scale A = Ag. Now we have known that, by changing the 
effective action Sint, x, |$] according to the Polchinski equation (2.3.18), correlation functions of 
low-momentum modes are unchanged. 

Even though the action might start with a simple form Spare, at lower scales Sint, contains 
the nonrenormalizable interactions and becomes very complicated. However, in the infrared 
region A}, < A?, the effective action Sint, K ie will be strongly attracted to the set of renormalized 
trajectories R, which is totally parametrized only by the relevant and marginal couplings, 
i.e., the renormalizable couplings In our case, by keeping the Zə symmetry ¢ +> —¢, the 
submanifold R is three-dimensional, where the dimension 3 is determined by the number of the 
renormalizable operators $7, ¢07¢, and ¢*. 

We parametrize the effective theory by effective vertices V2,, corresponding to the expansion 


(2.3.19): 


ae 
Sint, Ka 10] = `> (Qn)! Pp i Pron” (27)*8* (pı +++++ Pon) Van (A; P1,°+* Pon). 


n=1 Pils »P2n 


In order to define the perturbation theory, we introduce a formal variable g. First of all, Sint, is 
taken to be a formal power series in g, and hence 


Sint Ka lO] = Se oe [o]. (B.1.2) 
r=1 
At A = Aj, we set for r > 1 
7 Am?, AZo ae 
a= jas í 7 dg? 4 > (8,0)? + ao l (B.1.3) 


We should notice that the formal power series (B.1.2) starts from r = 1 at any scale A, and 


furthermore S E A 


is an even polynomial of degree < 2r +2. The RG flow is given by the 


‘Tn the perturbative argument, distinction between relevant and marginal is unnecessary. But the distinction 
will be important if we consider the continuum limit in nonperturbative arguments [6]. See also sec[2.5] 
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Polchinski equation (2.3.18). The boundary conditions for the RG flow are mixed type: at 
A = A, the theory is required to become ¢4-theory and at A = 0 the renormalization condition 
is given for the relevant parameters: From (B.1.3), we have 


(a) VS (A,) =0 (for 2n > 6), 
b.c. at A= A, (6) Vo (A.)(p1, P2, P3) = Aotr); (B.1.4) 
(c) VPA) (p) = AZ (Ao) p? + Am?) (Ao), 


which imply that TALA (Ao) = 0 for 2n+w > 5. Here we have omitted the last argument of the 
vertices since it is determined from the other arguments due to the momentum conservation. 
We should remind that Amọ) (Ao), AZin (^o), and Aor) have not been prescribed yet and they 
should be determined by the boundary conditions at A = 0, which are given by 


CR a m a p®) = 0%, 
b.c. at A=0 (0). G nP (A = o e =A, (B.1.5) 
(d VPUa=00 =m, 


where p? (i = 1,...,5) specify the renormalization points. The usual Bogoliubov-Parasyuk- 
Hepp-Zimmermann (BPHZ) renormalization condition can be obtained by setting p? = --- = 
pe = 0 and ri = 6,1 and AZ = Ain = 0. 

In the low energy scale A < m, the mass term changes the scaling of the fields drastically. By 
using another renormalization scale Ar > m, we need not consider such complicated situation 
since we can stop the RG flow at A = Ar. For that purpose, we rewrite the condition (B.1.5) 
into the boundary conditions at A = Ar: 


(a) Vi? (An) (Df, pF, p) = A + O(r), 
b.c. at A= Ap (b) Op, Op, VS (Ar) (PR) liu = 28u (AZẸ +O(r)),  (B.1.6) 
r) 


(c) VS? (Ar)(p2) = Amg? + O(r), 


where O(r) stands for the contribution of a finite number of connected, amputated Feynman 


graphs of order r whose external momenta are p?, pł, pł or pł or př and whose propagators 
K(p?/A?) 
“ae © 
In order to estimate effective couplings, we introduce a semi-norm on vertex functions. Let 


F(A, {pi}) = f(A; pı, -++ pn) be a sufficiently regular function of momenta {p;} defined in the 
theory with an effective cutoff A. We require the regularity of the function f(A; {p,}) so that 
the following definition of semi-norms makes sense: for any 7 > 0 


are 


Ə” F(A)||, = max max sup Oper ++ + Opie f(A; pi, Pn). (B.1.7) 


{iise iw} {May Hw} \pi|<max{2A,n} i1 


This definition is useful. Let us take the cutoff function Ka(p?°) = K(p?°/A°) so that it im- 
mediately damps above p? = A? so as to behave like the step function. Then K(p?/A?) = 
—A0, K (p?/A*) behaves as the delta function which has a peak around p? = A?, so that 


wh? K (p?/A*) 
p? + m2 


< C, and ||0 


4 A2K 2 / A2 
J dtp (p*/A*) < DA”, 


(Q27A)4 pP +m? 
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that is, the scaling behavior expected from the mass dimension is correct. 

We shall derive the inequalities which are necessary for the proof of perturbative renor- 
malizability. At first, let us rewrite the Polchinski equation (2.3.18) in terms of dimensionless 
quantities Ay, = A?”~*V2, (see Eq. (2.3.20): 


A?K (q?/A2) 


Pk me 2o) MA; a {Pi }) 


(- A + 2n — 4) AnA; {p}) =- X. Aow(A; {pa} 9) 


partition 


1 f dtg AK (@P/A) An+2(4; q, =q, {pi}). (B.1.8) 


2J (Q27A)* +m? 


We can readily get the following inequality for arbitrary n > 0 


(- AX +2n- 1) aman (B.1.9) 
n 
(r,n,w) 

Zeuas | OAS oly + S > AT ||" AL (A) [Ill AG Ray (Aly | 
ri kw; 


where the summation is taken over all the combinatorics for tree diagrams, i.e., 7’ runs over 
1,...,r—1, k runs over 1,...,n, and {w1, we, w3} runs over all the combinations for w + w2 + 
w3 = w. Here, Cyn, is independent of A and 7. We should notice that this is a very suitable 
form in terms of the perturbation theory. The first term in the right hand involves only lower 
orders with respect to the coupling constant and the second term only involves the term with 
larger n. Since 2n is the number of the legs, at each order of the perturbation An (^A) vanishes 
for large n. Thus we can expect that the 15 i will work. 

As we ig see in the next section, Eq.(B.1.9) is sufficient to prove the boundedness of the 
norms Ae (A R)\|(2A,n) as Ao —> œ. Peo in order to prove the convergence we need the 


following inequalities involving ôa, 0y AM (A (A). 
For irrelevant vertices, boundary conditions at A = A, is well-known for us. Therefore, for 
cases 2n + w > 5, apply ôa, to the equation: 


^o ds OY ae, oS wmf hf Od hie) 
| = (--2) s Aon(s; {pı}) = f dss o Pom An+2(5; q, —q, {p1 }) 


s°K (q?/s*) 
q + m2 


= Aah 


partition 


Aag(n—k-+1) (8; 9, (oo) , (B.1.10) 


which can be easily derived from Eq.(B.1.8). Taking the norm, we get for 2n + w > 5 


[Aan 0" ARA) (B.1.11) 
n 
(r,n,w) 
T | O"Amo(o)Iln + D> As MIO" Age (Ao) IInllO" a 
r! kw; 
(r,n,w) 


+h, n, ~ dss?" | [3a 3" AR l)la + YO s aa 0A OAOA day (9) ln 


r!’ kwi 
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Again, constants Ch, nr and Chn, are EL os of Ao, A, and 7. On the other hand, one may 


take the integration interval in Eq.(B.1.10) as [Ar, A], which will be more appropriate choice 
for the relevant vertices. Applying ar a On,, we obtain that 


AIA OPARA pi, +++ sPan)| — [A On 0" AR (Ans Pr, “++ Pan) (B.1.12) 
(r,n,w) 

i I dss?" | [3a 3” Arnio(s)llar + >) 573a 0A (llall Age llm | 
AR r! kwi 

where M = max{|p;|;i = 1,--- ,n — 1}. Now we have finished preliminaries and then we can 


prove perturbative renormalizability of the Euclidean massive @4-theory. 


B.2  Perturbative renormalizability of the ¢/-theory 


Here we follow the proof of the perturbative renormalizability of the ¢j-theory given by Keller, 
Kopper, and Salmhofer [94]. The rigorous statement of the theorem is as follows: 


Theorem B.2.1 (Perturbative renormalizability). Consider the field theory defined by 
the action and define the effective action Sim x, (^o) as the solution of the Polchinski 
equation with the boundary conditions and (B.1.5 (B.1.5). ee order by order with 
respect to a formal variable g in the perturbation theory, the limit lim VO(A = 0; p1, +++ ,Pn—1) 


exists for allr and n. Again, we have omitted the last argument of the vertices thanks to the 


momentum conservation. 
Moreover, lim yo) (A = 0; p1, +- ,Pn—1) 18 a smooth and polynomially bounded function of 
o— 00 
P1,*** ;Pn-1- 
Let us add a convenient notation. The symbol P In(z), where z = = iV with A € [AR, Ao] 
and hence z > 1, represents some possibly polynomial in ln(z). Also Py represents some 


polynomials in 7. Each time it appears, the coefficients of the polynomial can be different as 
long as they does not depend on 7 nor on A nor on Ag. 


Lemma B.2.1 (Boundedness). For any ņ > 0 and for A €E [Ap, Ao], we have 


[22 AS? (Ally < AT’ Pr- (Pn + | EPn x] ; (B.2.1) 
Proof. We will prove it by induction. The induction hypothesis is that Eq. (B.2.1) is true for 
r < To Or Nn > no for any w > 0. The induction starts since Eq. (B.2.1) holds trivially true for 
r= or 2n > arse 2, 
For the irrelevant vertices, i.e., for gr AÇ?) (A) with w + 2no > 5, we employ the boundary 
condition at A = A, to get 


òns | 


A anA j= Ve sd, (searag) | < a d |[sð; (820e AS)(s)) 
n 


n 


Apply Eq.(B.1.9) to the integrand, and then the induction hypothesis gives 


A 
° ds s 8 A 
4—2no QW Z 4—2Nno — w | o 
IA ð; ANA Illy < f P (Pn È | pmi), 
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k s Aa as ie 
The integration by parts gives dss ~P In — =A | Pin — + — | P, ln = 
Ar Ar Ar 
A b=1 


for a > 2, and then we can obtain the result and arrive at Eq.(B.2.1). 

Now let us consider the relevant vertices 2n, +w < 4. The boundedness of these vertices 
will be derived from the finiteness condition at A = Apr. Let us first discuss 2n, = 4 and w = 0. 
Then 


A A 
(AL? (A) — AY”) (Az) (tp? PS) < | ds|ð,A{” (s; př, py, p3)| < f ds||0,AY”'(s) lav, 

R R 
with M > max{|p?| : i = 1,2,3}. According to the induction hypothesis, Eq.(B.2.1) holds for 
AS) (A) with r =r, and 2n > 6 or with r < ro and n EN. 

Now we can discuss the detail for the term O(r,) in part (a) of the boundary conditions 
(B.1.6). Since it is constructed with the finite number of the diagrams of order r, with lower 
order vertices AW) (Ar), which can be already bounded by Eq. (B.2.1), and therefore we can 
obtain that 

|AY? (Ani pt, pF, p8)| < ¢, 


where c is a constant which is independent of A and A,. Therefore, 


A A 
T Ao 
IAL) (A: pE, pË, p®)| < e+ f ds||0.AS(s)Ila < c+ f dae (Pin ame oT, ) 
AR AR AR Ao AR 


< (Pin : + A p] e), 


Ar Ao AR 
Therefore, we have obtained the bound for the renormalization point (p1, p2, p3). For the general 
momenta, the Taylor theorem gives 


3 4 1 
AÇO (A; qr, q2: 93) = AÇ (A; pit, pF Ps) + YO YO PÈ, — ai) f dzôy „AY? (A; q + 2(p® — q)), 
0 


i=1 p=1 
and then taking the semi-norm ||- || we can obtain that 

A A A M A A A 
AQ(A)|] < Pi in (A EA p ti 
Ar). = a Pin Eon do Pa 

A A A 

< Pl In = 

< P a 


This is nothing but Eq. (B.2.1) for 7 = M, and hence it is true for any y > 0. 

Next, let us evaluate the case with 2n, = 2 and w = 2, i.e., estimate the bound for 
PAY (A). The assessment is analogous to that for the case 2n, = 4. By estimating the 
finite number of the Feynman diagrams contributing to O(r) in the part (b) of the boundary 


conditions (B.1.6), we get PAs) (A Ri DA < œ. The Taylor theorem gives 


1 
PAY? (Ar : pg) = BAS? (An: p=0)+ > pe, f dad, 02 AS? (Ar; xp). 
0 
H 
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We should notice that OP Al) (Ap: p = 0) only has the term proportional to 6,, due to the 
Lorentz covariance, so we can find the boundedness for the first term. Applying the bound 


a2 AS” (Az)|| we can obtain that Japs” (An; p4) 
dent of A,. Thus, 


< c" for some constant c” which is indepen- 


[A202 AS) (A; p) 


A 
< Ake" + | ds | 
AR 


ðs (5202A (8) 


M 


where M > |p|. Using Eq. (B.1.9), we can find the P In-type bound for AEA (Ae pit) 
Again the Taylor theorem gives the bound for any momentum and then we can arrive at the 


bound (B.2.1). 


For the case 2n, = 2 and w = 1, the Taylor theorem gives 
1 
BAPA p) = BAYA =0) + Dre | and, 0AL? CA: p) 
0 
H 


and the first term in the right hand side vanishes due to the Lorentz covariance. Therefore, we 
can find the desired bound for BASPA] 

Finally, we consider the case with 2n, = 2 and w = 0. By the renormalization condition, 
|AS? (Ap: pB)| < c” and then 


A 
PAM (Arp) < Ake + f as| 
AR 


Os (s°4y""(s)) | | 


’ 
M 


where M > |p|. We again obtain the P In-type bound for la2age? (A; pÈ) f The Taylor theorem 
gives the desired result. 


Now we can show the convergence of the effective actions in the limit A, — oo in the similar 
way. 


Lemma B.2.2 (Convergence). Assume that for n 2 0 and A E [Ap, Ao] 
| 
Ao 


< —Ww 
(SAPs 
OnI ARA|| < Asta Pn Fe (B.2.3) 


aL ASA) (B.2.2) 


R 
Then we can find that 


Proof. The statement will be inductively proved, and the induction scheme is precisely what 
is done in the proof of the previous lemma. Again the induction begins since the statement is 
trivial for r = 0 or 2n > 2r + 2 for any w > 0. 

Given no, consider separately the cases w + 2no > 6 and w + 2n, = 5. For w + 2n, > 6, 


Eq.(B.1.11) and application of Eqs.(B.2.2|B.2.3) to the right hand side of Eq.(B.1.11) give 


3—2no—w Ao Sa 3—2no A —2 .—wt+l1q/ Ao 
| < 32h wD Jy 22 4 | dss? A727 Ytp Jy = 
n Ar A A 


AA, 0 Amn (A) 


R 
A A 
— A 3—-2no-—w o —2 5—2no-—wW _ A 5—2no—w 1 o 
= As Pins + As? (Ag A E 
< hae In Ao 
AR 
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Therefore, we have obtained the result for 2n, + w > 6. For 2n, + w = 5, the first inequality 
above holds true but the following evaluation is different. So we get 


3—5 No ^o -2 3-5+1 py), Ao 
< Ay °PIn— + dsA,“s P In — 
n A A A 


[A0 aa 0" AS). (A) 
R 


R 


-2p Ao 
< AZ P” In ie 
Thus, we have obtained the desired bound for the irrelevant vertices. Notice that this proof is 
valid only for A < A, but due to the continuity the limit A —> A, justifies that the results hold 
true for any A €E [Ap, A]. 

Now let us consider the effective vertices with 2n, +w < 4. We should notice that 
[3a AY (An; pË, pË, p)|, [a024 (An; p?)| and |31, AS” (Ar; p®)| do not exceed A5?P In #2, 
since the contribution comes from the finite number of connected amputated diagrams contain- 
ing the vertex differentiated by ôa, and because the tree diagram vanishes due to ôa, AZ = 


Hj Ame,” = OLN =0. 
We can find that Eq.(B.1.12) gives the appropriate bounds for [3a AS? (Ap; p®, p®, p®)|, 
(3a, 02 AS) (Ap: p®)| and }O,, AS? (Ap; pË). The applications of the Taylor theorem gives 


Eq.(B.2.3) for these vertices. O 


Since the assumptions in the second lemma are ensured from the first one, we get the proof 
of the theorem of perturbative renormalizability. Let us give some concluding remarks on this 
appendix. 

In the proof, the notion of fixed points seems to be hidden although we have pointed out 
that the continuum limit of the quantum field theory is to control bare actions so that they 
approach to a fixed point along a renormalized trajectory. This is a trick of the perturbative 
renormalizability and it is related to the triviality problem of the 4 theory. 

In these proofs, we have implicitly assumed that we can use the Gaussian fixed point (GFP) 
which describes the massless free field theory as a UV fixed point of the ¢{ theory. Indeed, 
the theorem suggests that effective couplings obey the scaling law predicted from canonical 
dimensions up to some small logarithmic corrections, which give small anomalous dimensions. 
At each order of the perturbation theory the GFP can thus play a role of a UV fixed point and 
the irrelevant couplings are controlled by the renormalization condition. 

However, when we discuss the (Borel) summability of the perturbation theory the situation 
becomes totally different. The 1-loop RG flow of the ¢*-coupling suggests that it is marginally 
irrelevant around the GFP, and then the theory is not asymptotically free. The appearance 
of the Landau pole in the UV regions prevents us from taking the continuum limit in a non- 
perturbative way. This means that if we would like to use the GFP as a UV fixed point then 
the theory should be free and can contain no interaction terms. In the case of the ¢4 theory, it 
is believed that there exists no other UV fixed points and we should regard it as a low-energy 
effective theory with some UV cutoff. 

In the case of non-Abelian gauge theories, we know that they have asymptotic freedom. 
Therefore, it is widely believed that non-trivial fixed points of non-Abelian gauge theories can 
be realized with appropriate fine-tuning of the renormalizable couplings. 


Appendix C 


Angular momenta 


In this section, we give a brief review of the theory of angular momenta according to [95||96}. 
The purpose of this appendix is to specify the convention on the theory of angular momenta, 
since there exists several ones used in common. Let J = (Js, Jy, J.) be an operator of angular 
momenta, satisfying the commutation relation 


We denote the simultaneous eigenvector of J? and of J, by |j, m), with m = —j,—j +1,- ,j: 


Jlj,m) = jG + Dj, m), j-i, m) = mlj, m). (C.0.2) 
Here we have normalized the eigenvector |j, m) as 
(ji, moəlj2, m2) = Osis) ices (C.0.3) 


and the phase convention of states |j,m) with different m is chosen so that they satisfy 
Jlj m = /(j F m)(j m+ 1)|j,m + 1) for raising and lowering operators J} = Jy + iJy. 


C.1 Clebsch-Gordan coefficients and spherical tensors 


Let us consider the coupling of two different angular momenta Jı and J2, and denote the total 


angular momentum J = Jı + J2. We define the Clebsch-Gordan coefficients ar Biama DY 
|j, m) = ` Cn iame ™M1) 8 |j2; Mə}, (C.1.1) 
my ,™m2 
or, denoting |ji, j2; Mmi, M2) = |j; Mm) 8 |j2, M2), we have C”? jama Z (Ji; J2; Mi, Malj, Mm). 


There still exists ambiguity of sign, since the normalization (C.0.3) does not tell relative phases 
between states with different j’s. This ambiguity is resolved by requiring that (j||Ji||7 — 1) = 
—(j||J2||7 — 1) are real and positive, where (j||Jj||j’) are reduced matrix elements defined by 


Eq. (C.1.18). 


It should be noticed that the Clebsch-Gordan coefficients (oe jam-mı tj,m, are elements 
of a unitary transformation from the tensor product basis {|j1, j2; M1, M — ™1)}im, to the irre- 


ducible basis {|j, m) } 


jm jim — 
Ba Giimiam-na Cremen- = Õjj’. (C.1.2) 
mı 


1Regard these are the basis for the Hilbert subspace specified by the conditions J? = jı (j1 +1), JZ = jo(j2+1) 
and J, = Jj, + J22 = mM. 
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Figure C.1: The Euler angles 


Therefore, we can invert the formula (C.1.1): 


lji mı) 8 |J2,m — mı) = 2 Cis Jam-ma J; My- (C.1.3) 


In this phase convention, the Clebsch-Gordan coefficients satisfy following symmetries: 


gims = (— acne aom -m3 = (— oe js qisma -= (=p mm 253 +1 j2—=m2 
jımı;j2m2 —M13J2,—M2 Jj2m2;jımı 2ja +1 J1M1;J3 = M3 
(C.1.4) 


Before introducing spherical tensors, let us consider the rotational matrix and its transfor- 
mation properties. We denote the Euler angle ay as given in Fig/C.1]| The corresponding 
rotational operator exp(—i0n - J) is given by 


R(aBy) = exp(—iaJ,) exp(—iBJ,) exp(—iyJz). (C.1.5) 
We denote matrix elements of the operator R(aBy) by Di (apy): 
Dimlab) = (jm’|R(aBy)|jm). (C.1.6) 


The unitarity of the representation suggests that Dİ, (—7, —ß, —a) = D” (a, 6,7). Let us 
apply the rotational matrix to the both sides of Eq.(C.1.1), then we get the expansion 


moby) = `> C J2me2 Cs Foun ten, (aby) Di; ma (a By). (C.1.7) 


mı,m2 mi, ms 


Inverting this relation, or using Eq.(C.1.3) instead of Eq.(C.1.1), we obtain the relation 


Di, m1 (aby) DË, gma (aß) = 2 O j2mh C mis nae” ois (a8). (C.1.8) 
Schur’s orthogonal relation suggests that 


dQ Recs 
=p D? = J192 ee L 
/ Fg Pimms (181) D ym (OBI) = Fema g Imm Smm, (C.1.9) 


Here we denote the integration over the solid angle by dQ, expressed with the Euler angles as 
faQ = + ~ da fe sin 8d8 i dy. As a special case, we obtain 


20 T 
f ay f sin 0d0Y (0, P)Yrm (Oy) = Ôl mm’, (C.1.10) 
0 0 
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since D!.(a@80) = ,/ =Y; (0p). From (C.1.8), we get the composition rule: 


2+1 


2 a (2l + 1) 


Yn any (8, Q)Vin, ma ( (21 + 1) 


Cima Free T p). (C.1.11) 


The irreducible spherical tensor {Tr}. of rank L is the set of 20+ 1 operators Tgm with 
M=-L,-L-+1,--- , L satisfying 


L 
RTpmR = X` Dipulaby)Tim, (C.1.12) 


M'=-L 
where R = R(a,ß,y) is the rotational operator. Equivalently, we can characterize the irre- 
ducible spherical tensor using the commutation relations: 
[Ja Tem] = V (L F M)(L +M +1)Tt, m+, (C.1.13) 
[J., Ttm] = MTr m. (C.1.14) 
Due to the linearity of the condition (C.1.12), or of (C.1.13}C.1.14), the set of the irreducible 
spherical tensors of the same rank forms a vector space. 


We can form another irreducible tensor from two irreducible tensors Ty, m, (A1) and Tr, m (Az) 
as follows: 


Tim (Ai, A2) = `> CE M TaM taii (Ai) Trai An): (C.1.15) 


Mı,M2 
Here we have denoted all the other arguments of the irreducible tensors by A;. Especially when 
Lə = Lı, we can form the rotational invariants: 


(=)! 2L + IT oo Ar Aa) = J(=) Tran (4)Tr,, -m (A2). (C.1.16) 


As aspecial case of this formula, we can obtain the addition theorem of the spherical harmonics: 


P,(cos 0) = ma Ta FiO Yin (01, 91) Yim (02, p2). (C.1.17) 


The Wigner-Eckart theorem gives the eos quantum number dependence of matrix 
elements of an irreducible tensor: 


('m'|Trulim) = Chom 'lITeIId), (C.1.18) 


where the second factor (j7'||77||7) is called the reduced matrix element, which is independent 
of projection quantum numbers. 


Proof. Using the commutation relation (C.1.14), we find that (m — m — M)(j’m'|Tru|jm) = 0. 
Therefore, the matrix elements (j'm' |T; m| jm) vanish unless m’ = m + M. Another commuta- 


tion relation (C.1.13) gives 
SVG Em Fm! +I m FU Tralim) — JG FMG Em +1) Truli m1) 
= /(L=M)\(L+M + 1)(j'm'|Tr m| jm). 


The same recursion relation holds for the Clebsch-Gordan coefficients Gr > m Which satisfies 


the same constraint (m’ — m — M jan zm = 0. Therefore, the dependence of (j’m'|Trar|jm) 


+ 7 . jim! 
on projection quantum numbers is the same as that of Cim Lm 
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C.2 67 symbols and 9j symbols 


Here we consider the composition of three angular momenta Jı, J2 and J3. For the uncoupled 
representation, we can diagonalize the 6 operators J?, J2, J? and Jiz, Jaz, J3,. There are three 
different ways to couple those angular momenta with a resultant angular momentum J and its 
projection J,: 


(1) J, + J2 = Jia, Jig+J3= J, 


(2) Jot J3 = Jz, J3 + Si = J, 
(3) Jp +J3 = Jg, Jig + Jo = J. 


Let |j1j2(J12)j3; 7m) be the state vector corresponding to the composition (1) above. Us- 
ing the Clebsch-Gordan coefficients, we can relate these vectors with the uncoupled basis 


lji, M1; J2, M2; J3, mg): 


ljj2(ji2)j3; jm} = x. O uume ` ciem malji, Mi; j2, Ma} Js, ma). (C.2.1) 


™12,™3 ™m1,™2 


Introducing similar notations, we can readily find the following: 


\j2J3(Jo3)I13 JM) = 5y C TE ` CH malji Mai; Ja, Ma; j3, M3), (C.2.2) 


mM1,m23 m2,M3 

Eu a a = jm J13™13 aa oa 

lja (j13)j2; jM) E ; nee ree : Crimi Jama lJi, M1; J2, M2; j3, M3). (C.2.3) 
mM13,mM2 M1, m3 


Since states in each coupling choice form a complete set, there is a unitary transformation from 
one scheme to another one. We define the Wigner 6j symbols, or simply called the 6j symbols, 


[a a ey 
J3 32 J23 


lijaljia)js; imjods Joa) i139’) = djy Sam! (ETH S Oja + jz + 1) l ae \ 


J33 J JIB 
(Q24) 


The previous factor of the symbol is introduced so that the 67 symbols are invariant under any 
permutation of its columns or under interchange of the upper and lower arguments in each of 
any two columns. For examples, 

e an 

bef 


a er yeas fas A 


According to the definition (C.2.4), the 67 symbols may be expressed in terms of the Clebsch- 
Gordan coefficients: 


a a 


jm J12™M12 jim! j23M23 
C; C: C: G- 


J12™12,j3M3Z “jımı, jamo “~ jımı,j23M23 “ j2m2,j3 m3 
M1, M2,M3,M12,M23 


= 55S ram (—)2 +8 4345 \/(Qir9 aa Dat 4 Jı J2 J12 \. (C.2.5) 
J3 J J2 


Since the phase convention are given in the definition of the Clebsch-Gordan coefficients, the 
definition of the 67 symbols or the Racah coefficients is complete. 
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Let us next consider an addition of four angular momenta Jj, Jo, J3 and J4. There are three 
different ways to couple those angular momenta which are not directly related to the product 
of the two 67 symbols: 


1) J+J = Si, J3 + Ja = Js, J2 + Ju. = J, 


(2) J+J = Jg, Jz + Ja = Ja, J3 + Ju = J, 
(3) Jı + Ja = Ju, J2 + J3 = J233, Jig + J23 = J. 


The following argument is similar to that was done in the discussion for the 6j symbols. 
For each scheme of the composition, we introduce the state vectors |j1J2(J12)J3J4(j34);9™), 
jaja (ji3)j2ja (J24); jM), and |j1jal jia) j2j3(j23); jM), respectively. Using the Clebsch-Gordan co- 
efficients, we can explicitly write down those states as the tensor products. For instance, 


lji j2 (j12)j3ja (i34); jm) = eae ek ee |jımı; J2Mə; Jams; jama) > 
j j jimisj jams;j 
M1,- M4 M12,M34 
(C.2.6) 
In order to realize the unitary transformation in the different composition schemes, we introduce 
the Wigner 9j symbols by 


ji h2 jie 
(91.92(J12)I3.J4(934); jml jaja (jis) jajal joa); IM) = Sjj Omm Ujrojrsjoajaa js ja ja (C.2.7) 
J13 Ja 3J 


with Iab..c = y/ (2a + 1)(2b + 1) --- (2c +1). The 9j symbol does not change its value under the 
transposition and under the even permutations of its elements. Under the odd permutations, 
the phase factor (—)? appears, where R is the sum of all arguments of the 9j symbol. 


C.3 List of formulae 


Here we list up the useful formulae to calculate the summation of the products of the Clebsch- 
Gordan coefficients. We pick up them from the textbook [96]. 
First we should recall the symmetries of the Clebsch-Gordan coefficients (C.1.4): 


Crab = (=j e a = rn ais (C.3.1) 


Cc a—a 2c + 1 b,—8 b B 2c + 1 a,—a 
Ce = (=) \/ 2b + T Carai, = (—) a 2a 4 [Ce 1b,8" (C.3.2) 


Sums Involving One Clebsch-Gordan Coefficients 


So C88 ig = V2a +1 X (-) 7g = (2a + L)bi0. (C.3.3) 


and 


Sums Involving Two Clebsch-Gordan Coefficients 


Caa bpl aa 5a = Dec! Oy! (C.3.4) 
ap 
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> Cadog C aa pp" = Soa! gp", (C.3.5) 
cy 
NO (20+ 1)08 Ce, oy = (2b + 1) 8.0539 (C.3.6) 
cy 
S (=) 82a + LICOS pC, oo = v/(2b + I) (Be + 1)dy,-1'5,-9". (C37) 


aa 


Sums Involving Three Clebsch-Gordan Coefficients 


c a+b+e abe 
ee Chea Get eds 1 ow e f y (C.3.8) 
evar p (24+ 1)V2a+1 nec a bec 
ys Cop, cy Crp, ecCaa afo = ae ee note Conte ade [i-7 (C.3.9) 
aps V2e+1 e f 


a— Q Cc eE Cc eE a b C 
YO) C e A mt (Cas TATE (2c T 1)(2f Dee fo f TE \ : (C.3.10) 
aps 


Sums Involving Four Clebsch-Gordan Coefficients 


For the simplicity of expressions, let us put Ias.. = \/(2a + 1)(2b4 1)---(2c+ 1). 


cba 
aa dô j 
Ds Cog, eyes, foCee BC] fe, cy T Haags ae an, jul ier j do, (C.3.11) 
Bree jg k 
a c 
Soc bB. cy Cee E O = Heca >, V2k + 1O a0 TK d f (C.3.12) 
Byep J 
Sums Involving Clebsch-Gordan Coefficients and a 67 Symbol 
Again, we put Igo... = / (2a +1)(2b4+ 1)---(2c4+ 1). 
abc 
5o) Tea ys, aC Fo, cy { e f d \ = Cc. pee. fo (C.3.13) 
c+d+f fe bac pm ee 
DE ie? Cr. f od e [= CaO (C.3.14) 


fe 


Appendix D 


Properties of p in the vacuum 


In this appendix, we calculate some properties of the symmetric part of the four point vertex 
function ['? in the vacuum defined by Eqs. (4.2.7) and (4.2.9). Although we can calculate 
T'?(P) for general momenta P and for general cutoff scales k analytically, we do not here give 
full calculations to obtain the formula for general cases. Instead, we will give calculation for 
special but important cases to find properties of the vertex. 

We use the unit 2m = 1 for simplicity of calculations. 


D.1 Calculation of I°?(P°, 0) 


Let us calculate the vertex function T$ (P) when the spatial momentum P = 0. Before starting 
calculations, let us find an expression of [?(P) for general momenta and for general cutoff 


scales. Using the formula (4.2.7) with Eq.(4.2.9), we find that 


1 1 
TE e S 2) a 


2 


1 1 -fE 
r$(P0, P?) 8nas J (2r)? 


T | on 


where A, is an ultraviolet cutoff. Since the second line of this equation vanishes in the limit 
A, > œ, we can eliminate the UV cutoff from this formula. As a result, we find that 


1 1 
FR) Pa TR N E 


i. 1 o 1 1 
(P, PA Bras J Om)? |2 P0 2p t B+ OP + Rl +B) + Rl S) 


(D.1.1) 
By putting k = 0, we recover the formula (4.2.10). 
Now, let us substitute P = 0 into Eq.(D.1.1), then we get 
1 1 ql 1 1 
3 = = | == . (D.1.2) 
T?(P°,0) 8ras (27)? [22 iP? — 2u + 212 + 2R, (0) 


Since we use the regulator R,(l) = (k? — I7)0(k? — 1°), we find that 


1 1 i 1 1 a di [1 1 
T$(P°,0)  8ras Jy (27) [24 iP? — 2u + 2k? , (2n)3 [2U iP9— 2u re 
D.1.3 
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Integrations in Eq.(D.1.3) are quite simple, and we immediately obtain 


iP? 
1 1 1 /k k3 /3 oa Bo k 
= — — 1— Ź tan! | ———— 
S (po 2 i 
T(P, 0) 8ras 2r 8T T Poy, 
V2 


2 i(P0— 2+ 2k? 


(D.1.4) 
If ag < 0, there exists no bound states and we can put u = 0. On the other hand, if ag > 0, 
then there exists a bound state and for the zero density condition we should put u = —1/az in 


the vacuum, as we have already remarked in sec/4.2| 
If IR cutoff scales k are much larger than other physical scales P®, u, then we find that 
1 Z 1 k 


~ TEE E D.1.5 
IS(P0o,0)  8rag 6r? I (amro (Pkn 


as desired. Let us consider the opposite limit k — 0, then we get 


iP? 
1 1 k VFTOTH 2 k 1 
~ 1 . (k= 0) (D.1.6) 


TS(P°,0)  8rag 4r? 8T T fa, 7 T?_(P°,0) 
V2 


That is, linear terms in k cancel in the Taylor expansion around k = 0. We can also easily find 
that cubic terms with respect to k also cancel, and then the convergence of the limit k — 0 is 
very fast. Let us consider the BEC case ag > 0, and we plot ag/T$(P = 0) in Fig/D.1] using 
Ba P18. 

Later, we will find that it is useful to introduce the quantity 0,,(1/T'?) = i~10po(1/T?)(P°, 0). 
Indeed, this quantity is necessary in Eq.(D.3.3). 


1 k? /67? 1 2 “i k 
Ou75 P°.0 _ ;P0 — 2 D2)2 ; 1 — — tan - 
p (P°,0) (iP? — 2u + 2k?) 327 i y T iP? y 
1 k 


= D.1. 
872 iP0 — 2u + 2k? Det) 


In the next section, we will use this to consider effects of explicit breaking of the Galilean 
symmetry. 
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Figure D.1: Graph of as/T'?(P = 0) when ag > 0. 
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D.2 Spatial momentum dependence of ed 


Now, let us discuss the dependence of the vertex function T(P) on spatial momenta P. If the 
cutoff is removed so that k = 0, then the Galilean symmetry is satisfied and the dependence of 
P can appear only through the combination iP° + P?/2. However, the Galilean symmetry is 
explicitly broken and this restriction does not hold in general cases. 

Instead of calculating the full dependence of rT? on P, we try to find the leading order 
expansion of 1/['?(P) in terms of P. Due to the rotational invariance, the expansion starts 
from P? and then we will calculate its coefficient. 

We introduce an expansion parameter s and then we consider the expansion of 1/T'?(P°, s?P?) 
around s = 0. Since the linear term in s must vanish due to the rotational symmetry, we have 
to calculate the second derivative of 1/T$(P°, s? P?) at s = 0. Using Eq.(D.1.1), we obtain that 


o? 1 
ðs? T(P, s2P2)| o 
7 J dêl o? 1 
(27)? Os? (P? — Qu + 9? + 2? lla) Rel — s¥)| 
-- [EË ð sP?—P.-(l+s2)0(k? — (1+ s£)?) + P- (l— s£)0(k? — (l - s£)?) 
(27)? Os ee | s222 12 4 R, (1 s?) + R(t — 8B)” - 
dèl P?(1 — 6(k? — I?)) + 2(P - 1)?6(k? — 1°) 
=— (D.2.1) 
(27) iP? — 2u + 21? + 2Rz(L)]? 
Therefore we should calculate the following integration 
o? 1 __p i dl 1 o1 k? 
Os? T$(P9, 2P?) | o p (2r)? (iP? — 2u + 212)? ` 6r? (iP? — 2u + 2k2)2 j ` 
(D.2.2) 
The integration in Eq.(D.2.2) can be evaluated as follows: 
/ dl 1 o 1 el dl I? 
p (2T) (iP? — 2u + 212)2 872 J, CET n 
Pe 
TE az 
2 
= i= ~ tan! - : 
327 H — H i a H 
1 k 
f : D.2. 
8r? iP? — 2u + 2k? ous 
Combining Egs.(D.2.2/D.2.3) and Eq.(D.1.7), we find that 
o? 1 1 
= P’*),———.. D.2.4 
Os? TX (P?, s2 P?) | o “Te (P90) ( ) 
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Therefore,this means that up to the quadratic order the Galilean symmetry is respected and 
we get 
1 1 P? 1 1 


= Oo TE e, D.2.5 
EPA Eo TR O e (poo (Dias) 


where ellipses represent terms smaller than quadratic order of P?. The effect of the Galilean 
invariance breaking appears only in higher order terms. Neglecting these higher order terms 
and forcing the Galilean symmetry, we obtain the result: 


1 ol 1 (k kë /3 
Tr$(P0, P?) 8rags 2r? \2 iP? + P?/2-— 2u + 2k? 


a 2 i k 


(D.2.6) 
T T ER 


Before closing this appendix, let us evaluate spatial integrations of T$, which are necessary 
for calculating the atom-dimer scattering. At first let us consider quantities in Eq.(D.3.3). That 
is, we consider the case ag > 0 and calculate two following quantities: 


lo r3 (k? 1), [& Onl ei), (D.2.7) 


where the notations are defined in Eq.(D.3.4). From Eq.(D.2.6), we find that 


k2 
EE S ee k k3/3 )-4 at ae Dag k 


= = = J= 
T(k2,0) Snag m? |2 3+2 87 7 Etaj)’ 
(D.2.8) 

4 1 k? / 6r? 1 = tan” PEST k/8n? 
Ps"! .= = = a (D29) 


"TREO (8k? + 3)? 32m, [= 4 3 3k? + 2 
as 


We should emphasize that e vanishes quadratically in the limit k — 0 but E does not 
vanish. We now employ the leading order approximation in the expansion (D.2.5). Therefore, 


we get 
dl S(p.2 12 d*t S-1 | 729 pS,-l 
[Enen four rare pi 


1 ii ðS! 
= —____ [k-— tan! | ka| 2E . (D.2.10 
Tð r í oaks | or l ) 


Using ð r$ = —(T$)?0 r7}, we can similarly obtain that 


k q3] K „TS, Px k gl Da 14a, toa cen 20,5 Le 1 (D.2.11) 
(27)? (27)? 


S,-1 S,-1 S,-1 
ff ę Ay A E t ar 


MAS 4 pa rst Var 
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These are the formula to be used in the first and second attempts at calculations of the atom- 
dimer scattering. 

Next, let us consider quantities which may become necessary for the more complicated 
momentum dependence of g?4. Let us evaluate 


4 dl Tee) J dl. aT (k2, 1) (D.2.12) 

(2r)? [k2 + P + 1/a2)”’ (2r)? [k? + P + 1/a2]” 7 
for n = 0,1,2. The case of n = 0 is already calculated. Again using the same approximation, 
we find that 


k asl Poth) 
I (27)? [k? +12 + 1/a3] 
a KyO E VaR) has)? + 1] 
T kO k2(8 Tt /2r%71) — (ka2)/[(kas)2? + 1] 


Only = (kas)? (kas)? —1 oe 
x r ops tan tarki — mar k ops 


k gl Fok) 7 (kag)? /[(kag)? si] 
J (27)? [k2 + + 1/02) — 4r2krS! (D.2.14) 


(D.2.13) 


T (kas)?/(2(kag)? +1 
[A2(OsP R /2rF) = (kab)/ (kas)? + 1)| 


) 
| Vas Fas? +7) 
(POTS /a ge) — (ha)/ (kas)? + VP 


A {22k , (has?) a | (eas)? 
ary? (kag)? +1 (kas)? + 1 
Ot oy oa CE ee 
are’ \ (kas)? A 
j Bl A,PS(k?,P) _ (kas)?/[(kas)? +1] 
(Or) [k++ 1a] ers 
E Ol ECO eer) 
[A2(Oor 1/20 E7') — (ka8)/ ((kas)? + 1)] 


ky) RoR: 
Ore or) 1) — (ka2,)/ ((kas)? 


2k CAS E (kas)? r 
are! \ (kas)? ais (kas)? +1 


yOu a, cease Ore 
4 arg? * (eas)? F Je G a Hii 


(D.2.15) 
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(D.2.16) 


2 


frat BED - (PATE /OTÈ D) ((kas)?/[(bas)? + 11? 
(27)? [ke +P +1/a3] p250 78-1 [a2(OuPg*/2P 8) — (ka3)/ (kas)? + 1)| 


(KO, * + 207") (2(kag)? + 1) 
1 


[k28 7/202") — (ka3)/ ((kas)? + 1)] 


- (3(kas)? + 1)k? aLr + (kas)? 2r% 


D.3 Numerical approximate solution of atom-dimer scat- 
tering 


In this section, let us calculate atom-dimer scattering within an approximation neglecting feed- 
back effect from higher-order vertices. 

In order to perform a further simple analysis, let us make a crude assumption to neglect 
the momentum dependence of the vertex function g@@. In the integrand of the right hand side 
of Eq. (4.2.25), we put pı = p} = (i/a%,0) and poy3 = ph}3 = 0. Then we find an ordinary 
differential equation 


ad 5 ad P (Slp) Te (pi + 0) 
ede =e i (2s (s+ RD Gs + Rl) [GF RPGS + ncn) a Re 


As we have mentioned, we drop the feedback terms by hand. Performing the /° integration, we 
find that 


TSP? = 1? + R,(I2), 2) 
AP + 1/02 + R(E)P 


8kg = a, f del 2 ale = I? + R), 1?) 
k (27)? i 2/1? + 1/az + R,(L)| 


Recall that we adopt the regulator R;,(l?) = (k? — 1?)6(k? — I’), and then its derivative takes 
the form kO,R,(l?) = 2k?6(k? — 1°). Then Eq.(D.3.2) becomes 


k ql TS (k? TS (k2, 77 
kosa = 208 f sk? = | ? Ed yal J 
(27) (k2? + 1/as) k? + 1/a3 


k 43 S(p.2 72 S(k2 2 
+f al, (ee " AOR), (D.3.3) 


) . (D.3.2) 


(27)? k241/a2)3 © 2(k? + 1/a})? 
where aie 
r 
Pe P) =r (iP? = k’, P), orik, P) = TI (iP? = k? P). (D.3.4) 
1 


These quantities are calculated in Appendix [D] Especially, see Eqs.(D.2.10)D.2.11). 
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Figure D.2: k-dependence of the vertex function g? in the first attempt. 


By solving the ordinary differential equation (D.3.3) using Eqs.(D.2.10}D.2.11), we find that 


Ge fay —0.224. (D.3.5) 
The k-dependence of g@4 in our approximation can be found in Fig/D.2| As a result, the 1PI 
contribution to the atom-dimer scattering T-matrix at zero energy is given by 

Tet Jag =~ —22.5. (D.3.6) 
Therefore, the atom-dimer scattering length aq becomes 


Qaq/ag ~ 1.47. (D.3.7) 


The correct value of the atom-dimer scattering length was derived by G. Skornyakov and K. 
Ter-Martirosyan using the three-body Schrodinger equation, and it is given by aaa = 1.18a5. 
Our result (D.3.7) is different from this correct value, and it would come from lack of the 
feedback from higher vertex functions. However, it still improve the result of the tree 
diagram contribution significantly. 
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Appendix E 


Channel decomposition of the 
dipole-dipole interaction 


In this appendix, we explicitly calculate the decomposition of the dipole-dipole interaction 
(DDI) into channels labeled by angular momenta (J, Jz, L, S). 


E.1 DDI in the spherical coordinate 


At first, we rewrite the DDI in Eq. (5.2.3) in the spherical tensor basis. Since it has the quantum 
number L = 2 and J = 0 it should be proportional to 


Vir (q) x wee m3;2,m 5 C maS mı Jal Sma) B' D? ET ey ilo pr fp: (E.1.1) 


m mı, m2 H1,H2 


Here we have used the spherical basis for the vectors q = k — k’ and for the spin vector S: 
d+1 = (de + iqy)/V2, qo = qz and S41 = +(S, +48,)/V2, So = Sz. Indeed, we shall now show 
that 


2 
Q q a m Qa 2,m 
(q ` Sg&)(q 7 S$) L = al ` g, z SiC a X a. „mı; biriak (Srila mala À Ci ui, ll lua 
m mM, m2 H1,H2 


(E.1.2) 


Proof. Using the formula and symmetries (C.3.3) of the Clebsch-Gordan coefficients, 


we get 


b+ +y = — 
(=) (=) EOE ge ORB ay = by, (6.13) 
aQ 
Here we should notice that since Cy" 1p, = (—)'t**Cim 1, we need not consider a = 1 in the 


sum in Eq.(E.1.3) for our purpose. Now we easily find that (C°? = 
lm:l,—m J2i+1 


m m m Í -m 1 -mı —p1 
= a Ca erie luo — (= ee oa "Oma, — pa Om -m > “Cs (=) j Onna ,—ma 9 p11 pr 
3 


m 


+(anti-symmetric 41 + u2). 


Since >, (—)”£mY-m = «+ y for the vectors, we get the result. 
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Therefore, we have got the formula 


y? 
Veta) = = OP SS O a a O >) ae E 


q’ m M1, m2 H1,H2 


Substitute q = k — k’ with k = k'. Notice that k, = kY,,(k) and for k' < k 


= aE = aÈ AY. ( Yi4,(k) + const.6(k — k’), (E.1.5) 


l 
1 
where H; = J — is the harmonic number. The coefficient of the spherical delta function is 
n 


n=1 
constant in terms of the label l, and in the limit k’/k — 1 the coefficient diverges. 
Proof. Assume that k’ < k and take the limit k’ > k later. We can now use the sum rule 


1 k! 4r Á ‘ 
—— = —— Ym (kK) Y¥*,(K’). 
lk — k'| DDES im(k)Yim(k') 


Therefore, we get 


1 kn t2 (4r)? i a $ A a 
|k — k'|? — ` ` ķku+l2+2 (2l aL 1)(2l2 + pml Yim ( )Yiomal ae ). 


lmi l2,m2 


The formula (C.1.11) and the orthogonal relation (C.3.4) give 
1 pe 10 10 4r v> (L 
kke ` z Tira r2 Ch0,20010,120 inl) Yim (Re!) 
l1,l2 


lm 


Let us perform the summation over l4,l2 in order to simplify the formula and to extract the 
singularity in the limit k’ > k. 
This limit should be treated with some care. The infinite series 5% DOn 1,(Ci?o1.0)” diverges, 


2l+1 


3,71: I order to circumvent this artificial divergence, we write this series 


. 10 
since C9149 & bn. 
as 


4n 
(Choo) = 4m 2 (Choto) — (Choo) | +47 De (E.1.6) 
2l +1 TEA +1 


li,l2 11 ,l2 


and the first term becomes —27H; and the second term is a constant independent of l. 
Let us justify the above result. First we now use the formula (C.1.11): 


aap ua)? = | aePa (2) P(e), (E.1.7) 


In order to make the infinite series in Eq.(E.1.6) finite, we again multiply the convergence factor 
r= k'/k <1 so that 


co 


2 ly+l2 / l0 r= ly lə z P(x) 
2+1 Dt (Giogo) = KOJ P(x (e) Sor Pa (x) P(x) = Perey, 


11 ,l2=0 1,=0 l2=0 
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We would like to take the limit r — 1, but that procedure gives the divergence. The origin of 
the divergence is independent of l, so we here evaluate the difference of them: 


oo a 


xı := lim (E.1.8) 


r>ljJa 1 — 2rgz+r? 2 


In order to evaluate xı, we consider its generating function 


= 1 fa 1-t 
of aS f 2 ( — i) = In(1 — t). (E.1.9) 
= 2J/,l-a2\V1—-2i2#+t 


Therefore, we find that x; = —1/l for l > 1, and then we obtain the result. 


Since k, = kY.,(k), we can rewrite Eq.(E as 


Vf (k-k) =y 2e 27 ALY, ( ie + const.6(k — k’) 


B'or ( 
sO =p > Gin, to Ona) ona) 


™1,™2 


x DO OP ua (Yiya (R) — Yip (R) Yipak) — Yayo (k). (E110) 


H1,H2 


The spherical delta function gives zero, and we can now take the limit k’ = k. The spin matrix 


elements are given as 
(Sim)or = (sa| Smlsa') = /s(s + 1) saim 
which gives 


Vein (ke — k') = —2n7?s(s + 1) X` HY (kY, (É') (E.1.11) 
ile 
x Du) Ji ` D „mı;l, m2 a seven On Im2 
m1,m2 
x 3 Cu. „u1;l, u2 (Yim (k) = Yim (K)) (Yin (k) T: Yip (k')) 
H1, H2 
= —2r7?s(s + 1) XC MY. (RY (k') (E.1.12) 
lls 
xŠ (-) ye i ` E aa EN Sal iia e D is 
m mı,m2 H1,H2 
TOi i f n7 f al 
x |X ———— Inl + 5O tuina O10100¥re (k) T Yiu (k )) = 2Y tp, (k)Yin(k') 


an i 


E.2 Particle-particle channel 


Apply the channel decomposition (5.2.9) ) to the expression (E.1.12): 
yee mL be Jj ors SS. S's 
Vis, L's! = ` CiM, SS, Chm SS > > Coa sai sp! (E.2.1) 
M,Sz M',S¢ aß a’ Bp! 


x fam f dO Yřu(k k)Y ue (kl )[Vire (ke —k’)-V, Se + k’)). 
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There are two terms which we should compute, however the second term can 0 represented in 
terms of the first term. a us denote the first term as V given in Eq.(E and the second 
term as V given in Eq.(E so that V = V + V. Then 


VATE B =y y JJ Yr SS- oss 
Visas! ~~ Cruss.C ae oS! CasB? sa’,s3! (E.2.2) 


M,Sz M’,S! aß a'b! 


1 dO, f de Yrm (ÈY rm (ÈVE (k + k’) 


= JJz cr J; y y SS- Cs: 
<a = 2, Cim,ss.C L'M’ O'S! Ci. „SB ORTAN (E.2.3) 


M,Sz M',S!, aß a’ B! 
a dO, f dw Yiu (k)Yy m (~k) V (k — k’) 
1425—59" Yodel Je 
a ee ae (E.2.4) 


JJaJ J, 1L 2s- SI ed J, 
Therefore, it suffices to calculate V because Viszis = [L+ (-) Wists 


Now we find that 


Ide Is ITs JL SS. Ons: 
Visas! =) ) 5 ) CiM ss. C e SIS! ` > Keen 8B sa’,sf (E.2.5) 


M,S, M’,S! ab al pt 


x fam f advil \Yr m (k’) 
x —Qry*s(s +1 J> Yi (AM, (F) 


lilz 


J m 2, m 
x X (—) S Cn sm 31,m2 so Imi Cy Im2 ; Cy 1,uislue 
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m1,m2 H1,H2 


x |) es 4n(2U + 1) ate CD 10 (Yu (k) + Yiu. (k’)) — Yin (k)Yi (É | , 


AA 
= —2ry7*s(s + 1) dM yes Ciii.ss.Com, SS! (E.2.6) 


M,Sz M',S% 


— 
z 2, SS. S'S! 
x S- = Ep Ci ET ` > Cia, abe eal sp! OP te ins 
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m1,m2 aß a'p' 


x > Ona. — Yom (k Yin (kY (E 


lz,41,H2 


K Ar (2l + 1) $C i 1m C010 Vir, (k) + Yoru (k) — 2Y aye, (k)Y ina (k') 


RA 


Here we should notice that the summation over l, is done at the underlined region in Eq.(E.2.6). 
Let us perform the summation over the spin indices in the braced part in Eq.(E.2.6). Using 
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the formula (C.3.12) containing the 9j-symbols, we can find that 


S 

X SSz S'S; sa 8B _ 2 } DEL] kk SS- 7 
Coa OET sa’, lmi CB 1m2 = IT 557 2k + 1C, Im2 OSS kr S 

apa’ B’ k 


kr 
Performing the summation over m1, M2 gives J C im2C Imm = 92,k0—m,x from Eq. [34 


mı, m2 


and then the over-braced part becomes 


(E.2.8) 


Wm RAW 
rH WH 
rH WH 


V525 + 1)2(28' + 1)C EE om 


Next, let us perform the angular integration in the underlined part of Eq.(E.2.6). The key is 


the formula 
Db TO AY apa 
fæ Yima Yımı Yomo _ / * M 2 Jos 3 C ets (E.2.9) 


An(2l3 +1) Hmm 


which is an immediate result of Eq.(C.1.11). Then the integration in the underlined part 
becomes 


3 Vit 


vo 
JNa 1m4 po 10,10 


vU, 


2i+1)(2’4+1 20’ +1)(2' +1 
x N )( ) gim CiooS LOM" Le + isda M Let ar) 


Ar(QL +1) Vah 4r(2l+ 1) 


3(21 + 1) d 3(2L/ + 
n(2L a 1) SATI ope Cito m(2l + Ct uaCL'0, 10: (E.2.10) 


: : : 2; UU, : 
Again the summation over 41, 42 gives -p p Cig aps C ya 1u = 92m, in the first term, and 


in the second term the summation over /,, 41, and u2 gives 


1 2 
3 ORE ae Ct tun CHM Ay = V5(2E + 1) Jha | z a (E.2.11) 


lz,41,H2 


from Eq.(C.3.8), then the underlined part becomes 


3 (2L 1) ny 
Jz Cio.10 (Or + OL) (OL HA) CL M2m Cin, 20 


34, PU +1, 1 2 
=N Coil Po 0v 5(21+ 1) Jehan | F 7 i: (E.2.12) 
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Substituting Eqs. (E.2.8]E.2.12) into Eq.(E.2.6), we get 


at SLE i JJ 
Visits = —2n7y"s(s + 1) ` H, ` ` Cee ae S'S! (E.2.13) 
l M,Sz M’,S! 
S s & 
x X (—) 528s + 1P (28 +1) fho m4 5’ 8 8 
2 11 
3 (2L' + 1) l 
Jä (ôr + 61,1) TAL 4 1) Clim 2mC Toa 
3 2L 1 1 2 
E E41 Tog aC, 10v 521+ 1) C501 L'M' l LEI i 


Again, we calculate the summation over the projection quantum numbers M, M’,S,, and S. 
There are two terms in the square bracket in Eq.(E.2.13), however their dependence on the 
projection quantum numbers are similar. The first term becomes (from the formula (C.3.12)) 


JJ IS, LM SS. 
` Criss. Ch, S'S} Crim’ 2m Si S,:2, -m 


M,M',Sz,S!, 
T8 
= X VL +1) (28 + QI + 1) (2k + CH o, mCi p J’ LS! 
kK k 2 2 


= Ņ VQL4DQS+ VAI + 1) (2k + 1) COn -mC ino 4 T L S > (E-214) 
k k 2 2 


Performing the summation over m as in Eq.(E.2.14), we find that 


2—m JJ: JS, LM SS 
` (—) ` Crass. Ch „S'S, Crim! 2mC'sisi.2,—m 


m M,M’,S~,S! 
J L S 

= zyz y y/ (2L +1)(2S +127 +1) J L Sp. (E.2.15) 
0 2 2 


This 97 symbol can be reduced into the 67 symbol: 


J L S (—)S+h 424d S Sga 
J V S= 
0 2 2 ores ese r i} 
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Therefore, substitution of Eq.(E.2.15) into Eq.(E.2.13) gives 


~ 17 ; 3 
Visti’ =—(-)°* PM 6776 1.057 sls +1)(28+1)S M; (E.2.16) 
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x (ox, + ôr) Cio mO L'0,L0 + (-)" 100%, 10 Cho, 10 { 
We have now obtained that 


~ i FF 3 
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As a consistency check, let us confirm that a constant independent of l can be added to H, 
without any effect. Indeed, an additive constant term in H is canceled in Eq.(E.2.17)). To see 
this, we should notice that 


1 1 2 ~ 1 1 2 
E Chaho] LL 1 l= 5- aa ‘Cin, Cr T LIL 1 \ 
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) Cio, ME 20 


JOFO] 
—]1 _\L-0/_\L+L'-2 
m tox! ) l ) 


v5 v5 


Using the parity selection rule, we find that Eq.(E.2.19)) becomes (— 247 CH 10CH0,10/9; which 
implies cancellation of the constant part in Hj. 


C inr (E.2.19) 


As another consistency check, we can easily confirm that the expression in Eq.(E.2.18) sat- 
isfied restrictions coming from symmetry discussed in sec[5.2] Of course, it manifestly satisfies 
the rotational symmetry, and the parity is respected due to the third line of Eq.(E.2.18). 
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E.3 Channel decomposition of the contact interaction 


In the case of the contact interaction Vee) = gô% 80? (r), we can easily find the channel 
decomposition. For the immediate reference of the formula, we derive it here. In the momentum 
representation, the potential becomes 


V giy (Ki — k2) = 962,55. (E.3.1) 
Then, the BCS channel decomposition is 
ied IS 2 is: bye st 
Viss a, D Cini Sz Cru SE >» Se rere: sp" (E.3.2) 
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: f dOr ij dw Yiu (B)Y nab) [Veo palk — k’) — Viswa lk + BY] 
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a 4ngôr oL’ o Coo, 5s, C 00,5” S7 Ox Oe sp’ Daa’ dap" 60,8'5 Ba’ 
Sz, 5, aB ap 


= 4rg(1 — (—)°™= °) Szy ðJ, 7,0 Jsôssðr0ð1 0: (E.3.3) 
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